Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,349 C→G A193A (GCG→GCC)  thiQ ← thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,3490CG98.5% 226.3 / NA 65A193A (GCG→GCC) thiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (32/32);  minor base T (0/1);  total (32/33)
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

CCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGC‑GACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTT  >  NC_000913/72061‑72640
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                   
ccACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAc                                                                                                                                                                                                                                                                                          <  1:78434/300‑1 (MQ=255)
  gcccgttgaggcatggttctctCGATACTGTGCGTAAATGGGTTTGGGGTTAGGGGATTGTCTTGTGCGGGTTGCCTGCTGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACACACTGTCCGTAGGTCGGATAAGTTGTTCACGCCACATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGTCTCACCGCTCACCCACTCATTGGTCTTACCCTTCCCTCCGATGCGATTGTCTGCTGAC‑ACCCTAGAGCGCGTGGCGATCCGCGCGGCATCTTCCAcac                                                                                                                                                                                                                                                                                        <  2:374492/279‑3 (MQ=255)
  attctttgattcctggggatagcGATATTTTGCGTAAAGTGGTTTGGGTTTAGTGGATTGGCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTGATCAGGCCTACAAAACGCACTGTCCGTAGTTCGGATGAGGTGTTCACGCCGCAGCGGACAGTGCATACTAACCTGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCCTACCCTGCCATCCGATGCACCCGTCCGC‑TACCAACACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAcac                                                                                                                                                                                                                                                                                        <  2:341821/278‑3 (MQ=255)
  actcGTTGATGCCGGTTTATATCGATATTTTGCTGAAAGTGGTTGGGGTTTAGTGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAGGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCATGCGATGCACCCGTCTGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAcac                                                                                                                                                                                                                                                                                        <  2:97038/297‑3 (MQ=255)
  aCCCGTTGATGCCGGTGTATATCGATATTTTGCTTAAAGTGGTTGTGGTTTAGTGGATTGTCTTGTGCGGGTTGCCTGAGGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAcac                                                                                                                                                                                                                                                                                        <  1:455870/300‑3 (MQ=255)
  aCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGTTTGGAGTTTATGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCTTGGCGATCCGCGCGGCATCTTCCAcac                                                                                                                                                                                                                                                                                        <  1:261874/300‑3 (MQ=255)
  aCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTTGGGTTTAGTGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCATGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAcac                                                                                                                                                                                                                                                                                        <  2:809021/300‑3 (MQ=255)
  aCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTTGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAcac                                                                                                                                                                                                                                                                                        <  1:432515/300‑3 (MQ=255)
  aCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAcac                                                                                                                                                                                                                                                                                        <  1:829707/300‑3 (MQ=255)
  aCCCGGTGATGCTGGTGTATATCGATATTTTGCGTAGAGTGGTTGGGGTTTAGGCGATTTGCTTATGCGGGTTGCCCGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCCTACCCTGCCATTCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAcac                                                                                                                                                                                                                                                                                        <  2:810058/300‑3 (MQ=255)
    ccGTGGATGCCGGTGTATATCGATATTTTGCGTAAAGTGTTTGGGGTTTAGTGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTTAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCATGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCGTCCAcacgg                                                                                                                                                                                                                                                                                      <  1:852203/300‑5 (MQ=255)
          aTGCCGGTGTATATCGATATTTTGCGTAAAGTGTTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGac                                                                                                                                                                                                                                                                                <  1:679762/300‑1 (MQ=255)
          aTGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGTGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGac                                                                                                                                                                                                                                                                                <  2:548768/300‑1 (MQ=255)
          aTGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCACCAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACTACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACCCTGTGTGac                                                                                                                                                                                                                                                                                >  2:251278/1‑300 (MQ=255)
            tccGGTGTATATCGATATTTTGCGTAAAGTGTTTTGGGTTTATGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGGCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCACTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                              <  1:198783/299‑1 (MQ=255)
            tccGGTGTATATCGATATTTTGCGTAAAGTGGTTGGTGTTTAGGGGATTGTCTTGTGCGGGTGGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTNCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTTAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCATGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACtgtttgacac                                                                                                                                                                                                                                                                              <  2:210436/299‑8 (MQ=255)
            gCCGGTGTGTATCGGTATTTGGCGTAAAGTGGTTGGGGTTTAGGGGATTGTTTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTTCGTAGGTCGGATAAGGTGTTCACGCCGCATTTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTTTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                              <  1:244703/300‑1 (MQ=255)
            gCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                              <  1:127197/300‑1 (MQ=255)
            gCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACNGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                              <  1:574233/300‑1 (MQ=255)
            gCCGGTGTATATCGATATTTTGCGTAAAGGGGGTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAATCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                              <  1:208643/300‑1 (MQ=255)
               ggTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCAt                                                                                                                                                                                                                                                                           >  1:640449/1‑300 (MQ=255)
                       tCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTAATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTACACACTGTGTGACACCATCAATaacg                                                                                                                                                                                                                                                                   >  2:420854/1‑297 (MQ=255)
                                 gCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTtg                                                                                                                                                                                                                                                         >  1:644401/1‑300 (MQ=255)
                                               gggtttaggggcttatCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAAATCATTGGTCTTACCCTGCCAGGCGATGCGCACGTCGGC‑GACCACCACCGCGCGCGTGGCGATCCGCGCGGTATTTTCCACACTGTGTGAAAACATCAATAACGTCATTTTTTGCTGCTGGCAGCTCg                                                                                                                                                                                                                                           >  2:498481/16‑300 (MQ=255)
                                               gggTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGTGTGGCATCTTCCACACTGTGTGAGACCATTAATAAGGTTATTCTTTGCTGCTGGCAgatca                                                                                                                                                                                                                                           >  2:141256/1‑296 (MQ=255)
                                               gggTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCCGACAGTGCATACTAAACCGTAATCCACAATAATGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGACAGGCGATGCGCCCGTCGGC‑GCCCACCCACGAGTGCGTGGCGATCCGCGCGGCATTTTCCACATTGTGTGACACCATCAATAATGCCATTTTTTGCTGCTGGCAGCttg                                                                                                                                                                                                                                           >  2:142948/1‑298 (MQ=255)
                                                             gTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTc                                                                                                                                                                                                                             >  1:369136/1‑300 (MQ=255)
                                                                  gtgCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCAt                                                                                                                                                                                                                        <  1:739381/300‑1 (MQ=255)
                                                                          tGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGAc                                                                                                                                                                                                                >  1:415729/1‑300 (MQ=255)
                                                                           gcctaattcGACGCTGACCGCTTCTTAGCAGGCCGATAAAACGCACTGTCCTTAGGTTGGATAAGTTGTTCACGCCGCCTCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTTAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCCTACCCTTCCATGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACg                                                                                                                                                                                                               <  2:237294/292‑1 (MQ=255)
                                                                                 tGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGc                                                                                                                                                                                                         <  1:763132/300‑1 (MQ=255)
                                                                                          gaccGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACAATGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAg                                                                                                                                                                                                >  2:658328/3‑300 (MQ=255)
                                                                                                    tATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAAc                                                                                                                                                                                      <  1:381650/300‑1 (MQ=255)
                                                                                                         ggCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACAACATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTc                                                                                                                                                                                 >  2:62515/1‑300 (MQ=255)
                                                                                                                        ctgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTTCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAAt                                                                                                                                                                  >  1:498938/5‑300 (MQ=255)
                                                                                                                          gtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCg                                                                                                                                                                >  1:844637/3‑300 (MQ=255)
                                                                                                                          gtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCg                                                                                                                                                                >  1:550074/3‑300 (MQ=255)
                                                                                                                                     ggATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAATGGTTCATCGAGCAATAAAATCGGCTGTTcgcgc                                                                                                                                                     >  1:138864/1‑299 (MQ=255)
                                                                                                                                               ttCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAAc                                                                                                                                           <  1:408580/300‑1 (MQ=255)
                                                                                                                                               ttCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCACCCGTCGGC‑GACCACCATCGAGCGCGNGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTACTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAAc                                                                                                                                           <  2:211606/300‑1 (MQ=255)
                                                                                                                                                         atctgACAGTGCATACTAACCCGTAATCCCCAATAGTGCTTAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGc                                                                                                                                 <  2:821191/296‑1 (MQ=255)
                                                                                                                                                            tgACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCNCACCAGCGTCAGCATCTCCTGACGTAACGACGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCAc                                                                                                                              >  2:755259/2‑300 (MQ=255)
                                                                                                                                                                aGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGGCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACgcgc                                                                                                                          <  2:19319/300‑1 (MQ=255)
                                                                                                                                                                 gTGCATACTAACCCGTCATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACCGTGTGAAACCATCAATAACGCCATTTTTTGCTGCTGGCAGCTCGTGCTTAACAGCGTCAGCATCTACTGACGTAAAGCCGGATCGAGCGCAGAGAATGGTTCATCGAGCAATGAAGTCGGCTGTTCGCGCACCAGACAACGAGCTAACGCCACgcgcg                                                                                                                         >  2:486295/1‑299 (MQ=255)
                                                                                                                                                                   gCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCACCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCCACTGCCACgcgctgt                                                                                                                       >  2:355947/1‑300 (MQ=255)
                                                                                                                                                                     aTACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCGATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGANCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACgcnctgtcg                                                                                                                     >  2:714058/1‑300 (MQ=255)
                                                                                                                                                                     aTACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACgcgctgtcg                                                                                                                     >  1:98196/1‑300 (MQ=255)
                                                                                                                                                                            cccGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGAccg                                                                                                              >  1:805243/1‑300 (MQ=255)
                                                                                                                                                                                     cccAATAATGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGTCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGc                                                                                                     <  2:40667/300‑1 (MQ=255)
                                                                                                                                                                                                   aacACTCGCCTTACCGCTAAACAACTCATTGGTCTTACCCTGCCAGGCGAGGCGCCCGTCTGC‑GACCACCACCGAGCGCGTGTCGATCCGCGCGGCGTCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCGGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAAGAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACcg                                                                                       <  2:44780/298‑1 (MQ=255)
                                                                                                                                                                                                   aGCACTCGCCTTACCACTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCACCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACcg                                                                                       <  2:386070/300‑1 (MQ=255)
                                                                                                                                                                                                           ccTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTa                                                                               >  1:782453/1‑300 (MQ=255)
                                                                                                                                                                                                                          acTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATTTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCAACGCGCTGTCGCTGACCGCCGGAAAGCTTGCCAGGTAACCGCGCAATTAAATTATCAATCCCCa                                                                >  2:270357/1‑300 (MQ=255)
                                                                                                                                                                                                                            tCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATc                                                              >  1:134774/1‑300 (MQ=255)
                                                                                                                                                                                                                                 ggTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTggtg                                                         >  2:764490/1‑298 (MQ=255)
                                                                                                                                                                                                                                 ggTCTTACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTggcg                                                         >  2:764592/1‑300 (MQ=255)
                                                                                                                                                                                                                                  gtcttACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTggcgg                                                        >  1:315009/5‑300 (MQ=255)
                                                                                                                                                                                                                                  gtcttACCCTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCNGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGATGTCGCTGACCGCCGGAAAGCTCGCAAGGTAACCGTGCAATTAAATTATCAATCACAACCTggcgg                                                        >  2:204333/5‑300 (MQ=255)
                                                                                                                                                                                                                                  gtcttACCCTGCCAGGCGATGCGCCCGTCAGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATCTTTTGCTGCTGGCAGCTCGTGCTTACCAGCGTCAGCATCTCCTGACGTAACGCAGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTATTCGCGCACTCGACAACGCGCTAACGCCACGTGCTGTCGCTGATCGCATGGAAGCTCGACAGGTAACTGCGCCATTCAATTATCAATTCCCATTTggctg                                                        >  2:56244/5‑298 (MQ=255)
                                                                                                                                                                                                                                         ccTGCCATGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTTGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAg                                                 <  1:240626/300‑1 (MQ=255)
                                                                                                                                                                                                                                         ccTGCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACTTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAg                                                 <  1:240988/300‑1 (MQ=255)
                                                                                                                                                                                                                                            gCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGACGGAAAGCTCGCCAAGTAACCGCGCCACTAAATTATCAATCCCCATCTGACGGGCGCTAGCGt                                              >  2:139214/1‑300 (MQ=255)
                                                                                                                                                                                                                                            gCCAGGCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGAAACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTNAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGNCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCCGTCGCTGACCGCCGGAAAGCTCGCCACGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGt                                              >  2:682777/1‑300 (MQ=255)
                                                                                                                                                                                                                                                ggCGATGCGCCCGTCGGC‑GACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCAt                                          >  1:727579/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                         aTCCGCGCTGCATCTTCCACACTGTTTGACACCATCAATAACGTCATTTTTTGCTGCTTGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACTTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCNATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGTGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGAttt  <  2:194985/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                   
CCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGC‑GACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTT  >  NC_000913/72061‑72640

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: ‑