Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,877,617 T→A 6.5% intergenic (‑152/+88) yidB ← / ← gyrB DUF937 family protein/DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,877,6170TA 6.5% 646.1 / 4.2 191intergenic (‑152/+88)yidB/gyrBDUF937 family protein/DNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base T (105/71);  new base A (5/10);  total (110/81)
Fisher's exact test for biased strand distribution p-value = 5.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

AATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAA  >  NC_000913/3877469‑3877766
                                                                                                                                                    |                                                                                                                                                     
aaTCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                                      >  1:511723/1‑150 (MQ=255)
      cTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAt                                                                                                                                                <  2:1620949/150‑1 (MQ=255)
      cTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAt                                                                                                                                                >  1:142351/1‑150 (MQ=255)
        gTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                                                                                                              <  2:387539/150‑1 (MQ=255)
         tAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCAt                                                                                                                                             <  1:823940/150‑1 (MQ=255)
         tAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCAt                                                                                                                                             >  1:1260296/1‑150 (MQ=255)
                ttACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAtt                                                                                                                                      >  2:2742705/1‑150 (MQ=255)
                ttACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAtt                                                                                                                                      >  2:2123904/1‑150 (MQ=255)
                   cGTGTTGACCCCAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCaa                                                                                                                                   <  2:2484047/150‑1 (MQ=255)
                   cGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCTGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCaa                                                                                                                                   <  1:133618/150‑1 (MQ=255)
                       ttGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACa                                                                                                                               >  1:472783/1‑150 (MQ=255)
                       ttGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACa                                                                                                                               >  1:278208/1‑150 (MQ=255)
                       ttGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGATTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAAAACa                                                                                                                               >  1:1806058/1‑150 (MQ=255)
                        tGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAt                                                                                                                              >  1:876632/1‑150 (MQ=255)
                        tGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAt                                                                                                                              <  1:1005439/150‑1 (MQ=255)
                         gACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAtt                                                                                                                             >  2:2336505/1‑150 (MQ=255)
                          aCCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTg                                                                                                                            <  2:252969/150‑1 (MQ=255)
                          aCCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTg                                                                                                                            >  1:2666449/1‑150 (MQ=255)
                          aCCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTg                                                                                                                            >  2:676402/1‑150 (MQ=255)
                          aCCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTg                                                                                                                            <  2:599667/150‑1 (MQ=255)
                            ccAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCa                                                                                                                          >  2:308957/1‑150 (MQ=255)
                             cAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCaa                                                                                                                         <  1:1195926/150‑1 (MQ=255)
                             cAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCaa                                                                                                                         <  2:102889/150‑1 (MQ=255)
                             cAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGAAAAATCAATACAGTGCaa                                                                                                                         >  1:43976/1‑150 (MQ=255)
                               aaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGa                                                                                                                       >  2:1496172/1‑150 (MQ=255)
                                    tattGATTTTCTGGATGATATTTGCTTATTGGAGAAATTATAAGGGCTTAGATGATCTCCTCCGTTTTCTTATAGCTTCTTGCCGGATGAGGCGTGAAAGCCTTATCCGTCCTACAACATCATGAAAATTCAATACATTGCAAGATTTTc                                                                                                                  <  1:1468420/150‑1 (MQ=255)
                                    tataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTc                                                                                                                  <  2:272739/150‑1 (MQ=255)
                                         aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAgt                                                                                                             <  2:1209040/150‑2 (MQ=255)
                                          tttGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGc                                                                                                            <  2:2267289/150‑1 (MQ=255)
                                          tttGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGc                                                                                                            >  1:1709549/1‑150 (MQ=255)
                                             gCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCGATACATTGCTCGATTTTCGGTGGCCTg                                                                                                         >  1:1936105/1‑150 (MQ=38)
                                             gCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTg                                                                                                         <  1:1718572/150‑1 (MQ=255)
                                             gCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTg                                                                                                         >  1:1631941/1‑150 (MQ=255)
                                                   tgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGc                                                                                                   <  2:659111/150‑1 (MQ=255)
                                                   tgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGc                                                                                                   <  1:2572107/150‑1 (MQ=255)
                                                    gatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGAAAAGCg                                                                                                  >  1:11292/1‑150 (MQ=255)
                                                     atAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGt                                                                                                 <  1:2298502/150‑1 (MQ=255)
                                                     atAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGt                                                                                                 >  1:1986238/1‑150 (MQ=255)
                                                       aaTTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGTTAAGCGTAg                                                                                               <  2:1265136/150‑1 (MQ=255)
                                                       aaTTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAg                                                                                               >  1:1765459/1‑150 (MQ=255)
                                                         ttGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgcg                                                                                             >  1:2017521/1‑150 (MQ=255)
                                                         ttGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgcg                                                                                             >  1:1602860/1‑150 (MQ=255)
                                                         ttGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgcg                                                                                             >  1:1778257/1‑150 (MQ=255)
                                                           gCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCa                                                                                           >  1:1553561/1‑150 (MQ=255)
                                                           gCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCa                                                                                           >  1:1553234/1‑150 (MQ=255)
                                                           gCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCa                                                                                           >  2:2377574/1‑150 (MQ=255)
                                                            cGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCAt                                                                                          >  1:1814607/1‑150 (MQ=255)
                                                              gATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCa                                                                                        <  2:2804655/150‑1 (MQ=255)
                                                               aTTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                       >  2:818871/1‑150 (MQ=255)
                                                                   cAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGGTAAGCGTAGCGCATCAGGCAc                                                                                   >  1:1072279/1‑150 (MQ=255)
                                                                    ctataagaaaaGGGCGGAGATCATCTCCGCCCTTTTCTTATAGCGTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACg                                                                                  >  1:614395/8‑150 (MQ=255)
                                                                    aGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACg                                                                                  >  2:21668/1‑150 (MQ=255)
                                                                     gAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGc                                                                                 >  2:1456679/1‑150 (MQ=255)
                                                                        aagaaaaGGGCGGAGATCATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTAATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATATTCGGAGGCCTGATAAGCGTAGCGCAGCAGGCACGCTCg                                                                              >  1:120718/4‑150 (MQ=25)
                                                                               gggCGGAGATCATCTCCGCCCTTTTCTTATAGCTTCTTTCCGGATTCGGCGTGAAAGCCATATACGGCCAACAAAAAAAATAAAATTCAATAAAATTCAAGATTTTAGTATGCCTGATAATCGTAGCGCATAAGGCACGCTCGCAAGGtt                                                                       >  1:731370/1‑150 (MQ=255)
                                                                               gggCGGAGATCATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATTAAAATTCAATACATTTCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATAAGGCACGCTCGCATGGtt                                                                       >  1:2258973/1‑150 (MQ=255)
                                                                                ggCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCATATCCCGCCAAAAAAATCATGAAAATTCAATACATTGCAAGAATTTAGTATGCCAGATAATCGTAGCGCAAAAGGCACCCTCGCATGGTTa                                                                      >  1:306234/1‑150 (MQ=25)
                                                                                       atgatCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCATGAAAAATCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCAGGGTTAGCGccct                                                               >  2:610798/1‑148 (MQ=255)
                                                                                           tctcCGCCTTTTTTTTATAGCTTGTGGCCGGTGGGGGCGTGACCGCTTTATCGGGCAAACAAAAAAATTAAAATTCAATACATTGCAATTTTTTGGTTGGCATTATATGGGTAGCGAATCAGGCACGCTCGCATGGTTAGCGCCATTaaa                                                           <  1:2767758/150‑1 (MQ=255)
                                                                                           tcgccGCCTTTTTCTTATAGCTTCTTGCCGGTTGCGGCTTGAACGCCTTATCCGCCAAAAAAAACAATTAAAATTAAATACATTGCAATATTTTCTTTGGCCTTATAAGCTTCGCGAATCAGGCACGCTCGCATGGTTAGCGCCATTaaa                                                           <  1:2049155/147‑1 (MQ=255)
                                                                                           tcgccGCCCTTTTCTTATTGCGTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCATACAAATTCATGAATTTTCAATACATTGCACGATTTTCGTAGGCCTTATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTaaa                                                           <  2:1329464/147‑1 (MQ=21)
                                                                                            ctcCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTAAATACATTTCAAGATTTTCGTAGGCCTGATAAGCGTATCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAat                                                          >  1:2889742/1‑150 (MQ=255)
                                                                                                gCCCTTTTCTTATCGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGCACAACAAACTCCTGAAAATTCAATACATTGCAATTTTTTGTTTGGCCTTATATGCTTTGCGAATCAGCCACGCTCGCATGGTTAGCGCCATTAAATATCg                                                      <  1:1514537/150‑1 (MQ=25)
                                                                                                      ttCTTATATCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATTAAAATTCAATACATTGCAATATTTTCGTAGGCCTTATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTc                                                <  2:1638917/150‑1 (MQ=255)
                                                                                                        cTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTcgc                                              <  2:1661305/150‑1 (MQ=255)
                                                                                                        cTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTcgc                                              <  2:958148/150‑1 (MQ=255)
                                                                                                           ataGCTTCTTGCCGGATTCGGCGTGAACGCCATATCCGGCCAACAAAAACATGAAAATTAAATACAAAGCAAGATTTTCGTAGGCCAGATAAGCGTAGCGCATCAGGCACGCTCGCAAGGTTAGCGCCAATAAATATAGATAATcgcccc                                           >  1:1938262/1‑148 (MQ=21)
                                                                                                               cttcttGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTAGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTc                                       >  1:2659175/1‑150 (MQ=255)
                                                                                                                  cttGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAgg                                    >  2:117882/1‑150 (MQ=255)
                                                                                                                    tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGc                                  >  2:752234/1‑150 (MQ=255)
                                                                                                                     gCCGGTTGCGGCGTGAACGCCTTATCCGCCAAACAAAATCATTAAAATTCAATACATTGCAAGATTTTCGTTGGCCTGATAAGCTTAGCTCATCTGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCg                                 <  1:2512213/150‑1 (MQ=37)
                                                                                                                      ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGt                                >  2:1144788/1‑150 (MQ=255)
                                                                                                                      ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGt                                >  2:475960/1‑150 (MQ=255)
                                                                                                                       cGGATGCGGCGTGAACGCCTCCTCCGGCCTACAAAATCATTAAAATACAATACTTTGCAAGATTTTCGTTGGCCTTATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGtt                               <  1:250276/150‑1 (MQ=255)
                                                                                                                            gcggcgTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTTGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCtt                          <  1:1675417/150‑1 (MQ=255)
                                                                                                                            gcggcgTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCtt                          <  2:491138/150‑1 (MQ=255)
                                                                                                                            gcggcgTGAACGCCTTATCCGGCCAACAAAATCATTAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCtt                          <  1:989410/150‑1 (MQ=255)
                                                                                                                            gcggcgTGAACGCCTTATCCGGCCAACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTTATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCtt                          <  1:1873152/150‑1 (MQ=255)
                                                                                                                             cggcgTGAACGCCTTATCCGGCCTACAAAATCCTGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTc                         <  2:2869539/150‑1 (MQ=255)
                                                                                                                             cggcgTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTAAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTc                         <  1:2930925/150‑1 (MQ=255)
                                                                                                                               gcgTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCaa                       >  1:626618/1‑150 (MQ=255)
                                                                                                                                cgTGAACGCCTTATCCGGCCTACAAACTCATGAAAATTCAATACAGTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAAt                      <  2:1859267/150‑1 (MQ=255)
                                                                                                                                 gTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATa                     <  1:2123904/150‑1 (MQ=255)
                                                                                                                                   caaCGCATTAACCGGCATACAAAATCATGAAAATTCAATTCATTGCAAGATTTTCTTAGGCCTTATATGGGTTGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATaaa                   <  1:37882/149‑1 (MQ=255)
                                                                                                                                    aaCGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAAc                  <  2:1553561/150‑1 (MQ=255)
                                                                                                                                    aaCGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAAc                  <  2:1553234/150‑1 (MQ=255)
                                                                                                                                       gCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCa               >  1:278318/1‑150 (MQ=255)
                                                                                                                                       gCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCAGGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCa               >  1:1362410/1‑150 (MQ=255)
                                                                                                                                       gCCTTATCCGGCCAACAAAATCATTAAAATTCAATACTTTGAATTATTTTCTTAGGCCTTATAAGCGTTGCGAATCTGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCa               <  1:1007815/150‑1 (MQ=255)
                                                                                                                                        ccTTATCCTGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCGCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAACCGCAc              >  2:2042094/1‑150 (MQ=255)
                                                                                                                                        ccTTATCCGGCCTACAAAATCATTAAAATTCAATACATTGCAATATTTTCGTAGGCCTGATAAGCGTAGCGCATCTGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCAc              <  1:727354/150‑1 (MQ=255)
                                                                                                                                        cattATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCAc              <  1:1928462/148‑1 (MQ=255)
                                                                                                                                          ttATCCGGCAAAAAAAATCATGAAAATTCAATAAATTGCTATATTTTCGTAGGCCTGATAAGCTTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATAGTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACgg            <  1:397416/150‑1 (MQ=38)
                                                                                                                                             tCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGgcgc         >  1:2932735/1‑150 (MQ=255)
                                                                                                                                             tCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGgcgc         >  2:2051042/1‑150 (MQ=255)
                                                                                                                                              ccGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGcgcg        <  2:1128051/150‑1 (MQ=255)
                                                                                                                                              ccGGCCTACAAAAGCATGAAAATTCAATACATTGCAAGATTTTCGTTGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGcgcg        <  2:534153/150‑1 (MQ=255)
                                                                                                                                               cGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCgg       >  1:797777/1‑150 (MQ=255)
                                                                                                                                                    tACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCaa  <  2:35732/150‑1 (MQ=255)
                                                                                                                                                    tACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCaa  <  2:1126997/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAA  >  NC_000913/3877469‑3877766

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: