Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 989,945 A→G 21.4% S138P (TCA→CCA)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,9450AG21.4% 80.8 / 2.8 48S138P (TCA→CCA) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (14/21);  new base G (3/8);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCA  >  NC_000913/989711‑990175
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
ccgtCGGCGCGAGGGAAAGCCGACGCGCCAGTTCGTAGTTGCTTGTGCCCCGAAACCAGGATTCCTGTTGCGGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAACGCGGCGACGGGTGCCAATATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                       <  2:1051987/236‑1 (MQ=255)
ccacgggCGGGAGGGAAAGCCGACGCGCCAGAACGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTGGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGGCGACTAATTTGCGGTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                       <  2:336952/234‑1 (MQ=255)
ccATCGGCGTGACGGGAAGCCGGCGCGCCAGAGCGCAGTTGCTGGTGCTCACAAACCAGGATTCCTGTCGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                       <  2:168200/239‑1 (MQ=255)
     ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGGCGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACCCCCGGAAGACGGCTAATTTGCCTTTCCGCGGGGTGGGGGGTTGCGCCACGTCGGCCCCCGGGGagcgataacg                                                                                                                                                                                                                                  >  1:438034/1‑239 (MQ=255)
     ggCGTGAAGGAAAGCCGACTCGCCAGAACGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                  <  1:616190/239‑1 (MQ=255)
     ggCGGGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCCAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTAAGACAGCAAGGCGCCGAACTTGTTTTTCATAGCGTGGAGGGGTTGCGCCACGGCGCGCGGCGGTGggcgctgacg                                                                                                                                                                                                                                  >  2:1176973/1‑233 (MQ=255)
       cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcga                                                                                                                                                                                                                                <  1:976188/238‑1 (MQ=255)
        gggAGGGAAAGACGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATg                                                                                                                                                                                                                              <  1:673649/237‑1 (MQ=255)
          gAGGGAAAGCCGACGCGACAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAAGCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGtt                                                                                                                                                                                                                            <  1:496940/239‑2 (MQ=255)
           aGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGAGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTa                                                                                                                                                                                                                            <  1:75419/238‑1 (MQ=255)
                    cGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTACAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCAc                                                                                                                                                                                                                 <  2:831923/238‑1 (MQ=255)
                             agagCGTAGTTGCGGGTGCCCACAAACCAGGATTCCTGTGGCGAGCGCTTAACATTGGTTTCTTGTACTTCGGGAGAAAAACGGGGACGGGTGCCAAAATCCATAAGAGGGAAGCGCGACATATAAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGTTTGCGAAACGTCGGCCGGCGGGGAGCGATAGCGATGTCCCATTTCACTGATAACCg                                                                                                                                                                                                        <  2:95167/238‑1 (MQ=255)
                                           ggTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAACCGAGAAGGCGACTAATTTGCTTTCCAGCGGGTCGGGGGGTTGCGTGACGGCGGCCCGCCGTGGTCGGTAGGGGCGCACCATTTTCATGATAAACGCCAGCACGGCAACtt                                                                                                                                                                                         >  2:509553/1‑239 (MQ=255)
                                                      aCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgat                                                                                                                                                                                                                                      <  1:1129982/186‑1 (MQ=255)
                                                           gATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAACCCCAGAAAACGACTAATTTGCTTTTCAGCGCGGCCGGGGCTTGCGCAACGTCGGCCCGCGGTAAGCTGTAGCGGTGTACCCATTCCACGATGACCTCCTCCAACAGGACTCACCAGAGGAGGACtt                                                                                                                                                                          >  2:406889/1‑238 (MQ=255)
                                                                cTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGGCGAGGGCTTGCGCAACGTCGGCCCGCGGTGAGCGATAGCGATGGACCATTTCACTGATAACCCCCAGCAAAGGAACTTCCCAGAGGGAGACTTCACGc                                                                                                                                                                     >  1:934156/1‑238 (MQ=255)
                                                                  gTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGCCTAATTTTCTTTTCAGCGTGTCGGGGGGTTGCGCAAAGTCGGCCTGGGGGGGGCGATAGCGGTGTACCATTTCACTGATAACCGCCACCAAAGGAACTTCCCAAGGGATGACTTCCCGcccc                                                                                                                                                                  >  1:90289/1‑237 (MQ=255)
                                                                     gCAGACGCTAAACGATGGTTTCTGGTACTGCGCGAGAAAACCGGGGACGGGGGCAAAAATCCATCAGAGGGAGGCGCGACATATCAGGACGGGGGGTAAACGCGGAGAGGTCGTCTAATTTGCTTTCCGGCGTGTCGAGGGCTTGCCAACCGCCGGCCGGAGGGGGGAGATAGCGAGGTACATTTTCACGGAGAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACggg                                                                                                                                                               <  2:438034/239‑1 (MQ=255)
                                                                       aGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCTGTGagcgata                                                                                                                                                                                                                                     <  1:207612/170‑1 (MQ=255)
                                                                       aGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCTGTGagcgata                                                                                                                                                                                                                                     >  2:207612/1‑170 (MQ=255)
                                                                             ttAAGGTTGGTTTCTGGTACTAGGCGAAAAAAAGGGGGAGGGTGGCAAAAACCAATCAGAGGGAGGGGCTAAATAACAAGACGGGGGGTGACCCCGTAGAAGTCGACTATGTTGATTTCCGGGGGGCGGGGGGTTGGCGCAACGTCGGCCTGCGGTGGAGGTTAGCGTGGTACCATTTCACTAATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTa                                                                                                                                                        <  2:341389/238‑1 (MQ=255)
                                                                             taaaggATGGTTTCTTGACATCGGCGAGAAAAACGGGGCGGGGTGCAAAAATCCAACAGGTGGAAGAGCGAAATATCAAGACCGGGGGTTAACGCGGAAAAGGCGACTAATTTGTTTCCCAGCGTGTGGAGGGTTTGCGAAACGCCGGCCGGCGGTGGGCGATAGCGAGGTACCATTTCCCGAAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTa                                                                                                                                                        <  2:383272/233‑1 (MQ=255)
                                                                                     ggTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGGCGACTCAATTTCTTTCCAGAGGGGCGGGGGGCTGCGCAAAGCCGGCCTGCGGGGAGCGGAAGCGAGGTAACAATTTTCTCATAAACCCCAGCAAAGGGACTTCCCACAGGAGGGATTTACTCCACGGGCCCCTTTAAAGaaaa                                                                                                                                                >  2:981972/1‑236 (MQ=255)
                                                                                       tttCTTGTCCTCCGGAAAAAAACCGGCGCGGGGGGCAAAAACCCATCAGATGGAAGCGGGAAATACAACGACCGGCGTTTACCGCCGAGAAGCCGACTAATTTGCTTCCCAGGGTGTCGAGGGTTTGCGAACCGCCGGCCTGCGGTGGGCGTTAGCGATTTACCATTTCCCGGAAACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATc                                                                                                                                              <  2:1059599/238‑1 (MQ=255)
                                                                                         tgttgcacggcgcgAAAAAACCGGCGAGGGGTGCAAAAACCCATAAGATGGAAGCGCGACATACCAAGACGGGGGTTTAACGCCGAGAAGGCGACTAATTTGCTTCCCGGCGTGTCGGGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACGGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAg                                                                                                                                           <  2:1030287/229‑1 (MQ=255)
                                                                                               aCTTCGCGAGAAAAACGGCGACGGGTCCCAAAACCCATCAGATGGAAGCGCGACATATCAAGCCCGGCGTTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTGGCGAAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGc                                                                                                                                      <  1:838957/238‑1 (MQ=255)
                                                                                                   cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGGCGACTAATTTGCTTTTCAGCGGGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGAGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGATTGCCATTa                                                                                                                                 >  1:858265/1‑239 (MQ=255)
                                                                                                   cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCCAGAAAGGGGATAAATTGCTTTCCAGAGGGTCGGGGGGTTGCCCCAACGCGGCCCGCCGGGGGGCATAGCGGTGGACCATTTTCCCTATCAACCCCCCGAAAAGAACCTTCCCGGGGATTGATTCCCGCCCCCGGCCCCCTAAACGAATATCCaccctggccttt                                                                                                                                 >  2:129666/1‑228 (MQ=255)
                                                                                                   cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGCTAACGCCCAGAAGTCGACTAATTTGCTTTCCAGCGCGTCGGGGGCTTGCGCAAAGTCCGCCTGCCGTGAGCGATAGCGGTTGACCATTTCACTGATAACCGCCCGCAAAGGAACTTCCCCGAGGATGACTTCACACCCCCGGCCCCCTTAACGAATATCCCACCTTCCATTa                                                                                                                                 >  2:176576/1‑239 (MQ=255)
                                                                                                       gataaagagcgcggcgGGTCCACAAATCCATGCTATGGAAGGTCGACATAGCAAGCCCGGTGGTTAACGCAGGAAAGCCGTCTTCTTCTCTTTCCAGGGTGTGGGGGTTGTGGCCACCGTCGGGTCGGGGGGGCCGATAACTATGTCCCATTACACGGATATCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTGGCCATTATCg                                                                                                                              <  2:742607/229‑1 (MQ=255)
                                                                                                                 ggggcGGGTGAAAAAACCCATAATAGGGAAGCGCGACATACAAAGCCTGGCGGTTAGCGCCGAGAGTCCGATTAGTGTGCTTTCCGGCGTGCCGAGGGGTGGCGCAACTCCGGCCGGGGGGGAGCGATGCGAATGTCCCTTTTCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACg                                                                                                                     <  2:779208/234‑1 (MQ=255)
                                                                                                                                            ggAAGCGCGACATATCAAGACCGGCGGATAAAGCCGCGAAATCGAATAATTTTCTTTCCAGCGTGTCGGGGGCCTGCGCAACGGCGGCCTGCGGTGAGCGATAGCGATGTACCATTTTACTTATAACCGGCAGCAAAGGGACTTTCCAGAGGAAGACTTCAAGCCACGGCCCGACTAAACGAATATCCAACTTGCCCATATAGTTGGACCCTGTGGCTTTTTCCGGGTTAAAGCggc                                                                                            >  2:897926/1‑236 (MQ=255)
                                                                                                                                               aGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGGCGACTAATTTGCTTTCCAGCGGGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAAACGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAAAATCCAGCTTGCCCTTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTc                                                                                        >  1:1072813/1‑238 (MQ=255)
                                                                                                                                               aGCGCGACATATCAAGACCGGCGGTTAAAGCCCGGAAGTCGGCTAATTTTCTTTCCCGCGTGTCGAGGGCTTTCGCCACGTCGGCCTGCGGGGTGCGCTAGCGGTGTACCATTTCCCTGATAACCGCCCGCAAAGGAACTTCCCCGAGGAGGACTTTACCCCACCGGCCGCCTTAACGAAAACTCCGCCTGCCCTTATCTTTGTGCACGGGGACTTGTTCCAGGGTAAAGCgggaacc                                                                                        >  2:943205/1‑233 (MQ=255)
                                                                                                                                                     ccATATAAGGAGCGGTGGTAAACGCGGAAAATGCGCTTAATGTGCTTCCCGGGGTTTGGGGGCCTGGCGCAGCGTCGGCCGGCGGTGGGCAATGGCGATGTACCATTTCCCGGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                  <  2:1087425/237‑1 (MQ=255)
                                                                                                                                                     caaagaacaagccGGGGGGTAACGGCCGAGATGCCGAATAATTTGCTTTCGAGGGGTTCGAGGTCTTCGGCAACGTCGGCCGGGGGTGAGCATTAGGTGTGTACCTTTCCCCGGATAGCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                  <  2:576046/228‑1 (MQ=255)
                                                                                                                                                     aCATATAAATCCCGGCGGTTACGGCCGAGAAGTCGACTAATTTGCTTTCCGGGGTGTCGGGGGCTTGCGCAACGGCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                  <  1:34652/238‑1 (MQ=255)
                                                                                                                                                        tatCAAGCCGGGCGGTTACCGCCGAGAAGTCGCTTAATTTGTTTTCCGGCGTGTCGAGGGTTTGCAAAACGTCGGCCTGCGGTGAGCAATAGCAAGGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAg                                                                              <  2:1072813/239‑1 (MQ=255)
                                                                                                                                                                   gcggTAAACGCCGAGAAGTCGACTAATTTGCGTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGTTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                    <  1:347951/238‑1 (MQ=255)
                                                                                                                                                                                    agtgtgatagtttGATGTCCAGGGTTTCGAGGGTTTGCGCAACGCCGGCCGGCGTGGAGCAAAAACAATGTACCATTTCCCTAAAAGCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggc                                                   <  2:1060354/229‑1 (MQ=255)
                                                                                                                                                                                                tGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAg                                        >  1:490804/1‑238 (MQ=255)
                                                                                                                                                                                                ggCTTTCCCGCTTGTCGAGGGCTTGCCGACCAGCGGCGTGGGCAGCGAGATAGAGATGTCCCATTTCTATGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc                                       <  1:329344/238‑1 (MQ=255)
                                                                                                                                                                                                   ttcccggCGTGTGGTGGGCTTGCCCAACGTCGCCCGCGGGTGAGCGATAGCAAGTTACATTTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCa                                     <  1:689224/232‑1 (MQ=255)
                                                                                                                                                                                                   ttccaagggtgcCGAGGTTTTGCACAACGTCGGCCTGCGGTGGGCGATACCATTTCACCATTGCACTTATAGCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCAc                                    <  1:1031769/231‑1 (MQ=255)
                                                                                                                                                                                                    ttccacgtggtCGAGGGTTTGCACAACGTCGGCCTGGGGTGAGCGATAGCAATTATCCATTACACTTAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACt                                   <  1:553835/230‑1 (MQ=255)
                                                                                                                                                                                                    ttCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGcc                                      >  1:898767/1‑236 (MQ=255)
                                                                                                                                                                                                       cggcggttCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGat                                <  1:934470/232‑1 (MQ=255)
                                                                                                                                                                                                             gtCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAACGATGTACCATTTCACTGCTAACCGCCAGCAAAGGCACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGCATATCCAGCATGCCATTATCGTTGGACACGGTGACTTGTTACGGGTTAAAGCGGAACTCGAGTAACCAGTTAAGATAGTCGGCCTTACAGACAGGCAGGGCGGAAAGGCACTGAACTtcat                          >  2:666170/1‑239 (MQ=255)
                                                                                                                                                                                                                     cTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGc                  >  1:1067453/1‑239 (MQ=255)
                                                                                                                                                                                                                      ttGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCa                 >  2:353335/1‑239 (MQ=255)
                                                                                                                                                                                                                        gcgcAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                >  1:904207/1‑238 (MQ=255)
                                                                                                                                                                                                                               cGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGCAAGCCACTGATATTCATCATCCTGTAGGCGCAgg         >  1:1061032/1‑239 (MQ=255)
                                                                                                                                                                                                                                      cggCGGTGAGCGATTGCTATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  1:1079509/237‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGGGGTGGGCGCTAGCGATGTACCTTTTCACTGCTACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  2:257671/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGGGGTGAGGGATAGCGATGTCCAATTTCGCGAATAACCGCTAGCAAAGGAACTTCCAAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  1:323069/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGCGGTGAGCATTAGCAATATCCCATTTCACTTATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  1:174618/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGCGGTGAACGATAGCGATGTACCATTTCACTGATAGCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  2:858265/239‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCGGTGAGAGAGAGCGAGGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  1:589659/237‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
CCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCA  >  NC_000913/989711‑990175

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: