Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 989,890 T→G 22.2% D156A (GAC→GCC)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,8900TG22.2% 84.7 / 3.0 56D156A (GAC→GCC) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (11/25);  new base G (9/4);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 2.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGGCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGT  >  NC_000913/989662‑990109
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                
tCCGGGCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAt                                                                                                                                                                                                                          <  2:1172103/237‑1 (MQ=255)
   gggCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCt                                                                                                                                                                                                                     <  1:234480/239‑1 (MQ=255)
         aTTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTAAGCCGGGGAATCGGCTAATTTTCttttcc                                                                                                                                                                                                                >  1:1030287/1‑235 (MQ=255)
         aTTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGCTAACGCCACGAAGTCGGCTAATTTTCTTTccc                                                                                                                                                                                                                >  1:59029/1‑237 (MQ=255)
           ttGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCg                                                                                                                                                                                                            <  2:158107/239‑1 (MQ=255)
                     gCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGc                                                                                                                                                                                                  <  1:91808/239‑1 (MQ=255)
                     gCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGc                                                                                                                                                                                                  <  1:813202/239‑1 (MQ=255)
                      ccTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGACAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGc                                                                                                                                                                                                  <  1:92137/238‑1 (MQ=255)
                                ttCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGAAATATCAAGACCGGCGGTTAAAGCCGAGAAAGCGGCTAATTTGGTTTTCAGCGGTGGCGAGGGTTACGCAAccg                                                                                                                                                                                      >  1:496772/1‑237 (MQ=255)
                                ttCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGGTTCTGGTGCCCACAAACCAGGATTCCTGTTGCACACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCACGACCGGCGGTTCAAACCCAGACATACACTACTTTTCTTTCCAGGGTGTCGCGGCCTTtcgccccgt                                                                                                                                                                                      >  2:166748/1‑234 (MQ=255)
                                           gtgtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTAAACCGCCAAGTGGGAAAAATTGTTTTCCCACGCGTCGGGGGCATgcggcaccgccgcccgtggg                                                                                                                                                                           >  2:838957/1‑222 (MQ=255)
                                            tgtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAATCGACTAATTTGCTTTCCCGCGTGTCGGGGGGTTGGGGAACGCCGGccccggggg                                                                                                                                                                          >  1:257671/1‑232 (MQ=255)
                                             gtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTg                                                                                                                                                                          <  2:774788/238‑1 (MQ=255)
                                              ttCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAt                                                                                                                                                                                                                          <  1:637580/191‑1 (MQ=255)
                                                 ccgtCGGCGCGAGGGAAAGCCGACGCGCCAGTTCGTAGTTGCTTGTGCCCCGAAACCAGGATTCCTGTTGCGGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAACGCGGCGACGGGTGCCAATATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                     <  2:1051987/236‑1 (MQ=255)
                                                 ccacgggCGGGAGGGAAAGCCGACGCGCCAGAACGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTGGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGGCGACTAATTTGCGGTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                     <  2:336952/234‑1 (MQ=255)
                                                 ccATCGGCGTGACGGGAAGCCGGCGCGCCAGAGCGCAGTTGCTGGTGCTCACAAACCAGGATTCCTGTCGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                     <  2:168200/239‑1 (MQ=255)
                                                   aTCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAAAGCCCAGAAAGCGACAAATTTTCTTTTCAGCGTGGCGGGGGGTTGGGCCAAGTCgggcctcggggaggggta                                                                                                                                                                   >  1:383272/1‑223 (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)
                                                                                                                                                    cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGATAACGCCGGCAAGGCGACGAATTTGATTTCCACCGCGGCGGGGGGTCGCGGAAAGCCGGCCTGCGGTGTGGGGTAGAGGGGTAGCCTTTCACCGATAACCCCCCCGCAGGGAAAGTACCAAAGGGTGACTTTACGCCACCGACCGGCTAAACGAaaataactacaggcaatg                                                               >  2:589659/1‑225 (MQ=255)
                                                                                                                                                       cgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAACCCGCGAAGACGACTACTTTGGTTTTCAGCGGGTCGGGGGCTTGCGCAACATCGGCCCTCGGGGAGCGATAGCGAAGGACCATTTTACTGCTAACCACCAGCCAAAGAACATCCCCCAGGAGGAATTCACGCCCACGGCAGCATAAAcaacgaatcatgtctcccacacg                                                             >  2:1191393/1‑217 (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)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                
TCCGGGCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGT  >  NC_000913/989662‑990109

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