Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,218,770 C→A 3.1% L301I (CTA→ATA)  metA → homoserine transsuccinylase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,218,7700C→A 3.1% 5.8 ‑1.2 0.0233L301I (CTA→ATA) metAhomoserine transsuccinylase
Reads supporting (aligned to +/- strand):  new base (2/6):  ref base (102/123):  total (104/129)
Fisher's exact test strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

GATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTT  >  NC_007779/4218671‑4218870
                                                                                                   |                                                                                                    
 aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTa                                                                                                    >  1:7921223/1‑101 (MQ=255)
 aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTa                                                                                                    >  1:6368514/1‑101 (MQ=255)
  tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTAc                                                                                                   <  1:1281922/101‑1 (MQ=255)
  tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTAc                                                                                                   <  1:11629299/101‑1 (MQ=255)
   ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTCCg                                                                                                  >  1:1707586/1‑101 (MQ=255)
   ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACg                                                                                                  <  1:8193883/101‑1 (MQ=255)
        aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGcaca                                                                                             <  1:10012619/101‑1 (MQ=255)
         aaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACAt                                                                                            <  1:2924426/101‑1 (MQ=255)
          aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATg                                                                                           <  1:4885528/101‑1 (MQ=255)
           tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGa                                                                                          >  1:5944160/1‑101 (MQ=255)
            ncacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGaa                                                                                         >  1:13066052/2‑101 (MQ=255)
            acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGaa                                                                                         <  1:12774254/101‑1 (MQ=255)
             cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAAt                                                                                        >  1:6564598/1‑101 (MQ=255)
             cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAAt                                                                                        <  1:8726257/101‑1 (MQ=255)
              acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                                                       >  1:3147330/1‑101 (MQ=255)
               ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATcc                                                                                      >  1:13691986/1‑101 (MQ=255)
                  cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAAc                                                                                   >  1:5584593/1‑101 (MQ=255)
                  cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAAc                                                                                   <  1:9691139/101‑1 (MQ=255)
                   gagcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACg                                                                                  <  1:9304668/101‑1 (MQ=255)
                   gagcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACg                                                                                  <  1:3825141/101‑1 (MQ=255)
                      cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg                                                                               >  1:6200778/1‑101 (MQ=255)
                       gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTGCCAGAGCACGCCATACGATCGACGGCACATGAATGCAACGCTgg                                                                              >  1:10579308/1‑101 (MQ=255)
                          cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAtt                                                                           <  1:10085310/101‑1 (MQ=255)
                          cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAtt                                                                           >  1:9448076/1‑101 (MQ=255)
                          cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAtt                                                                           >  1:7872662/1‑101 (MQ=255)
                            ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTaa                                                                         <  1:14345737/101‑1 (MQ=255)
                              cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtc                                                                       >  1:5335063/1‑101 (MQ=255)
                               gTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtct                                                                      >  1:10829821/1‑101 (MQ=255)
                               gTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtct                                                                      >  1:14463559/1‑101 (MQ=255)
                                tAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtctt                                                                     <  1:7248191/101‑1 (MQ=255)
                                 aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                                    >  1:14067099/1‑101 (MQ=255)
                                  gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttct                                                                   >  1:4057156/1‑101 (MQ=255)
                                    cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                                 >  1:1612828/1‑101 (MQ=255)
                                      cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                               <  1:5622389/101‑1 (MQ=255)
                                      cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATAAGGCCATACGATATACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                               >  1:3072665/1‑101 (MQ=255)
                                       ggTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                              <  1:2746129/101‑1 (MQ=255)
                                        gTAATTTACTGTTTCCCAACTGGCTCTACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATa                                                             <  1:10492758/101‑1 (MQ=255)
                                         tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                            <  1:6895217/101‑1 (MQ=255)
                                         tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                            <  1:10647115/101‑1 (MQ=255)
                                         tAATTTACTGTTTACCAACTGGCTCAACTAATACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                            >  1:13904026/1‑101 (MQ=255)
                                         gaaTTTAATGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATCCGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                            <  1:2989762/100‑1 (MQ=255)
                                           aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                          >  1:9177841/1‑101 (MQ=255)
                                           aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                          <  1:57575/101‑1 (MQ=255)
                                           aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                          >  1:13144787/1‑101 (MQ=255)
                                           aTTTACTGTTTACCAACTGGCTCAACTATTAAGTCTACCAGATCACTCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                          <  1:10377690/101‑1 (MQ=255)
                                            tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                         <  1:7548875/101‑1 (MQ=255)
                                              tACTGTTTACCAACTGGCTCATCTATTTCTTCAAGCAGAAGACGGCATACGAGATACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                                                       <  1:8237636/101‑1 (MQ=255)
                                              tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                                                       <  1:1758303/101‑1 (MQ=255)
                                              tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                                                       <  1:7669321/101‑1 (MQ=255)
                                                 tGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt                                                    >  1:11398354/1‑101 (MQ=255)
                                                  gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                                   >  1:5836292/1‑101 (MQ=255)
                                                  gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                                   <  1:4615262/101‑1 (MQ=255)
                                                   tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                                  >  1:2176723/1‑101 (MQ=255)
                                                   tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                                  <  1:4836321/101‑1 (MQ=255)
                                                    ttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                                 <  1:13297291/101‑1 (MQ=255)
                                                    ttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                                 >  1:4488037/1‑101 (MQ=255)
                                                      aCCAACTGGCTCAACTATTTTTTCAAGCAGAAGACGGCATACGAGCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                               <  1:11114879/101‑1 (MQ=255)
                                                      aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                               >  1:13341307/1‑101 (MQ=255)
                                                        cAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa                                             >  1:8936674/1‑101 (MQ=255)
                                                         aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                            <  1:9584059/101‑1 (MQ=255)
                                                            tGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                         <  1:12873992/101‑1 (MQ=255)
                                                            tGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                         >  1:3952092/1‑101 (MQ=255)
                                                            tGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                         <  1:844115/101‑1 (MQ=255)
                                                               cTCAACTATTACGTCAACCAGATCACGCCATACGAGCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                                      <  1:11436338/101‑1 (MQ=255)
                                                                tCAACTATTACGTCTACCAGATCACGCCATACGATATACGGCACATGAAGCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAc                                     >  1:9067238/1‑101 (MQ=255)
                                                                 cAACTATTTCGTCAAGCAGATGACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                                    <  1:13567683/101‑1 (MQ=255)
                                                                 cAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                                    >  1:9468745/1‑101 (MQ=255)
                                                                 cAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                                    <  1:6590552/101‑1 (MQ=255)
                                                                  aaCTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCt                                   >  1:2011164/1‑101 (MQ=255)
                                                                   aCTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCtt                                  <  1:8092222/101‑1 (MQ=255)
                                                                    cTATTACGTCTACCAGAACACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTa                                 <  1:14139100/101‑1 (MQ=255)
                                                                         tcGTCAAGCAGAAGACGCCATACGAGCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                            <  1:7724996/100‑1 (MQ=255)
                                                                         aCGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                            <  1:232276/101‑1 (MQ=255)
                                                                         aCGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                            <  1:6687062/101‑1 (MQ=255)
                                                                         aCGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                            <  1:1904517/101‑1 (MQ=255)
                                                                          tgTCTAGCAGAAGACGGCATACGAGCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                           <  1:685438/100‑1 (MQ=255)
                                                                          cGTCTACCAGATCACGCCATACGATCTACGGCACATGATTCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                           >  1:11918822/1‑101 (MQ=255)
                                                                           gTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                          >  1:9495333/1‑101 (MQ=255)
                                                                           gTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                          >  1:11004003/1‑101 (MQ=255)
                                                                                  aGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc                   >  1:11983559/1‑101 (MQ=255)
                                                                                   gATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCg                  >  1:2259426/1‑101 (MQ=255)
                                                                                       aCGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg              <  1:3892927/101‑1 (MQ=255)
                                                                                         gCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc            >  1:11271331/1‑101 (MQ=255)
                                                                                         gCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc            >  1:11252689/1‑101 (MQ=255)
                                                                                          ccATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt           >  1:773544/1‑101 (MQ=255)
                                                                                           cATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCtt          >  1:4120088/1‑101 (MQ=255)
                                                                                           cATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTAAGCACCTTACCTCAGGCACCTTCGGGTGCAtt          >  1:3163446/1‑101 (MQ=255)
                                                                                             tACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCtttt        <  1:6999694/101‑1 (MQ=255)
                                                                                             tACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCtttt        <  1:531311/101‑1 (MQ=255)
                                                                                              tcGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCttttt       <  1:12967062/100‑1 (MQ=255)
                                                                                              aCGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCttttt       >  1:12468946/1‑101 (MQ=255)
                                                                                               cGATCTACGGCACATGACTCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCtttttt      >  1:9113788/1‑101 (MQ=255)
                                                                                               cGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCtttttt      <  1:2590354/101‑1 (MQ=255)
                                                                                                gATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTa     >  1:5137744/1‑101 (MQ=255)
                                                                                                 aTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTAt    <  1:8353146/101‑1 (MQ=255)
                                                                                                  tCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTAtt   <  1:8598470/101‑1 (MQ=255)
                                                                                                   cTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTAttt  <  1:10261231/101‑1 (MQ=255)
                                                                                                   cTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTAttt  >  1:8292505/1‑101 (MQ=255)
                                                                                                   cTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTAttt  <  1:347485/101‑1 (MQ=255)
                                                                                                   cTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTAttt  <  1:11132618/101‑1 (MQ=255)
                                                                                                   |                                                                                                    
GATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTT  >  NC_007779/4218671‑4218870

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.