Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 71,551 2 bp→TC coding (201‑202/765 nt) yabI → DedA family protein YabI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,5510CT94.4% 201.4 / ‑0.2 73C67C (TGC→TGTyabIDedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (16/51);  total (18/55)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
*NC_00091371,5520TC91.7% 218.2 / 1.3 73L68L (TTG→CTG) yabIDedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base C (15/51);  total (18/55)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCG  >  NC_000913/71333‑71779
                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                       
ccATTTGATAGGCACATTATGCATGCATTGCAGTAACTATTTATTACCCAATCTACCGTGTATTCAGTGATGGCGGTGGTGTTGGTGGCCGTTCTGGAGGCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGATGGTAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGTCTGGATCTCCTTctggatgg                                                                                                                                                                                                         <  2:517905/250‑1 (MQ=255)
ccATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTctggctgg                                                                                                                                                                                                         <  1:507518/250‑1 (MQ=255)
     tgtcAGGCACATTATGCAAGCATGTCTGGAACACTTTATTACCCAACCTACCGTGTATTCATGGATGGCGGTGGTGTTGGTGGACTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATGGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGgctgggctgg                                                                                                                                                                                                    <  2:266521/247‑1 (MQ=255)
        cAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCGCCTTCTGGCTGGGCTGGCGtt                                                                                                                                                                                               <  1:279793/251‑1 (MQ=255)
              caTTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTTCCGTGTATTCATTGATGGCGGTGGTGTTTGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTaaaa                                                                                                                                                                                         <  1:649164/251‑1 (MQ=255)
                     cAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGtt                                                                                                                                                                                  <  1:545384/251‑1 (MQ=255)
                      aaGCATTGCTGGAACACTGTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGGGCGCTGATTGGCAGCGGCGGGGTTAAGTTTTGGGCCGCCTGGCTGGGAGGGGTTGGCGGGTTTCTTCTGGGGGACTGGATTTCCTTCTGGCTGGGGCGGCGGTTTTAAAAACCGTTg                                                                                                                                                                                 >  1:542225/1‑251 (MQ=255)
                           ttgcgtgAACACTTTTTTACCCAAGCTAGCTTTTATTCATTGATGTCTGTGGTTGTGGTGGCCTTTCTGGTGTCGCTGGCGCTGGTCGGTTTTTGTCTTCCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGTCAGCGGCGAGTTAAGTTTCTGGCACACCTGGCTGTCAGGGATTGTCGGTTGTCTGCTGGTCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCAt                                                                                                                                                                              <  2:787465/243‑1 (MQ=255)
                             gCTGGAACATTTTATTACCCAATCTACCGTGTATGCATTGATGGCGTGGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGTTTTGATTCTACCCGGTGCGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTGGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCTTTTTAAAAAGCCGTTGCATCGCt                                                                                                                                                                          <  2:24023/251‑1 (MQ=255)
                                ggataacgttatttgcagatctcccGTGTAGTCATTGGTGGCGGTGGTGTTGGTGGCCTTTCTGGAGGCTCTGGCGCTGGTCGGTTTGATTCTGCCCGGGACGGGGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACTCATGGCTGGCAGGGATTGTCGGTTGTCTTCTGGGCGACTGGATCGCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGt                                                                                                                                                                       <  2:422161/228‑1 (MQ=255)
                                   atacttgttaaCCCAATCGTCCGTGTATTCATTGATGGCGGTGGGGTTGGGGGCCTGTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAATTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCTCTGGTCAt                                                                                                                                                                    <  2:383256/241‑1 (MQ=255)
                                    cacTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCAtt                                                                                                                                                                   <  1:665083/251‑1 (MQ=255)
                                        ttattaCCCAATCTACGGTGTAGTCGTTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCt                                                                                                                                                                <  2:361317/250‑1 (MQ=255)
                                          attaCCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTgaa                                                                                                                                                             <  1:559130/251‑1 (MQ=255)
                                           ttacccatactaCCGTGTATTCATTGATGTTGGTGGTGTTGGTGGCCTGTCTGGAGTCGCGGGCGCTGGGCGGTTTGATTCTACCCGGTACGGTGCTGAGGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTgaa                                                                                                                                                             <  2:914917/239‑1 (MQ=255)
                                                      aCCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGc                                                                                                                                                 <  1:390593/251‑1 (MQ=255)
                                                        cGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCat                                                                                                                                               <  1:292347/251‑2 (MQ=255)
                                                         gtgtATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTGCGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCatt                                                                                                                                              <  1:390956/251‑3 (MQ=255)
                                                          tgtATTCGTTGATGGCGGTGGTGTTGGTGGCCTTTCTGTGGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGAGTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCatt                                                                                                                                              <  2:838870/250‑3 (MQ=255)
                                                           gtATTCATTGATGGTGGTGGTGTTGGTGTCCTTTCTGGAGTCGCGGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAGCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCatta                                                                                                                                             <  1:978895/250‑4 (MQ=255)
                                                             agtCATTGATGGCTGTGGGGTTGGGGGCCTTGCGTGAGTCGCGGGCGCTGGTCGGTTTGATTCCACCCGGTACGGTTCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGAGGCGAGTTAAGTTTCTGGGACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCTTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCattactc                                                                                                                                          <  2:312744/249‑2 (MQ=255)
                                                               tCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGGTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCattactcg                                                                                                                                         <  1:774660/250‑3 (MQ=255)
                                                                cATTGATGGCGGTGGTGTTGGTGGTCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGAGGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATATCTGAAGAAAAACAAAGCATTACTCGAt                                                                                                                                       <  2:693328/251‑1 (MQ=255)
                                                                cATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGAt                                                                                                                                       <  1:299936/251‑1 (MQ=255)
                                                                 aTTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATa                                                                                                                                      <  1:370362/251‑1 (MQ=255)
                                                                   ggagggggggtgtgtTGGTGGCCTTGCTGGAGTCGCTGGCTCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCTACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATaa                                                                                                                                     <  2:935048/239‑1 (MQ=255)
                                                                        gCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCGGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGa                                                                                                                               <  1:967714/251‑1 (MQ=255)
                                                                          ggtggtGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTTGGCGACTGGATCTCCTTCTTGCTGGGCTTGCGTTTTTAAAAACCGTTGCATCGCTTGTCCTTTTTTAAGAAAAACAAAACATTACTCGATAAAACTgcac                                                                                                                             >  1:431060/1‑248 (MQ=255)
                                                                             ggtgTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGGTTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACAt                                                                                                                           <  1:141123/250‑1 (MQ=255)
                                                                             ggtgTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATg                                                                                                                          <  1:106690/251‑1 (MQ=255)
                                                                                   ggtggCCTTTCTGGGGTCGCGGGCGCTTGTCGGTGTGATTCTACCCGGGACGGGGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGGTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGttt                                                                                                                     <  2:410135/250‑2 (MQ=255)
                                                                                       gCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCa                                                                                                                <  1:698445/251‑1 (MQ=255)
                                                                                        ccTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGGTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGGGTTTTAAAAAGCCGTTTCATCGCTGGTCATTTTTGAAGAAAAAAAAAGCATTACTCGATAAAACTGAACAAGCGTTGCATCa                                                                                                                >  1:876822/1‑250 (MQ=255)
                                                                                              gggAGTCGCGGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGc                                                                                                         <  1:50472/250‑1 (MQ=255)
                                                                                                      ctggcgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGGTTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGGTGGCGTTTTTAAAAGCCGTTGCATCGCTGGGCATTTCTGTAGAAAAAAAAAGCATTACTCGATAAAAATGAACATGCGTTTCATCAACAAAGCATGTTCAc                                                                                                >  1:203310/1‑251 (MQ=255)
                                                                                                             tggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGGGGTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGGTCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGACGTTGCATCGCTGGTCATTTCTTAAGAAAAACAAAGCATTACTCGGTAAAACTGAACATGCGTTGCATCCAAACAGCATGTTTACCATTCTg                                                                                         >  1:394733/1‑251 (MQ=255)
                                                                                                              ggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCt                                                                                          <  1:766355/249‑1 (MQ=255)
                                                                                                               gTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTggt                                                                                       <  1:129597/251‑1 (MQ=255)
                                                                                                                 cGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTggtc                                                                                      <  1:410903/250‑1 (MQ=255)
                                                                                                                  ggTTTGATTGTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGTCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCGTCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTTATCATGCGTTGCATCATCACAGCATTTTCACCATTCTggtcgg                                                                                    <  1:778830/251‑1 (MQ=255)
                                                                                                                       gATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGgtcggtcgc                                                                                <  1:917478/250‑2 (MQ=255)
                                                                                                                          tCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGcttt                                                                             <  1:903995/250‑1 (MQ=255)
                                                                                                                             aCCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTttgtt                                                                          <  1:843308/250‑1 (MQ=255)
                                                                                                                                         gCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGt                                                              <  1:704680/250‑1 (MQ=255)
                                                                                                                                          gtGATGGCGGGGCTGGGGGCGCTGTGTGGCAGCGGCGAGTGAAGTTTCGGGCACTCCTGGCTGGCAGGGATTGTCGGTTGGCTGCTGGGTGACTGGACCTCCGTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACAGGCGTTGCATCAACCCAGCATGTTCACCATTCTGGTCGGGGGCTTTGTTGGCCCGTCGCGTcc                                                            <  2:926307/250‑1 (MQ=255)
                                                                                                                                          cgtagggCGGGACTTGGAGCCCTGATTGGCAGCGTCGAGTTAAGTTTCTGGCACGCCTGGCTGGGAGGGAGTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTCCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATAACTCTATAAAACTGATCATGCGTTGCATCAACACAGCATGTTCCCCATTCTGGTCGGTCCCTTTGTTGGCCCGACGTGTcc                                                            <  2:475029/246‑1 (MQ=255)
                                                                                                                                                  ggggCTGGGAGCGCTGATTGGCAGCGGCGCGTTAAGTTTCTCGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGc                                                    >  1:756836/1‑251 (MQ=255)
                                                                                                                                                  ggggCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGc                                                    >  1:882472/1‑251 (MQ=255)
                                                                                                                                                    ggCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGcc                                                   >  1:344918/1‑250 (MQ=255)
                                                                                                                                                    ggCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCa                                                  <  1:839365/251‑1 (MQ=255)
                                                                                                                                                         ggAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACACCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCCTTTCTGAAGAAAAACCACGCATTAATCGATAAAAATGAACATGCGTTGCTTCAACCCAGCATGTTCCACCTTCTGGTCGGTCGCTTTTTTGTACCGACGCGACCCCTGGAGACAAtggt                                             >  2:865995/1‑251 (MQ=255)
                                                                                                                                                           aGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAAtggtgg                                           <  1:358915/251‑1 (MQ=255)
                                                                                                                                                                tGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGGTTGTCGGTTGTCTTCTGGGCGACTGGATATCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCCGGTCATTTCTGAAGAAAAACAAAGAATTACTAGATAAAACTGAACATGCGTTGCATAAACACACCCTTTTCACCATTCTGGTCGGTCGCTTTTTTGGCCCGCCGCGGCCGCTGGTGCCAATGGTGGCgggg                                      >  2:572563/1‑250 (MQ=255)
                                                                                                                                                                 gATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCggg                                       <  1:84598/249‑1 (MQ=255)
                                                                                                                                                                   ttggtagtggtgAGTTAAGTTTTGGGCACGACTGGCTGGCAGGGAGTGTCGGTTTTTTGCTGGGCGACTGGATCGCCTTCTGGCTGGGCTGGAGTTTTAAAAAGCCGTTGCATAGCTGGTCATTTCTGAAGAAAAACAAAGTATTTCTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAAt                                    <  2:314580/239‑1 (MQ=255)
                                                                                                                                                                    tGGCAGCGGCGAGTTAAGTTTCTTGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACAGGCGTTGCATCAACACAGCATGTTCACCATGCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATg                                   <  1:909459/250‑1 (MQ=255)
                                                                                                                                                                     ggCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTGCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATGACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCt                                 <  2:935128/251‑1 (MQ=255)
                                                                                                                                                                       gaGCGGCGATTTAAGTTTCTGGCACGCCTGGCTGGCGGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTgg                               <  2:671891/250‑1 (MQ=255)
                                                                                                                                                                         gcggcgAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGAt                             <  1:154821/251‑1 (MQ=255)
                                                                                                                                                                         gcggcgAGTTAAGATTCTGGCACGCCTGGCTGGCAGGGATTGTAGGTTGTATGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACACGATAAAACTGAACCTGCGATGCCTCAACACAGCCTGTTCACCCTTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGAGGCGGGaagggtggg                              >  2:920128/1‑243 (MQ=255)
                                                                                                                                                                          gggtgAGTTAAGTTTCTGGCACGCCTGGTGGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATc                            <  2:588197/247‑1 (MQ=255)
                                                                                                                                                                            gcgAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTg                          <  1:502289/251‑1 (MQ=255)
                                                                                                                                                                            gcgAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTg                          <  1:271930/251‑1 (MQ=255)
                                                                                                                                                                              gAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGcc                        <  2:225136/251‑1 (MQ=255)
                                                                                                                                                                               aGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAAGGCTGGAtattccg                       >  2:553288/1‑245 (MQ=255)
                                                                                                                                                                                         tCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAAt                                                <  2:55119/217‑1 (MQ=255)
                                                                                                                                                                                         tCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAAt                                                >  1:55119/1‑217 (MQ=255)
                                                                                                                                                                                             gCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGGTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAAttca          <  2:876822/251‑3 (MQ=255)
                                                                                                                                                                                               aCGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGGATGATGGGTTTGCCGGGGGccaacttcatt        >  2:266385/1‑241 (MQ=255)
                                                                                                                                                                                               aCGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATtcatt        >  1:799652/1‑251 (MQ=255)
                                                                                                                                                                                                 gCCTGGCTGGCAGGGCTTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATTCTGGATCTGCCGGGGGCAAAATTCATTAc      >  1:735934/1‑251 (MQ=255)
                                                                                                                                                                                                   ctggctggCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGGATGCTGGATCTGCCGGGGGCAAAATTCATTACGc    >  2:663762/1‑251 (MQ=255)
                                                                                                                                                                                                     ggctggcAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTAGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGACCGACGCGTCCGCTGGTGCCAATGGCGGCGGGAATGCTGGGTCTGCCGGTGGCAAAAttcaatccgccc  >  1:919264/1‑241 (MQ=255)
                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                       
CCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCG  >  NC_000913/71333‑71779

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