Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,539,988 T→A 3.8% T194S (ACC→TCC)  narY ← nitrate reductase 2 (NRZ) subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,539,9880T→A 3.8% 18.3 ‑0.7 ‑0.1235T194S (ACC→TCC) narYnitrate reductase 2 (NRZ) subunit beta
Reads supporting (aligned to +/- strand):  new base (6/3):  ref base (114/112):  total (120/115)
Fisher's exact test strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.52e-01

Only 100 of 242 total aligned reads displayed.

GGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCG  >  NC_007779/1539888‑1540088
                                                                                                    |                                                                                                    
 gACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTc                                                                                                     <  1:6373652/101‑1 (MQ=255)
     aCCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCa                                                                                                 <  1:6002486/101‑1 (MQ=255)
      ccGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCAc                                                                                                >  1:1798132/1‑101 (MQ=255)
        gCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGc                                                                                              >  1:4930966/1‑101 (MQ=255)
        gCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGc                                                                                              >  1:7943514/1‑101 (MQ=255)
         cTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCa                                                                                             >  1:12189963/1‑101 (MQ=255)
         cTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCa                                                                                             <  1:11611288/101‑1 (MQ=255)
          ttATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAg                                                                                            >  1:5083546/1‑101 (MQ=255)
          ttATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAg                                                                                            >  1:7029797/1‑101 (MQ=255)
            aTGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCt                                                                                          <  1:4813748/101‑1 (MQ=255)
            aTGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCt                                                                                          <  1:11736948/101‑1 (MQ=255)
             tGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTg                                                                                         <  1:175243/101‑1 (MQ=255)
                 caAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGAt                                                                                     <  1:1454804/101‑1 (MQ=255)
                  aaaCGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGAtt                                                                                    >  1:2710152/1‑101 (MQ=255)
                   aaCGTCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTg                                                                                   >  1:13450433/1‑101 (MQ=255)
                   aaCGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTg                                                                                   >  1:6586137/1‑101 (MQ=255)
                      gccagccaCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAgg                                                                                >  1:12767125/1‑101 (MQ=255)
                           ccaCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAgtg                                                                           <  1:1688829/101‑1 (MQ=255)
                            caCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAgtgt                                                                          >  1:2009337/1‑101 (MQ=255)
                            caCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAgtgt                                                                          <  1:5904604/101‑1 (MQ=255)
                             aCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGtt                                                                         >  1:12426703/1‑101 (MQ=255)
                              cGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTc                                                                        <  1:12704321/101‑1 (MQ=255)
                               ggCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCg                                                                       <  1:2646529/101‑1 (MQ=255)
                               ggCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCg                                                                       <  1:12914192/101‑1 (MQ=255)
                                    ttATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAg                                                                  >  1:14100777/1‑101 (MQ=255)
                                    ttATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAg                                                                  >  1:10139357/1‑101 (MQ=255)
                                    ttATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAg                                                                  <  1:8127075/101‑1 (MQ=255)
                                    ttATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAg                                                                  <  1:10873123/101‑1 (MQ=255)
                                     tATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAgg                                                                 >  1:1839040/1‑101 (MQ=255)
                                     tATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAgg                                                                 >  1:5647913/1‑101 (MQ=255)
                                     ggtCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCAcagatg                                                                 >  1:10397538/3‑97 (MQ=255)
                                     ggtCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCAcagatg                                                                 >  1:5833734/3‑97 (MQ=255)
                                      attgttgtcaatcaGCACAATGCCGTCTTCCGAGCGTTTGTAGATAGCGCCGCTTGGTCCCGGCGCCGCGCCGCTGGGATTTAGGCAGTGTTCGCACCGgc                                                                <  1:6563747/94‑1 (MQ=255)
                                        ccTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGgcgc                                                              >  1:7446666/1‑101 (MQ=255)
                                        ccTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGgcgc                                                              <  1:3542877/101‑1 (MQ=255)
                                         cTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGcgcg                                                             <  1:13827113/101‑1 (MQ=255)
                                          tGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCgg                                                            >  1:4399707/1‑101 (MQ=255)
                                              cgatcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAgg                                                        <  1:12091224/99‑1 (MQ=255)
                                                atcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTa                                                      <  1:14177417/101‑1 (MQ=255)
                                                atcaGCACAATGCAGTCTTCTTCGCATTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTa                                                      <  1:8845816/101‑1 (MQ=255)
                                                accaGCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCACAGATGCGAGGCAGGTa                                                      >  1:6054340/3‑101 (MQ=255)
                                                 ccaGCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCACAGATGCGAGGCAGGTag                                                     <  1:9668545/100‑2 (MQ=255)
                                                  caGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTAca                                                    <  1:8627963/101‑1 (MQ=255)
                                                  caGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTAca                                                    >  1:4553588/1‑101 (MQ=255)
                                                  caGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTAca                                                    <  1:9090275/101‑1 (MQ=255)
                                                   aGCACAATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGGTCAAGCAGTGCTCACAGATGCGAGGCAGGTagaa                                                   >  1:9578467/1‑98 (MQ=255)
                                                     cacaATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTAcatca                                                 <  1:13780431/101‑1 (MQ=255)
                                                     cacaATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTAcatca                                                 <  1:6747214/101‑1 (MQ=255)
                                                     cacaATGCCGTCTTCTGTGCGTTTGTACATCGCGCCCGATGGGCAGGAAGACACACAGGTCGGGTTCAAGCAGTGCTCACAGATGCGAGGCAGGTagaaca                                                 >  1:5763854/1‑96 (MQ=255)
                                                      acaATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTAcatcat                                                >  1:13330479/1‑101 (MQ=255)
                                                        aaTGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGa                                              <  1:10072406/101‑1 (MQ=255)
                                                          tGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAg                                            <  1:10430436/101‑1 (MQ=255)
                                                           gCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGGACATCATGAAgg                                           >  1:9939854/1‑101 (MQ=255)
                                                             cGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGgtg                                         <  1:126731/101‑1 (MQ=255)
                                                             cGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGgtg                                         <  1:2391872/101‑1 (MQ=255)
                                                              gTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGgtgt                                        <  1:8465321/101‑1 (MQ=255)
                                                               tcttctTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGtt                                       <  1:8288919/101‑1 (MQ=255)
                                                                cttcttCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGttt                                      <  1:6187733/101‑1 (MQ=255)
                                                                  tcttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTc                                    >  1:13293416/1‑101 (MQ=255)
                                                                  tcttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTc                                    <  1:10574353/101‑1 (MQ=255)
                                                                   cttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCa                                   >  1:13947317/1‑101 (MQ=255)
                                                                   cttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCa                                   >  1:11407525/1‑101 (MQ=255)
                                                                    ttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCaa                                  >  1:13662721/1‑101 (MQ=255)
                                                                    ttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCaa                                  <  1:13248015/101‑1 (MQ=255)
                                                                    ttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCaa                                  >  1:4369245/1‑101 (MQ=255)
                                                                     tcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCaaa                                 <  1:3357558/101‑1 (MQ=255)
                                                                      cgcgTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAAc                                <  1:2800199/101‑1 (MQ=255)
                                                                      cgcgTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAAc                                <  1:7371297/101‑1 (MQ=255)
                                                                       gcgTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACt                               <  1:5448577/101‑1 (MQ=255)
                                                                       gcgTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACt                               <  1:12871548/101‑1 (MQ=255)
                                                                       gcgTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGGACATCATGAAGGTGTTTTCAAACt                               >  1:12047930/1‑101 (MQ=255)
                                                                       gcgTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGGACATCATGAAGGTGTTTTCAAACt                               >  1:9798908/1‑101 (MQ=255)
                                                                        cgTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTg                              <  1:3458025/101‑1 (MQ=255)
                                                                         gTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGc                             <  1:7771315/101‑1 (MQ=255)
                                                                         gTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGc                             >  1:7880382/1‑101 (MQ=255)
                                                                          tttGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGcc                            >  1:4116349/1‑101 (MQ=255)
                                                                           ttGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGccc                           >  1:6014126/1‑101 (MQ=255)
                                                                            tGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCg                          <  1:6890935/101‑1 (MQ=255)
                                                                            tGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCg                          >  1:7488468/1‑101 (MQ=255)
                                                                                gATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACa                      <  1:4501930/101‑1 (MQ=255)
                                                                                gATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACa                      <  1:1948170/101‑1 (MQ=255)
                                                                                 aTAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACAt                     >  1:652337/1‑101 (MQ=255)
                                                                                   aGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACAtct                   >  1:11100025/1‑101 (MQ=255)
                                                                                      gccgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCt                <  1:994850/101‑1 (MQ=255)
                                                                                      gccgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCt                >  1:9256686/1‑101 (MQ=255)
                                                                                       ccgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCtt               >  1:824501/1‑101 (MQ=255)
                                                                                       ccgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCtt               <  1:9004847/101‑1 (MQ=255)
                                                                                       ccgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGGACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCtt               >  1:10882319/1‑101 (MQ=255)
                                                                                        cgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCttt              >  1:6822537/1‑101 (MQ=255)
                                                                                         gcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCtttt             <  1:2332349/101‑1 (MQ=255)
                                                                                         gcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCtttt             <  1:3917040/101‑1 (MQ=255)
                                                                                            tGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCa          <  1:3187903/101‑1 (MQ=255)
                                                                                             gggCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCAt         >  1:12545602/1‑101 (MQ=255)
                                                                                              ggCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATg        >  1:5955303/1‑101 (MQ=255)
                                                                                               gCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATgg       >  1:3050486/1‑101 (MQ=255)
                                                                                               gCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATgg       >  1:12575204/1‑101 (MQ=255)
                                                                                                   gTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTAGATCTCCTTTTGCATGGCCTc   <  1:3351981/101‑1 (MQ=255)
                                                                                                   gTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTc   <  1:7585844/101‑1 (MQ=255)
                                                                                                   gTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTc   <  1:12897351/101‑1 (MQ=255)
                                                                                                    tCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCg  <  1:10201356/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCG  >  NC_007779/1539888‑1540088

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