Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,539,968 G→C 6.9% I200M (ATC→ATG narY ← nitrate reductase 2 (NRZ) subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,539,9680G→C 6.9% 45.6 ‑1.1 ‑2.1246I200M (ATC→ATGnarYnitrate reductase 2 (NRZ) subunit beta
Reads supporting (aligned to +/- strand):  new base (10/7):  ref base (106/123):  total (116/130)
Fisher's exact test strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.83e-03

Only 100 of 249 total aligned reads displayed.

TGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTAC  >  NC_007779/1539868‑1540067
                                                                                                    |                                                                                                   
tGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAg                                                                                                     >  1:5194211/1‑101 (MQ=255)
 gAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGa                                                                                                    >  1:6496252/1‑101 (MQ=255)
   agtttttttttttGTAGGGACACCCTTTTATGCACAACCGCCAGCCACGTCTTTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATa                                                                                                  <  1:9198226/95‑1 (MQ=255)
   aGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATa                                                                                                  <  1:1506551/101‑1 (MQ=255)
    gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg                                                                                                 <  1:1314175/101‑1 (MQ=255)
    gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg                                                                                                 <  1:1563382/101‑1 (MQ=255)
      aGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgcg                                                                                               <  1:7471357/101‑1 (MQ=255)
        aTTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGcgcc                                                                                             >  1:5103381/1‑101 (MQ=255)
          ttttttGTACGGACAACCGCTTATGCACAAACGCCGGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGcgccgc                                                                                           >  1:9333625/1‑101 (MQ=255)
           tttttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCt                                                                                          <  1:9754391/101‑1 (MQ=255)
           tttttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCt                                                                                          >  1:13847037/1‑101 (MQ=255)
            ttttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCtt                                                                                         <  1:4650319/101‑1 (MQ=255)
            ttttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGCCTTCTTCGCGTTTGTAGATAGCGCCGCtt                                                                                         <  1:4285866/101‑1 (MQ=255)
             tttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTg                                                                                        >  1:13542449/1‑101 (MQ=255)
             tttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTg                                                                                        >  1:9683719/1‑101 (MQ=255)
             cttGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGcgcccgatg                                                                                        <  1:8572484/100‑6 (MQ=255)
              ttGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGcgcccgatgg                                                                                       <  1:8981800/101‑7 (MQ=255)
               tGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTggg                                                                                      >  1:794569/1‑101 (MQ=255)
                gTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGc                                                                                     >  1:6529080/1‑101 (MQ=255)
                gTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGc                                                                                     <  1:12603235/101‑1 (MQ=255)
                  aCGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAg                                                                                   <  1:168492/101‑1 (MQ=255)
                   cGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAgg                                                                                  >  1:3861324/1‑101 (MQ=255)
                   cGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAgg                                                                                  <  1:14269778/101‑1 (MQ=255)
                   cGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAgg                                                                                  <  1:3741152/101‑1 (MQ=255)
                     gACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTc                                                                                <  1:8612613/101‑1 (MQ=255)
                      aCAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCg                                                                               <  1:955654/101‑1 (MQ=255)
                      aCAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCg                                                                               <  1:10903644/101‑1 (MQ=255)
                       cAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGc                                                                              >  1:6421710/1‑101 (MQ=255)
                          ccGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCAc                                                                           <  1:9420912/101‑1 (MQ=255)
                          ccGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCAc                                                                           >  1:7995659/1‑101 (MQ=255)
                          ccGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCAc                                                                           <  1:12259979/101‑1 (MQ=255)
                          ccGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCAc                                                                           >  1:11368976/1‑101 (MQ=255)
                           cGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACg                                                                          <  1:10498978/101‑1 (MQ=255)
                             cTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCa                                                                        <  1:1863562/101‑1 (MQ=255)
                             cTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCa                                                                        <  1:8108694/101‑1 (MQ=255)
                               tATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGc                                                                      >  1:11915170/1‑101 (MQ=255)
                                  gCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTgg                                                                   <  1:13791981/101‑1 (MQ=255)
                                   cacaAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTggg                                                                  >  1:6655268/1‑101 (MQ=255)
                                     caAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGAt                                                                >  1:4773415/1‑101 (MQ=255)
                                     caAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGAt                                                                <  1:10078746/101‑1 (MQ=255)
                                      aaaCGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGAtt                                                               >  1:7719428/1‑101 (MQ=255)
                                       aaCGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTg                                                              >  1:6112883/1‑101 (MQ=255)
                                        aCGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGa                                                             <  1:1675496/101‑1 (MQ=255)
                                          gccagccaCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAgg                                                           >  1:10737130/1‑101 (MQ=255)
                                          gccagccaCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAgg                                                           <  1:3639028/101‑1 (MQ=255)
                                            cagccaCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCa                                                         >  1:11802719/1‑101 (MQ=255)
                                             agccaCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAg                                                        >  1:3207043/1‑101 (MQ=255)
                                             agccaCGGCATTTATCCTGATCAACCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAg                                                        >  1:13477661/1‑101 (MQ=255)
                                                caCGGCATTTATCCTGATCAATCAGGACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAgtgt                                                     >  1:2308263/1‑101 (MQ=255)
                                                 aCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGtt                                                    >  1:4277330/1‑101 (MQ=255)
                                                    gCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCTCGTTTGTAGATAGCGCCGCGTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGc                                                 <  1:2844985/101‑1 (MQ=255)
                                                       tttATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGcaca                                              <  1:4701806/101‑1 (MQ=255)
                                                        ttATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAg                                             <  1:12210802/101‑1 (MQ=255)
                                                        cggtCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCAcagat                                             <  1:10425116/98‑4 (MQ=255)
                                                         tATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGACGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAgg                                            >  1:697239/1‑101 (MQ=255)
                                                          aTCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGgc                                           <  1:12737125/101‑1 (MQ=255)
                                                          aTCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGgc                                           <  1:39744/101‑1 (MQ=255)
                                                           tCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCACAGATgcg                                          <  1:6045846/101‑1 (MQ=255)
                                                           tCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGgcg                                          <  1:5415649/101‑1 (MQ=255)
                                                            ccTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGgcgc                                         >  1:3496291/1‑101 (MQ=255)
                                                             cTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGcgcg                                        <  1:7537619/101‑1 (MQ=255)
                                                              tggtcgaccaGCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCACAGATGCGAgg                                       >  1:4054613/9‑101 (MQ=255)
                                                              tggtcgaccaGCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCACAGATGCGAgg                                       >  1:2538130/9‑101 (MQ=255)
                                                                atcaatcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCa                                     <  1:5800973/101‑1 (MQ=255)
                                                                  cgaccaGCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCACAGATGCGAGGCAgg                                   >  1:7006802/5‑101 (MQ=255)
                                                                   aatcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGt                                  <  1:14194838/101‑1 (MQ=255)
                                                                   aatcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGt                                  <  1:12340212/101‑1 (MQ=255)
                                                                     tcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTAc                                <  1:9086068/101‑1 (MQ=255)
                                                                      caGCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCACAGATGCGAGGCAGGTaga                               <  1:10095139/101‑3 (MQ=255)
                                                                        gCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTAcatc                             >  1:11688216/1‑101 (MQ=255)
                                                                        gCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCACAGATGCGAGGCAGGTagaac                             <  1:2594691/101‑5 (MQ=255)
                                                                         cacaATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCACAGATGCGAGGCAGGTagaaca                            >  1:12965243/1‑96 (MQ=255)
                                                                          acaATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTAcatcat                           <  1:5594412/101‑1 (MQ=255)
                                                                          acaATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTAcatcat                           <  1:1368380/101‑1 (MQ=255)
                                                                            aaTGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGa                         <  1:3438691/101‑1 (MQ=255)
                                                                            aaTGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGa                         <  1:2696578/101‑1 (MQ=255)
                                                                             aTGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGaa                        >  1:593867/1‑101 (MQ=255)
                                                                             aTGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGaa                        <  1:13267301/101‑1 (MQ=255)
                                                                              tGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAg                       <  1:9939154/101‑1 (MQ=255)
                                                                               gCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAgg                      >  1:10663180/1‑101 (MQ=255)
                                                                                ccGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAggc                     >  1:6193553/1‑100 (MQ=255)
                                                                                 cGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGgtg                    <  1:9825024/101‑1 (MQ=255)
                                                                                  gTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGgtgt                   <  1:1340170/101‑1 (MQ=255)
                                                                                  gTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGgtgt                   >  1:3519289/1‑101 (MQ=255)
                                                                                  gTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGgtgt                   <  1:4086596/101‑1 (MQ=255)
                                                                                    cttcttCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGttt                 >  1:12712623/1‑101 (MQ=255)
                                                                                     ttcttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGtttt                <  1:641176/101‑1 (MQ=255)
                                                                                     ttcttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGtttt                >  1:2877391/1‑101 (MQ=255)
                                                                                        ttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCaa             >  1:10060243/1‑101 (MQ=255)
                                                                                        ttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCaa             <  1:8203621/101‑1 (MQ=255)
                                                                                        ttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGCGTTTTCaa             >  1:6819777/1‑101 (MQ=255)
                                                                                         tcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCaaa            <  1:4146347/101‑1 (MQ=255)
                                                                                            cgTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTg         <  1:8055793/101‑1 (MQ=255)
                                                                                            cgTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTg         <  1:7853291/101‑1 (MQ=255)
                                                                                            cgTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTg         <  1:7386867/101‑1 (MQ=255)
                                                                                             gTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGc        <  1:5671068/101‑1 (MQ=255)
                                                                                             gTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGc        <  1:2513094/101‑1 (MQ=255)
                                                                                              tttGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGcc       >  1:6443512/1‑101 (MQ=255)
                                                                                                 gTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGt    <  1:13991987/101‑1 (MQ=255)
                                                                                                  tAGATAGCGCCGCTTGGGCAGGTAGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTa   <  1:13092308/101‑1 (MQ=255)
                                                                                                    |                                                                                                   
TGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTAC  >  NC_007779/1539868‑1540067

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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.