Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,788,475 C→T G125G (GGC→GGT yibQ → divergent polysaccharide deacetylase domain‑containing protein YibQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,788,4750CT98.2% 184.4 / NA 57G125G (GGC→GGTyibQdivergent polysaccharide deacetylase domain‑containing protein YibQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/36);  minor base A (1/0);  total (21/36)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCGATGCCTTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCGTTGAGTAAACAGCCGCTGGAGAAAAATACGCTACGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTATTCCTCGACGATTCGCAAAATGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGACCTGGCGCGTCGCAACGGTTCGACC  >  NC_000913/3788230‑3788712
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
cTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGgtagca                                                                                                                                                                                                                                           >  2:274655/1‑250 (MQ=255)
    cGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATg                                                                                                                                                                                                                                      <  2:329874/251‑1 (MQ=255)
     gATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGTATCAACAAGCACATGGGTAGCAAGATGa                                                                                                                                                                                                                                     <  2:269375/251‑1 (MQ=255)
      aTGCCCTCCGCTATCTCCGGCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGAc                                                                                                                                                                                                                                    <  2:153591/251‑1 (MQ=255)
         ccctCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGAAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCATATGCCGTGGGGATCAACAACCACATGGGAAGCAAGATGACctc                                                                                                                                                                                                                                 >  1:96977/4‑251 (MQ=255)
          cctCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACctct                                                                                                                                                                                                                                <  1:178606/249‑1 (MQ=255)
             ccGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTaa                                                                                                                                                                                                                              <  2:171976/250‑1 (MQ=255)
                    ctcCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGCCCTCTAACCTGTTTg                                                                                                                                                                                                                      >  1:25908/1‑251 (MQ=255)
                      ccGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGttt                                                                                                                                                                                                                       <  1:39872/248‑1 (MQ=255)
                                     ccGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATg                                                                                                                                                                                                     <  1:275195/251‑1 (MQ=255)
                                      cGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGc                                                                                                                                                                                                    >  1:284931/1‑251 (MQ=255)
                                           cACCGCACGCCAGAGGGAGGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCATCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCTTCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGGAGTGCGGTCAAGAACGTGCCTTATGCCGTGGGAATCACCAGCCACATGGGTAGCAAGATGGCCTCTAGCCTGTTTGGTATTCAGAAAGTGATACAGgc                                                                                                                                                                                                <  2:413293/250‑1 (MQ=255)
                                                     cAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGtt                                                                                                                                                                                     <  2:240877/251‑1 (MQ=255)
                                                     cAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGAGTCATCTCCCGAGGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGtt                                                                                                                                                                                     <  2:375220/251‑1 (MQ=255)
                                                     cAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGGTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGAGCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGtt                                                                                                                                                                                     <  2:58687/251‑1 (MQ=255)
                                                      agagagATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTa                                                                                                                                                                                    <  2:400136/251‑1 (MQ=255)
                                                      agagagATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTa                                                                                                                                                                                    <  1:32234/251‑1 (MQ=255)
                                                       gagagaTGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGTGTTAc                                                                                                                                                                                   <  2:284931/251‑1 (MQ=255)
                                                              ggCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCttt                                                                                                                                                                            <  1:14222/251‑1 (MQ=255)
                                                                cGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCAGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTAc                                                                                                                                                                          >  1:101628/1‑251 (MQ=255)
                                                                   ccAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTc                                                                                                                                                                       >  2:16452/1‑251 (MQ=255)
                                                                     aaaGCTCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCt                                                                                                                                                                     <  2:25908/251‑1 (MQ=255)
                                                                     aaaGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATGCGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTGTACTTCCt                                                                                                                                                                     <  2:144163/251‑1 (MQ=255)
                                                                      aaGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTc                                                                                                                                                                    <  1:409862/251‑1 (MQ=255)
                                                                          gcATAACCGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCACGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAACGTGATGCCGGCACTGGAGCGTTACAATCTTTACTTCCTCGAc                                                                                                                                                                 >  2:266222/1‑250 (MQ=255)
                                                                                cAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTaa                                                                                                                                                          <  1:289611/251‑1 (MQ=255)
                                                                                           aaGTGTTGCTTCATCTCCAGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATACGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGCCCTCTACCCTGTTTTGTCTGCAGAAAGTGATGCAGGCACTGTAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAAt                                                                                                                                               >  1:271234/1‑251 (MQ=255)
                                                                                               gtTGATTCCTCTCCCGCCGGCAACACTGAGTAAAAAGCCGCTGGAGAAAAATACGCTGCGCCCGGCGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGCAAGTGATGCAGGCACTGGAGCGTTACAATCTTTCCTTCCTCGACAGCGTAACCATCGGTAATAccc                                                                                                                                           >  1:299311/1‑251 (MQ=255)
                                                                                                              gATGGCACCACTGAGTAAACAGCCGCTGGAGAAAACTACGCTGCGCCCGGAGATGAGCAGCGACGCAATTGAGCGCATTATTCGTAGCGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGcgcgccg                                                                                                                            >  1:20887/1‑251 (MQ=255)
                                                                                                                    accactGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTc                                                                                                                                                                       <  1:152494/197‑1 (MQ=255)
                                                                                                                    accactGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTc                                                                                                                                                                       >  2:152494/6‑202 (MQ=255)
                                                                                                                                 cAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGa                                                                                                          >  1:397569/1‑250 (MQ=255)
                                                                                                                                  aGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGa                                                                                                          >  2:171257/1‑249 (MQ=255)
                                                                                                                                         ggAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGAGCACCAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCGTTACTTCCTCGACAGCGGAACCATCGGTAATACCCAGGCGATGAGCGCCGCGCTAGGTACTGGCGTGaattggatca                                                                                                 <  2:385750/251‑9 (MQ=255)
                                                                                                                                           agaAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCaa                                                                                                >  1:55652/1‑250 (MQ=255)
                                                                                                                                            gaAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAAc                                                                                              <  2:335228/251‑1 (MQ=255)
                                                                                                                                              aaaaTACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACgg                                                                                            <  2:63326/251‑1 (MQ=255)
                                                                                                                                                 aTACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAg                                                                                         <  1:274655/251‑1 (MQ=255)
                                                                                                                                                 aTACGCTGCGCACGGCGATGAGCAGCGACGAAATTGCGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGACGGTGATCAAACGGa                                                                                           >  2:275531/1‑249 (MQ=255)
                                                                                                                                                    cgctgcgcCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGtg                                                                                      >  2:82789/1‑250 (MQ=255)
                                                                                                                                                      ctgcgcCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGGGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTACCCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTAATGGCGTGAAGGTGATCACACGGCAggggtt                                                                                    >  2:11314/4‑247 (MQ=255)
                                                                                                                                                         cgcCCGGCGATGAGAAGCGACGAAATTGAGCTCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGGGGGGATCAACAACCACAGTGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATGCCCAGGCGATGCGCGCCGCGCAAGGTACGGGCGTGAAGGTTATCAACCGGAAGTTGTTCCt                                                                                 <  2:210554/251‑1 (MQ=255)
                                                                                                                                                         cgcCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCt                                                                                 >  2:312149/1‑251 (MQ=255)
                                                                                                                                                                  attaggagcGGCGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGGGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTc                                                                        <  2:384996/245‑1 (MQ=255)
                                                                                                                                                                    gAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGc                                                                      <  1:312149/251‑1 (MQ=255)
                                                                                                                                                                            acgaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCagaacgaa                                                              <  1:171257/251‑8 (MQ=255)
                                                                                                                                                                                    gAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACAt                                                      <  1:357701/251‑1 (MQ=255)
                                                                                                                                                                                       cgcATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCg                                                   <  1:16452/251‑1 (MQ=255)
                                                                                                                                                                                          attattCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCgtgt                                                <  1:152713/251‑1 (MQ=255)
                                                                                                                                                                                                cGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAGCCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAAt                                           <  1:52138/250‑1 (MQ=255)
                                                                                                                                                                                                cGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTGACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAAtt                                          <  1:275531/251‑1 (MQ=255)
                                                                                                                                                                                                 gTAGTGCGGTCAATAACGTGCACTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGGGGGCAAttt                                         >  2:258657/1‑251 (MQ=255)
                                                                                                                                                                                                   agtgcgagcAATAACGTGCCCTATGCCGTGGGGATCAACAGCCACAGGTGTAGCAAGAGGACCTCTAACCTGTGTGGTATGCATAAAGTGATACAGGCACTAGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTGCTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGATCGAAGCGGACCTCCGTGTTCAATTTaa                                       <  2:96977/243‑1 (MQ=255)
                                                                                                                                                                                                    gTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAAAGAAGAGGACATCAGTGTGCAATTTAAt                                      >  1:364967/1‑251 (MQ=255)
                                                                                                                                                                                                           cAATAATGTGCCCTATGCCGTGGGGATAAACAAACACATGGGTCGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAAGGCCCAGTCGATGCGCGCCGCGCAAGGTACTGGCGTTAAGGTGATTAATCGGAAGTTGTTCGTCGACTATTCGCAGAACGAAGCTGACAGCCGTGTGCAATTTAATCTCGCaa                               <  2:437333/251‑1 (MQ=255)
                                                                                                                                                                                                                  gTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGatct                        <  1:258657/251‑4 (MQ=255)
                                                                                                                                                                                                                     cccTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGgc                     >  2:216387/1‑251 (MQ=255)
                                                                                                                                                                                                                                        acaacCACAAGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGCTCTGGCGCGGCGCAACGGTTCAAcc  >  2:132737/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
CTGGCGATGCCTTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCGTTGAGTAAACAGCCGCTGGAGAAAAATACGCTACGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTATTCCTCGACGATTCGCAAAATGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGACCTGGCGCGTCGCAACGGTTCGACC  >  NC_000913/3788230‑3788712

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