Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,346,592 C→T A51T (GCT→ACT)  pdeR ← cyclic di‑GMP phosphodiesterase PdeR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,346,5920CT98.4% 214.7 / NA 63A51T (GCT→ACT) pdeRcyclic di‑GMP phosphodiesterase PdeR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (40/22);  total (41/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTT  >  NC_000913/1346296‑1346882
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                  
gtgtTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTg                                                                                                                                                                                                                                                                                                 >  2:373218/1‑300 (MQ=255)
     cTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGATGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCgg                                                                                                                                                                                                                                                                                            >  2:206174/1‑300 (MQ=255)
      tttCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGc                                                                                                                                                                                                                                                                                           <  1:517304/300‑1 (MQ=255)
        tCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAg                                                                                                                                                                                                                                                                                         <  2:265712/300‑1 (MQ=255)
             ccTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAAc                                                                                                                                                                                                                                                                                    >  2:339063/1‑300 (MQ=255)
              ctatgtAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACt                                                                                                                                                                                                                                                                                   >  2:122432/4‑300 (MQ=255)
               tgtgTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTNCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTg                                                                                                                                                                                                                                                                                  >  1:464767/1‑300 (MQ=255)
                gtgtAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACNAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGc                                                                                                                                                                                                                                                                                 >  2:597877/1‑300 (MQ=255)
                gtgtAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGTAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGc                                                                                                                                                                                                                                                                                 <  1:152501/300‑1 (MQ=255)
                gtgtAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAAccgc                                                                                                                                                                                                                                                                                 >  1:617603/1‑297 (MQ=255)
                gtgtAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGc                                                                                                                                                                                                                                                                                 <  1:152132/300‑1 (MQ=255)
                         tCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTNTTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCATGATGCGAGCAGTCTGTTCGGCAGATAACCGCAATGTACGa                                                                                                                                                                                                                                                                        >  2:165389/1‑300 (MQ=255)
                          cacaTAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGAt                                                                                                                                                                                                                                                                       <  1:407408/300‑1 (MQ=255)
                                          gCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTc                                                                                                                                                                                                                                                       >  1:254397/1‑300 (MQ=255)
                                                 atatTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTAAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTa                                                                                                                                                                                                                                                >  1:610270/1‑300 (MQ=255)
                                                               gttgagaATAACTGTTGCGGAATGTGTTTTAGAAACTACCTGCTCAGCAAAGGGAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATACAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCATGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCaaaat                                                                                                                                                                                                                                  <  2:308168/297‑2 (MQ=255)
                                                                          cTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTAACAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGATGCttt                                                                                                                                                                                                                       >  2:433767/1‑300 (MQ=255)
                                                                              ggctgAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGGGTGATAAATCACGAACAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGACCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTctgt                                                                                                                                                                                                                   <  2:26489/296‑4 (MQ=255)
                                                                                               aaattACCTGCGCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGGCATGCCATGCAGACACGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCGACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGAt                                                                                                                                                                                                  <  1:39565/296‑1 (MQ=255)
                                                                                               aaaCTACTTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGAt                                                                                                                                                                                                  <  2:531891/300‑1 (MQ=255)
                                                                                                      cTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGCGTTTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCNACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCCTCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGa                                                                                                                                                                                           <  2:598309/300‑1 (MQ=255)
                                                                                                        gCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTCCGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGTGATc                                                                                                                                                                                         >  2:66470/1‑300 (MQ=255)
                                                                                                                 aGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGt                                                                                                                                                                                >  1:308703/1‑300 (MQ=255)
                                                                                                                      aGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTAc                                                                                                                                                                           <  2:396328/300‑1 (MQ=255)
                                                                                                                        cccGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACaa                                                                                                                                                                         <  2:415974/300‑1 (MQ=255)
                                                                                                                           gtttgATAAACCACGGGCAACCGCCGTTGTGTCATGCCAGGCAGACGCGTTGCCAGCCCATTCCCGTTGATTGATTCTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGTTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGGTGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGt                                                                                                                                                                      <  2:91343/297‑1 (MQ=255)
                                                                                                                           gtgtGATAAATCACGAGCAACCGCCGGTGTGTCATGCCAGGCAGACGCGTTGCCCGCCCATTCCCGTTTATTGATTTTTCGTTCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGGCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGAGACGAGAAAGTTGTACAACGt                                                                                                                                                                      <  2:594435/300‑1 (MQ=255)
                                                                                                                                aTAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTgg                                                                                                                                                                 <  1:61358/300‑1 (MQ=255)
                                                                                                                                  aaaccACGAGCAACCGCCGGTGTGTCATGCCAGGCAGACGCGTTGCCACCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTNAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGNAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGAc                                                                                                                                                               <  2:547912/296‑1 (MQ=255)
                                                                                                                                  aaaTCACGAGCAACCGCCGTTGTGTCATACCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTTGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCATACTCATCATCAGGCTGGAGTTGACGACCTTCATTTCCCTGATGCGAGCATTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGAATTGTGCGATCCGAGAAAGTTGGACAACGTTGTGGAc                                                                                                                                                               <  2:34598/300‑1 (MQ=255)
                                                                                                                                            cAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACgg                                                                                                                                                     <  1:516344/300‑1 (MQ=255)
                                                                                                                                                  ccgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGCCGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTACCTGATGCGCGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAACGCGCAAAACATCGCTGCTTTCAGTAAACCGACAGTATGGATCGTGCGATCCGCGAAAGTTGTCCAACGTTGTGGACTCCCTAACGGTCTTCa                                                                                                                                               >  2:9468/1‑300 (MQ=255)
                                                                                                                                                  ccgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCa                                                                                                                                               >  1:17735/1‑300 (MQ=255)
                                                                                                                                                     gtgtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCatga                                                                                                                                            >  2:127237/1‑300 (MQ=255)
                                                                                                                                                       gtgtCATGCCATGCAGACGCGTTGCCAGCCCCTCCCCGTTTATTGATTTTTCGTCCTCCCAGATGCACAGAAAGATCGCTTTCACCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGCATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTACCGATTTTCctgatg                                                                                                                                          >  1:592317/1‑300 (MQ=255)
                                                                                                                                                       gtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATTGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAACGCGCAAAACATCGCTGCTTTCCGTCACCCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGCCTCACTAACGGTTTTCatgatg                                                                                                                                          >  2:256479/1‑300 (MQ=255)
                                                                                                                                                       gtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCatgatg                                                                                                                                          >  1:664433/1‑300 (MQ=255)
                                                                                                                                                       gtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGCTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGTGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAAAGGTTTTCatgatg                                                                                                                                          >  2:355520/1‑300 (MQ=255)
                                                                                                                                                       gtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCGACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCatgatg                                                                                                                                          >  2:624651/1‑300 (MQ=255)
                                                                                                                                                         gtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCa                                                                                                                                        >  2:240307/1‑300 (MQ=255)
                                                                                                                                                         gtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCa                                                                                                                                        >  1:380607/1‑300 (MQ=255)
                                                                                                                                                         gtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTAAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCa                                                                                                                                        >  2:271704/1‑300 (MQ=255)
                                                                                                                                                         gtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTAACAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGATGCTTTCCGTCAACCGCCAGTATGGATTGTGTGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCa                                                                                                                                        >  2:336938/1‑300 (MQ=255)
                                                                                                                                                             tgccatgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGTTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtc                                                                                                                                    <  1:639073/300‑1 (MQ=255)
                                                                                                                                                                        cgcgTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTACGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAAc                                                                                                                         >  2:268255/1‑300 (MQ=255)
                                                                                                                                                                                  gCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACCGAAAGATCGCTTTCNTCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTNGTGGACTCCATCACGGTTTTCATGATGCATATCCTCAGGTTAACAGCGGAGTcc                                                                                                               >  2:644068/1‑300 (MQ=255)
                                                                                                                                                                                  gCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTcc                                                                                                               >  1:71776/1‑300 (MQ=255)
                                                                                                                                                                                       ttCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGt                                                                                                          <  1:343071/300‑1 (MQ=255)
                                                                                                                                                                                                           gTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCAGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATg                                                                                      >  1:512859/1‑300 (MQ=255)
                                                                                                                                                                                                              cTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACt                                                                                   >  1:518389/1‑300 (MQ=255)
                                                                                                                                                                                                                    gATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATg                                                                             <  1:597346/300‑1 (MQ=255)
                                                                                                                                                                                                                          caGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTa                                                                       >  1:278734/1‑300 (MQ=255)
                                                                                                                                                                                                                                               cGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTCACGGCGGCTAAACGCATTGtt                                                  >  2:409824/1‑300 (MQ=255)
                                                                                                                                                                                                                                               cGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtt                                                  >  1:376010/1‑300 (MQ=255)
                                                                                                                                                                                                                                               cGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGATGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTNAATAAAAGGAACAACATGACTTACATGAAATTAGCGGCGGCTAAACGNATTGtt                                                  >  2:576217/1‑300 (MQ=255)
                                                                                                                                                                                                                                                 gTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATNGCTGCTTTCCGTCAACCGCCAGTATGNATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGTATTGtttt                                                >  2:465214/1‑300 (MQ=255)
                                                                                                                                                                                                                                                     gacgcgtcgtcaGGCTGGAGGTGGTGATCGTCATTTCCCTGATTCTAGCAGTCTGTTTGGTAGATAGCTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAGCATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCAATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTGACATGAAATGAACGGCGGCTAAACGCATTGtttttttt                                            <  2:42280/291‑1 (MQ=255)
                                                                                                                                                                                                                                                     gACTCATCATCAGGCTGGAGGTGATGACCGTNATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTNCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttt                                            >  2:490176/1‑300 (MQ=255)
                                                                                                                                                                                                                                                         catcatcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCGACCGCAAGTATGGATTGTGCGATCCGAGAAAGTTGTACCACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGATTAAATAAACAGAACAACATGACTTACATGAAATTAACGGCGGCTAAAAGCATTGTTTTTTTtatat                                        >  2:146655/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                 ggCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATatttattt                                >  1:388376/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                        tCCCTGATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGACAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGANTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGc         >  2:610902/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                             gATGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGACAGtt    >  2:521982/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                               tGCGAGCAGTCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCCACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGtttt  >  1:536557/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                  
GTGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTT  >  NC_000913/1346296‑1346882

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