Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 708,774 G→A A147A (GCG→GCA chiP → chitooligosaccharide outer membrane channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913708,7740GA98.0% 157.9 / NA 51A147A (GCG→GCAchiPchitooligosaccharide outer membrane channel
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/27);  minor base C (1/0);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCTCACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACG  >  NC_000913/708532‑709018
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
cTTTCTCACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAgg                                                                                                                                                                                                                                              >  1:518015/1‑251 (MQ=255)
          cccccTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTc                                                                                                                                                                                                                                      <  2:13941/246‑1 (MQ=255)
            cccTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGcc                                                                                                                                                                                                                                  <  1:637953/250‑1 (MQ=255)
             ccTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTAAGGCGATTGAAATGGCGGAAAACGGCGACATCTCCCACCCGAACGAAATCGCGTTTTAAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATCAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGGTTGGGCAAGGGCAGGTTACATTCAGcc                                                                                                                                                                                                                                  >  1:646786/1‑250 (MQ=255)
                   aTGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGt                                                                                                                                                                                                                           >  1:542433/1‑251 (MQ=255)
                             cGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCt                                                                                                                                                                                                                   <  1:623330/249‑1 (MQ=255)
                                ttttCAGTCCGGTTAGGATGCTGATGTGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAACTGGCGGTAAAAGCCGCCAGCTCCCACCCGTACGAACTCGCGTTTTAACAAAGTAATAAAGCCTATGACAAAGACTGGTCCGGCAGCAACAGCAGTCTAAGACTGTATAAAGCTGCGCACACATTTAAATAAGTGACGGTTTGGACCAGGGCAGGTTACAtgacgccaaacagtaaaacactatt                                                                                                                                                                                                                >  2:39479/1‑225 (MQ=255)
                                              aTGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGc                                                                                                                                                                                                <  1:601313/251‑1 (MQ=255)
                                                 ctgctgATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTGATAAAGCCTAAGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCttt                                                                                                                                                                                             >  1:43691/1‑251 (MQ=255)
                                                  tgctgATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCtttt                                                                                                                                                                                            >  1:183778/1‑250 (MQ=255)
                                                   gttgATATGTGCGGCCTTGATATGGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTAt                                                                                                                                                                                           <  2:542433/249‑1 (MQ=255)
                                                              cGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTAAGCCAACTGGTCAAACGCTGTTAGCGCCGCCCTGGAGCTTTTTGCCAGGTACtt                                                                                                                                                                                >  2:11318/1‑251 (MQ=255)
                                                               ggCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAACACGGCGACATATCCCACCCGAACGAAATCGCGTTTTCAAACAGTAATAAAGCCTATGACTAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGTTGCGGCCAAATTTAAAAACGGTCCGGTTTGGGCAAGGGCAGGGTGCATTCAACCACCTGGTCAAACGCTGTTAGCGCCGCACTGGTGCTTTATGACAGGtccctt                                                                                                                                                                               >  2:632877/1‑246 (MQ=255)
                                                                    tGATATTGCTGCTTGTACGGCGAATGAAATGGCAGAAAACGGCGACAGCTCCCTCCCTAACGAAATCTCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTTGTCCGGCGAAAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAGTTTAAATACGTTCCGGTTTGGGCAAGGGCAGGTTACATCCAGCCAACTGGTCAAACTGAGTTAGCGCCGCACTGGTGCTTTATGCCAGATAGTTATCAg                                                                                                                                                                           <  2:254276/250‑1 (MQ=255)
                                                                                gTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCg                                                                                                                                                              >  1:615510/1‑251 (MQ=255)
                                                                                 tttACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAACGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCgg                                                                                                                                                             >  1:649857/1‑251 (MQ=255)
                                                                                 gttACGGCGATTGATATCGAGGAAAACGGCGACAGCTACAACCCGAACGAAATCACGTTTTCAAAAAGTAAAAAAACCTAGGACGAGGACTGGTCAGGCGACAAAAGCGGTATAAGACTGTATAAAGCTGCGGCCAAATCTAAATACTGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATTCCAGGTACTTATCAGGGGGCGGAAGcc                                                                                                                                                               <  2:211232/248‑1 (MQ=255)
                                                                                  ttACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCggg                                                                                                                                                            >  1:94247/1‑251 (MQ=255)
                                                                                         gATTGAAATGGCGGAAAGCGGCGACAGCTCCCACCCGAACGAAATCGCGGTTTAAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGGATAAGCCTGTATAAAGCTGCGGCCAAATTGAAATACGGTCAGGTTTGGGCAAGGGCAGGTTACCTTCAGCAAACTGGTCAAAAGCTGTTAGCGCCGAACTAGGGCTTTATGCCAGGGACTTATCAGGGGGCGGAAGCCGGGGCGAAtt                                                                                                                                                     >  2:624193/1‑251 (MQ=255)
                                                                                         gATTGAAATGGCAGAAAACGGCAACAACTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCTGTATAAGCCTGTATAAAGCTGCGGCCCAATTTAAATAAGTTCCGGTTGGGGCAAGGGCAGGTTACATTAAGCCAACTGGTCAAACGCTTTTATCGACGCACTGGAGCTTTTTGCCAGGGACTGCTCCGTGGGCGGACGCCGGGGGGAAtt                                                                                                                                                     >  2:362938/1‑251 (MQ=255)
                                                                                              aaaTGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGAt                                                                                                                                                <  1:362938/251‑1 (MQ=255)
                                                                                                aTGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTa                                                                                                                                              <  1:120984/251‑1 (MQ=255)
                                                                                                  ggCGGAAAACGGCGACAGCACCCACCCGAACGAAATCGCGTTTTCAAAAATTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTAc                                                                                                                                             <  2:352364/250‑1 (MQ=255)
                                                                                                                    cTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTg                                                                                                                           >  1:131803/1‑250 (MQ=255)
                                                                                                                     ccccACCCGAACGAAATCTCGTTTTCAAAAAGTAATAAAGCCTATGATGAAGTCTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGTTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTGAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAg                                                                                                                         <  2:518015/250‑1 (MQ=255)
                                                                                                                         aCCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGa                                                                                                                                                                   <  1:305199/205‑1 (MQ=255)
                                                                                                                         aCCCGAACGAAATAGCGTTTTAAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGACAAATTTAAATAAGGTACGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGAGACGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGa                                                                                                                                                                   >  2:305199/1‑205 (MQ=255)
                                                                                                                           ccaaactaaagcgcgTTTTCAAAAAGTAATAAAGCCTATGACGAAGCCTGGGCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATGCGGTCCGGTTGGGGCAAGGGTATGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCTCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGTTGCCGGGGCGAATTTTGATTATGGCGATGCTGGTGCGTTGAGTTtctc                                                                                                                   <  2:377133/240‑1 (MQ=255)
                                                                                                                                 gaaATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACAt                                                                                                             >  1:65045/1‑251 (MQ=255)
                                                                                                                                    aTCGCGTTTTCAAAAAGTAGTAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGATGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACAtttg                                                                                                          >  2:552933/1‑248 (MQ=255)
                                                                                                                                                          aaGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGtt                                                                                                                       <  2:93749/216‑1 (MQ=255)
                                                                                                                                                          aaGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGtt                                                                                                                       >  1:93749/1‑216 (MQ=255)
                                                                                                                                                                  gACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATc                                                                             <  1:11318/250‑1 (MQ=255)
                                                                                                                                                                    cGAAGACTGGTCAGGCGCCAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTTCTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACAAACGAATACAAAGCGCCGTGGcttcggt                                                                          <  2:679276/251‑7 (MQ=255)
                                                                                                                                                                      aaGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGaa                                                                        <  2:615510/251‑1 (MQ=255)
                                                                                                                                                                         aCTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATg                                                                     <  1:632877/251‑1 (MQ=255)
                                                                                                                                                                            ggTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGAt                                                                  <  1:426087/251‑1 (MQ=255)
                                                                                                                                                                                       aaaaGCGGTATAAGCATGTATAAAGCGGCGGCCAAATTTAAATACGGACCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTATGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTAtaa                                                       <  2:183778/251‑3 (MQ=255)
                                                                                                                                                                                        aaaGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAACTAAGGTCCGGTTTGGGCAAGGGCAGGTTACCTTCAGCCAACTGGTCAAACGCTGTTAga                                                                                                                                                                                                             >  2:331196/1‑99 (MQ=255)
                                                                                                                                                                                        aaaGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAgc                                                                                                                                                                                                             <  1:331196/100‑1 (MQ=255)
                                                                                                                                                                                             ggTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACg                                                  >  1:555617/1‑250 (MQ=255)
                                                                                                                                                                                                     ccTGTATAAAGCTGCGGCCAAACTGAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAacca                                         <  1:196265/251‑1 (MQ=255)
                                                                                                                                                                                                                   cGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTAt                           >  1:675179/1‑251 (MQ=255)
                                                                                                                                                                                                                       cAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATcagc                       >  2:570863/1‑248 (MQ=255)
                                                                                                                                                                                                                                     gtacGGTTTGGGAAAGTGCTGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCTCGGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGTCCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCg         <  2:498055/248‑1 (MQ=255)
                                                                                                                                                                                                                                       ccGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTgg                                                                                                         <  2:495761/153‑1 (MQ=255)
                                                                                                                                                                                                                                       ccGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTgg                                                                                                         >  1:495761/1‑153 (MQ=255)
                                                                                                                                                                                                                                        cGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATc                                                                             <  1:480130/180‑1 (MQ=255)
                                                                                                                                                                                                                                        cGGTTTGGGCAAGGGCAGGGTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGCAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATTTGGACCAACGAATACAAAGCGCCGTGGCATc                                                                             >  2:480130/1‑180 (MQ=255)
                                                                                                                                                                                                                                          gTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATa    <  2:555617/251‑1 (MQ=255)
                                                                                                                                                                                                                                             tGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACg  <  2:65045/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
CTTTCTCACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACG  >  NC_000913/708532‑709018

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