Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,536,873 G→C 5.9% P146R (CCG→CGG)  narW ← nitrate reductase 2 (NRZ), delta subunit (assembly subunit)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,536,8730GC 5.9% 531.4 / 19.9 171P146R (CCG→CGG) narWnitrate reductase 2 (NRZ), delta subunit (assembly subunit)
Reads supporting (aligned to +/- strand):  ref base G (78/83);  new base C (3/7);  total (81/90)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

Only 100 of 176 total aligned reads displayed.

CAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCT  >  NC_000913/1536724‑1537020
                                                                                                                                                     |                                                                                                                                                   
cAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACg                                                                                                                                                     >  2:100413/1‑150 (MQ=255)
cAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACg                                                                                                                                                     >  2:2203585/1‑150 (MQ=255)
 aaTAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                                                                                                                                    <  1:1912886/150‑1 (MQ=255)
 aaTAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                                                                                                                                    >  1:1603088/1‑150 (MQ=255)
 aaTAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGATGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                                                                                                                                    >  2:706456/1‑150 (MQ=255)
  aTAAACTTCACCTGTTCCTCTTCCCAAACTGGATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTTACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCTAACAACGCGTACCACGGt                                                                                                                                                   <  1:399602/150‑1 (MQ=255)
  aTAAACTTCACCTCTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGt                                                                                                                                                   <  2:1762920/150‑1 (MQ=255)
   tAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGAGTACCACGGTg                                                                                                                                                  >  1:1687695/1‑150 (MQ=255)
    aaaCTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGc                                                                                                                                                 <  2:1408064/150‑1 (MQ=255)
           aCCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGAGACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgt                                                                                                                                          <  1:883344/150‑1 (MQ=255)
           aCCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgt                                                                                                                                          <  2:1627370/150‑1 (MQ=255)
               gTTCCTCTTACCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGATGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgttgtt                                                                                                                                      <  2:1508842/150‑1 (MQ=255)
                ttCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGttt                                                                                                                                     <  2:64369/150‑1 (MQ=255)
                ttCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGttt                                                                                                                                     >  1:57156/1‑150 (MQ=255)
                 tCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGAAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGtttt                                                                                                                                    >  1:1639616/1‑150 (MQ=255)
                  ccTCTTCCCACTCTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttga                                                                                                                                   <  2:1743877/150‑3 (MQ=255)
                  ccTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTa                                                                                                                                   >  1:830724/1‑150 (MQ=255)
                  ccTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTa                                                                                                                                   >  2:1894577/1‑150 (MQ=255)
                  ccTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttga                                                                                                                                   <  2:2154983/150‑3 (MQ=255)
                   ctctTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttgat                                                                                                                                  <  1:13180/150‑4 (MQ=38)
                   ctctTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTATTGTCACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttgat                                                                                                                                  <  1:1054942/150‑4 (MQ=38)
                     ctTCCCACACTGCATCAAGCGCCTGGCGGTAGTCATCGCGCTCTTCTCTGTTTACTTGTATTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGAGTACCACGGTGCCTCGCGTTGTTTTAAGc                                                                                                                                >  1:1474378/1‑150 (MQ=255)
                     ctTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttgatgc                                                                                                                                <  1:485446/150‑6 (MQ=38)
                       tCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGa                                                                                                                              >  2:757537/1‑150 (MQ=255)
                         ccACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGAcc                                                                                                                            >  1:1242909/1‑150 (MQ=255)
                              cTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCaa                                                                                                                       >  1:1856606/1‑150 (MQ=255)
                                gCATCAAGCGCCTGGCGGGTGTCATCGCGCTCATCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGc                                                                                                                     <  1:404518/150‑1 (MQ=255)
                                 cATCAAGCGCCTGGCGGGAGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACACCGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCa                                                                                                                    >  1:429167/1‑150 (MQ=255)
                                  aTCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTAACACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg                                                                                                                   >  1:298255/1‑150 (MQ=255)
                                  aTCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTCTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg                                                                                                                   <  1:170849/150‑1 (MQ=255)
                                   tCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAgg                                                                                                                  <  2:1529276/150‑1 (MQ=255)
                                   tCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAgg                                                                                                                  <  1:1736245/150‑1 (MQ=255)
                                      aGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCg                                                                                                               >  2:1358683/1‑150 (MQ=255)
                                        cgcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAg                                                                                                             >  1:1031404/1‑150 (MQ=255)
                                         gcCTGGCGTGTGACATCGCGATCTTCGCTTTTTAATTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCATCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAgg                                                                                                            <  2:377078/150‑1 (MQ=255)
                                         gcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAACCAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAgg                                                                                                            <  1:1460625/150‑1 (MQ=255)
                                         gcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAgg                                                                                                            <  1:1830034/150‑1 (MQ=255)
                                         gcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAgg                                                                                                            >  2:901143/1‑150 (MQ=255)
                                          ccTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGa                                                                                                           >  1:1797069/1‑150 (MQ=255)
                                                  gtgtCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCa                                                                                                   <  1:983841/150‑1 (MQ=255)
                                                    gtCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAAc                                                                                                 <  1:1312862/150‑1 (MQ=255)
                                                    gtCATCGCGCTCTTCGCTGATTACTTGTTTTGTGACACTGTCACTTGACAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAAc                                                                                                 <  1:1771486/150‑1 (MQ=255)
                                                    gtCATCGCGCGCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAAc                                                                                                 <  1:1787053/150‑1 (MQ=255)
                                                        tCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCCCGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTg                                                                                             >  2:1631713/1‑150 (MQ=255)
                                                          gcgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAg                                                                                           >  2:1697542/1‑150 (MQ=255)
                                                           cgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGc                                                                                          >  2:425889/1‑150 (MQ=255)
                                                             ctctTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCaa                                                                                        >  2:87568/1‑150 (MQ=255)
                                                                   gCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCcttc                                                                                  <  2:269585/150‑1 (MQ=255)
                                                                     tGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCcttctt                                                                                <  1:206574/150‑1 (MQ=255)
                                                                       tttACTTGGTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTc                                                                              <  2:144869/150‑1 (MQ=255)
                                                                        ttACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGACAAGCAGGGCGAGGATCGGCGCAACAATGAGCAATCCTTCTTTCg                                                                             >  2:1939596/1‑150 (MQ=255)
                                                                         tACTTGTTTTGTGACACTGTCACTCGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGCTCGGCGCAACATTGCGCAATCCGTCTTTCGc                                                                            >  2:631160/1‑150 (MQ=255)
                                                                          aCTTGTTTTGTGACACTTGCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGcc                                                                           <  1:1397840/150‑1 (MQ=255)
                                                                          aCTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGcc                                                                           <  2:1286907/150‑1 (MQ=255)
                                                                             tGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGa                                                                        <  1:706456/150‑1 (MQ=255)
                                                                                 ttGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGt                                                                    <  1:1358683/150‑1 (MQ=255)
                                                                                   gtgACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCa                                                                  <  1:673634/150‑1 (MQ=255)
                                                                                   gtgACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACAACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCCACATTGAGCAAACCTTCTTTCGCCGGATCGTcc                                                                  >  2:1101630/1‑149 (MQ=255)
                                                                                     gACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAgg                                                                >  1:2203809/1‑150 (MQ=255)
                                                                                       cacTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCa                                                              <  1:1100849/150‑1 (MQ=255)
                                                                                        acTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAg                                                             <  1:2143641/150‑1 (MQ=255)
                                                                                          tGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCCGGCAGTa                                                           >  2:2002742/1‑150 (MQ=255)
                                                                                          tGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGACGATCAGGcgaaa                                                           <  1:373263/150‑5 (MQ=255)
                                                                                            tCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACg                                                         <  1:757537/150‑1 (MQ=255)
                                                                                              aCTTGAAAGGGAGCTTCCCGCCAGTTGCGGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCt                                                       <  2:246725/150‑1 (MQ=255)
                                                                                                  gAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTaaa                                                   <  1:2024438/150‑1 (MQ=255)
                                                                                                  gAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTaaa                                                   <  1:1433171/150‑1 (MQ=255)
                                                                                                  gAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTaaa                                                   >  2:1065295/1‑150 (MQ=255)
                                                                                                   aaaGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAAt                                                  <  2:845313/150‑1 (MQ=255)
                                                                                                    aaGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTCCCTCGCGTTGTTTTAAGCGCCCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGCTCGTCAGGCAGTACGCTtatata                                                 >  2:469295/1‑150 (MQ=255)
                                                                                                      gggAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCCCGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACt                                               >  2:52353/1‑150 (MQ=255)
                                                                                                      gggAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACt                                               >  2:505398/1‑150 (MQ=255)
                                                                                                      gggAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACt                                               <  2:57156/150‑1 (MQ=255)
                                                                                                         aGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCa                                            <  2:2093820/150‑1 (MQ=255)
                                                                                                          gCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCCTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGCCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAg                                           >  2:1925738/1‑150 (MQ=255)
                                                                                                           cTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGCCCGCCAAGAAGTGCGAGGATCGGCGCATCATTGAGCAATACTTCTTTCGCCAGATCGTCAGGCACTACGCTTAAATACTCCAGa                                          >  2:1217111/1‑150 (MQ=255)
                                                                                                            ttCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGAt                                         >  1:508438/1‑150 (MQ=255)
                                                                                                             tCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATa                                        >  2:266839/1‑150 (MQ=255)
                                                                                                             tCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCGTAAATACTCCAGATa                                        >  2:628635/1‑150 (MQ=255)
                                                                                                                cGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAg                                     >  1:288980/1‑150 (MQ=255)
                                                                                                                 gCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATAGAGCAATCCTTCATTCGCCAGATCGTCAGCCACTACGCTTAAATAATCCAGATACAGc                                    >  1:263623/1‑150 (MQ=255)
                                                                                                                        gcagcaGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGAt                             <  1:867662/150‑1 (MQ=255)
                                                                                                                        gcagcaGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGAt                             <  2:1084290/150‑1 (MQ=255)
                                                                                                                        gcagcaGAGCATCAAACAACGCGGACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGAt                             <  1:1639230/150‑1 (MQ=255)
                                                                                                                           gcagAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAAt                          >  2:1028858/1‑150 (MQ=255)
                                                                                                                                 cttaaaaCAACGCGAGGCACCGTGGTACGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCa                    >  1:747002/5‑150 (MQ=38)
                                                                                                                                 cATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCa                    <  2:1153616/150‑1 (MQ=255)
                                                                                                                                  aTCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCATTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAg                   <  1:1894577/150‑1 (MQ=255)
                                                                                                                                  aTCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAg                   >  2:197802/1‑150 (MQ=255)
                                                                                                                                    cAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGtt                 >  1:1304883/1‑150 (MQ=255)
                                                                                                                                    cAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCATTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGtt                 >  1:1233952/1‑150 (MQ=255)
                                                                                                                                    aaaaCAACGCGAGGCACCGTGGTACGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTTAATACTCCAGATACAGCGGTAGATAATCGGGCAGtt                 >  2:1291716/2‑150 (MQ=255)
                                                                                                                                     aaaCAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTc                >  1:1014066/1‑150 (MQ=255)
                                                                                                                                       acaacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTccc              <  2:953832/150‑1 (MQ=255)
                                                                                                                                       acaacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTccc              <  2:2151417/150‑1 (MQ=255)
                                                                                                                                         aacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGa            <  2:1687695/150‑1 (MQ=255)
                                                                                                                                          acGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGAc           <  1:169827/150‑1 (MQ=255)
                                                                                                                                          acGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGAc           <  2:1639616/150‑1 (MQ=255)
                                                                                                                                           cgcgTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACa          >  1:1256585/1‑150 (MQ=255)
                                                                                                                                                 ccACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAg    >  2:2113125/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCT  >  NC_000913/1536724‑1537020

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: