Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,536,868 A→G 5.8% Y148H (TAC→CAC) ‡ 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,8680AG 5.8% 523.6 / 14.6 173Y148H (TAC→CAC) ‡narWnitrate reductase 2 (NRZ), delta subunit (assembly subunit)
Reads supporting (aligned to +/- strand):  ref base A (82/81);  new base G (3/7);  total (85/88)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

CTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAAT  >  NC_000913/1536721‑1537014
                                                                                                                                                   |                                                                                                                                                  
cTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTAcc                                                                                                                                                  >  1:845313/1‑150 (MQ=255)
   cAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACg                                                                                                                                               >  2:100413/1‑150 (MQ=255)
    aaTAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                                                                                                                              <  1:1912886/150‑1 (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:2110069/150‑1 (MQ=255)
           ttCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcg                                                                                                                                       <  1:58893/150‑1 (MQ=255)
              aCCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGAGACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgt                                                                                                                                    <  1:883344/150‑1 (MQ=255)
               ccTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgtt                                                                                                                                   >  2:867662/1‑150 (MQ=255)
                  gTTCCTCTTACCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGATGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgttgtt                                                                                                                                <  2:1508842/150‑1 (MQ=255)
                   ttCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGttt                                                                                                                               <  2:64369/150‑1 (MQ=255)
                   ttCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGttt                                                                                                                               <  1:828772/150‑1 (MQ=255)
                   ttCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGttt                                                                                                                               >  1:57156/1‑150 (MQ=255)
                    tCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGAAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGtttt                                                                                                                              >  1:1639616/1‑150 (MQ=255)
                      ctctTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttgat                                                                                                                            <  1:13180/150‑4 (MQ=38)
                      ctctTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTATTGTCACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttgat                                                                                                                            <  1:1054942/150‑4 (MQ=38)
                        ctTCCCACACTGCATCAAGCGCCTGGCGGTAGTCATCGCGCTCTTCTCTGTTTACTTGTATTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGAGTACCACGGTGCCTCGCGTTGTTTTAAGc                                                                                                                          >  1:1474378/1‑150 (MQ=255)
                          tCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGa                                                                                                                        >  2:757537/1‑150 (MQ=255)
                           cccACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGAc                                                                                                                       >  2:1820437/1‑150 (MQ=255)
                           cccACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGAc                                                                                                                       >  2:295660/1‑150 (MQ=255)
                            ccACACTGCATCAAGCGCCTGGCGGGTGTCATCGCTCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGAcc                                                                                                                      >  2:1630405/1‑150 (MQ=255)
                                 cTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCaa                                                                                                                 >  1:1856606/1‑150 (MQ=255)
                                  tGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAg                                                                                                                <  1:2115663/150‑1 (MQ=255)
                                  tGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTCCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAg                                                                                                                >  1:428677/1‑150 (MQ=255)
                                    cATCAAGCGCCTGGCGGGAGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACACCGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCa                                                                                                              >  1:429167/1‑150 (MQ=255)
                                     aTCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTCTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg                                                                                                             <  1:170849/150‑1 (MQ=255)
                                      tCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAgg                                                                                                            >  1:2172957/1‑150 (MQ=255)
                                       cAAGCGCCTGGCGGGTGTCAGCGCGCTCTTCGCTGTTTCCTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAggg                                                                                                           <  1:1413529/150‑1 (MQ=255)
                                         aGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCGGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCg                                                                                                         >  1:246725/1‑150 (MQ=255)
                                         aGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCg                                                                                                         >  2:1358683/1‑150 (MQ=255)
                                           cgcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAg                                                                                                       >  1:1031404/1‑150 (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)
                                              cTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGAt                                                                                                    <  2:1603088/150‑1 (MQ=255)
                                                   gggTGGCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGgcg                                                                                               <  1:1309764/150‑1 (MQ=255)
                                                    ggTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGACAAGCAGGGCGAGGATCGgcgc                                                                                              >  1:175172/1‑150 (MQ=255)
                                                     gtgtCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCa                                                                                             <  2:1223497/150‑1 (MQ=255)
                                                     gtgtCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGGCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACCACGCGTACCTCGGTGCCTCGCGGTGTTTTAAGCGACCGCCAAGCAGTGCGAGGATCGGCGCa                                                                                             <  1:361743/150‑1 (MQ=255)
                                                       gtCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGAGCGGCGCAAc                                                                                           <  1:1778203/150‑1 (MQ=255)
                                                       gtCATCGCGCGCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAAc                                                                                           <  1:1787053/150‑1 (MQ=255)
                                                          aTCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACAtt                                                                                        >  1:1253773/1‑150 (MQ=255)
                                                              cgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGc                                                                                    >  2:425889/1‑150 (MQ=255)
                                                              cgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGc                                                                                    <  2:2004869/150‑1 (MQ=255)
                                                                ctctTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCaa                                                                                  >  2:87568/1‑150 (MQ=255)
                                                                     cGCTGTATACTTGTTTTGTGACCCTGTCACTTGAAAGGGCGCTTCCCGCCAGTTGCAGCAGAGCATCAAACACCGCCTACCACGGTGCCTCGCGTTGTTTTATGCGACCGCCCAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCctt                                                                             >  2:1071623/1‑150 (MQ=255)
                                                                        tGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCcttctt                                                                          <  1:206574/150‑1 (MQ=255)
                                                                        tGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCACGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCcttctt                                                                          <  1:1478103/150‑1 (MQ=255)
                                                                             aCTTGTTTTGTGACACTTGCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGcc                                                                     <  1:1397840/150‑1 (MQ=255)
                                                                             aCTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGcc                                                                     <  1:901143/150‑1 (MQ=255)
                                                                             aCTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGcc                                                                     <  1:1131566/150‑1 (MQ=255)
                                                                                    ttGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGt                                                              <  1:1358683/150‑1 (MQ=255)
                                                                                      gtgACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACAACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCCACATTGAGCAAACCTTCTTTCGCCGGATCGTcc                                                            >  2:1101630/1‑149 (MQ=255)
                                                                                        gACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGTCTGATCGTCAgg                                                          >  2:1524979/1‑150 (MQ=255)
                                                                                        gACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAgg                                                          <  1:707640/150‑1 (MQ=255)
                                                                                          cacTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCa                                                        <  1:1100849/150‑1 (MQ=255)
                                                                                            cTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCCCGGTGCCTCGCGTTGTTTTAAGCGACCGCAAAGCAGGGCGAGGATCGGAGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGt                                                      >  1:1611780/1‑150 (MQ=255)
                                                                                               tCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACg                                                   <  1:757537/150‑1 (MQ=255)
                                                                                                 aCTTGAAAGGGAGCTTCCCGCCAGTTGCGGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCt                                                 <  2:246725/150‑1 (MQ=255)
                                                                                                 aCTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCt                                                 <  2:1917725/150‑1 (MQ=255)
                                                                                                 aCTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCt                                                 <  1:295660/150‑1 (MQ=255)
                                                                                                 aCTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCt                                                 >  2:1445209/1‑150 (MQ=255)
                                                                                                   ttGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTa                                               <  1:485967/150‑1 (MQ=255)
                                                                                                    tGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTaa                                              >  2:1696854/1‑150 (MQ=255)
                                                                                                     gAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTaaa                                             >  2:1928230/1‑150 (MQ=255)
                                                                                                     gAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTaaa                                             <  2:389660/150‑1 (MQ=255)
                                                                                                      aaaGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAAt                                            >  2:119120/1‑150 (MQ=255)
                                                                                                       aaGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTCCCTCGCGTTGTTTTAAGCGCCCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGCTCGTCAGGCAGTACGCTtatata                                           >  2:469295/1‑150 (MQ=255)
                                                                                                         gggAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCCCGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACt                                         >  2:52353/1‑150 (MQ=255)
                                                                                                         gggAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACt                                         <  2:57156/150‑1 (MQ=255)
                                                                                                           gAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTcc                                       <  2:222196/150‑1 (MQ=255)
                                                                                                             gCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCTGATCGTCAGGCAGTACGCTTAAATACTCCAGa                                    <  2:1031404/150‑1 (MQ=255)
                                                                                                                tCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATa                                  <  2:1727786/150‑1 (MQ=255)
                                                                                                                tCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATa                                  >  2:266839/1‑150 (MQ=255)
                                                                                                                tCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCGTAAATACTCCAGATa                                  >  2:628635/1‑150 (MQ=255)
                                                                                                                   cGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAg                               >  1:288980/1‑150 (MQ=255)
                                                                                                                      cAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCgg                            <  1:688896/150‑1 (MQ=255)
                                                                                                                           gcagcaGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGAt                       <  1:867662/150‑1 (MQ=255)
                                                                                                                              gcagAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAAt                    >  2:1028858/1‑150 (MQ=255)
                                                                                                                              gcagAGCATCAAACAACGAGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACTCTTAAATACTCCAGATACAGCGGTAGATAAt                    >  2:2120557/1‑150 (MQ=255)
                                                                                                                                    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)
                                                                                                                                         aacaacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGAAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTcc         >  2:1396611/1‑150 (MQ=255)
                                                                                                                                          acaacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTccc        <  2:953832/150‑1 (MQ=255)
                                                                                                                                          acaacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTccc        <  2:2151417/150‑1 (MQ=255)
                                                                                                                                          acaacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTccc        <  2:1655301/150‑1 (MQ=255)
                                                                                                                                             acGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGAc     <  2:508438/150‑1 (MQ=255)
                                                                                                                                             acGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGAc     <  2:1966225/150‑1 (MQ=255)
                                                                                                                                             acGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGAc     <  2:1639616/150‑1 (MQ=255)
                                                                                                                                             acGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGAc     <  1:169827/150‑1 (MQ=255)
                                                                                                                                                cgTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGTGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAAt  >  2:2202217/1‑150 (MQ=255)
                                                                                                                                                cgTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAAt  >  2:636564/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAAT  >  NC_000913/1536721‑1537014

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: