Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,786,936 C→T H33H (CAC→CAT envC → murein hydrolase activator EnvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,786,9360CT95.6% 206.6 / ‑6.3 68H33H (CAC→CATenvCmurein hydrolase activator EnvC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/37);  minor base G (1/2);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAG  >  NC_000913/3786690‑3787179
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
cgggAAAGCGGTTGCAGGCGGCAAACTTTCTGACAGCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCgg                                                                                                                                                                                                                                                 <  1:211660/248‑1 (MQ=255)
 gTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCgg                                                                                                                                                                                                                                                 >  1:217138/1‑250 (MQ=255)
   gAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATg                                                                                                                                                                                                                                              >  1:81551/1‑251 (MQ=255)
   gAAAGCGGGTGCAGGCGGCAACCTTTCTGACATCGCGCCTACCTTGTTTTCGCTGATGGGTATGGAAATCCCGCAAGAGATTACTGGTAAGTCGCGGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATCCCATGACACGGGCCGTGAAACCGCGCAGGTTTGTAATCAGGCCCATCATCTACGCCAGCGTGCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGTCCATGCGGATg                                                                                                                                                                                                                                              <  2:276338/251‑1 (MQ=255)
            tgcaGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGAc                                                                                                                                                                                                                                      <  2:407307/247‑1 (MQ=255)
               aGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCAGGCGGATGAGCGTGACCAAc                                                                                                                                                                                                                                  <  2:255655/251‑1 (MQ=255)
                  cggcAAACTTTCTGACATCGAGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCCTCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTc                                                                                                                                                                                                                                >  1:152920/1‑250 (MQ=255)
                   ggcAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACTGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCaa                                                                                                                                                                                                                              <  1:18440/251‑1 (MQ=255)
                      aaaCTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAATCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCAGCTACGCCAGCGTTCTGAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGTATGAGCGTGACCAACTCAAATc                                                                                                                                                                                                                           <  2:182047/251‑1 (MQ=255)
                       aaCTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAAGCCCGCAAGAGATGACTGGTAAGCCGCTGTTCAGCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCt                                                                                                                                                                                                                          <  2:216212/251‑1 (MQ=255)
                               gACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGc                                                                                                                                                                                                                  >  1:374610/1‑251 (MQ=255)
                                aCATCGCGCCGACCATGTTGTCGCTGATGGGGATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTGCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCAGGCGGATGAGCGTGATCAACTCAAATCTATTCAGGcc                                                                                                                                                                                                                 <  2:249538/251‑2 (MQ=255)
                                       gcCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATcg                                                                                                                                                                                                          >  2:14737/1‑251 (MQ=255)
                                                     cGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAgcg                                                                                                                                                                                             <  1:391933/250‑1 (MQ=255)
                                                         gATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCgg                                                                                                                                                                                        <  2:80870/251‑1 (MQ=255)
                                                          aTGGGTATGGAAACCCCGCAAGAGATGACGGGTAGGCCGCTGTTCATCGTAGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAAGACCATGCCACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCGTCATCGACGTCATCGTGCTTAGCTCTGGCGTAGTGTTGTGCGCGTTTTCCGCCCATGCGGATGAGTGTTACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGTGCGCGCGGt                                                                                                                                                                                       <  1:200199/251‑1 (MQ=255)
                                                          aTGGGTAGGGAAATCCCGCAAGAGATGACGGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGt                                                                                                                                                                                       <  2:138432/251‑1 (MQ=255)
                                                            gggTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAAGCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGTCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTAc                                                                                                                                                                                     <  2:292046/251‑1 (MQ=255)
                                                            gggTATGGAAAGCCCGAAAGAGAGGACTGGTAAGCCGATGTTCATCATGGAATAAGCCCTCCCCATGAGGGGAAAGGCTATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAATAGCGCGCGGTa                                                                                                                                                                                      <  2:17725/250‑1 (MQ=255)
                                                                aTGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGcc                                                                                                                                                                                  <  1:224347/250‑1 (MQ=255)
                                                                     aaTCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCCCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGc                                                                                                                                                                            <  2:121676/251‑1 (MQ=255)
                                                                          cGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGcaacaa                                                                                                                                                                       <  2:315424/251‑1 (MQ=255)
                                                                                  aTGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACgcgc                                                                                                                                                               <  1:136670/251‑1 (MQ=255)
                                                                                   tGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATcgccg                                                                                                                                                                                                       >  1:94382/1‑210 (MQ=255)
                                                                                   tGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATcgccg                                                                                                                                                                                                       <  2:94382/210‑1 (MQ=255)
                                                                                            aGCCGCTGTTCATCGTGGAATAATCCCGCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGGAATCAGGCCCATCATCTACGCCAGCGTTCGTAGCGCTGGCGTATGGTTGTGCGCCTTGTCCGCCCATACGGATGAGCGTGACCAAGTCAAAACTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTc                                                                                                                                                     <  1:59313/251‑1 (MQ=255)
                                                                                             gCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCg                                                                                                                                                    <  2:374610/251‑1 (MQ=255)
                                                                                                gCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGGGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCACCAACAACGCGCAAGCCTGCTCGccc                                                                                                                                                 >  2:370714/1‑249 (MQ=255)
                                                                                                 cTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCGGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCTCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCCGGACGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAAAAACCCGCAAGCCTGCGCGcac                                                                                                                                                 >  2:299451/1‑250 (MQ=255)
                                                                                                  tGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGcaca                                                                                                                                                >  1:386754/1‑250 (MQ=255)
                                                                                                   gTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCTCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAAt                                                                                                                                              <  1:374059/251‑1 (MQ=255)
                                                                                                      cATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTg                                                                                                                                            <  2:102222/250‑1 (MQ=255)
                                                                                                       aTCGTGGACTAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCAGGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGa                                                                                                                                           >  1:333125/1‑250 (MQ=255)
                                                                                                                      ctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGgaa                                                                                                                                >  1:173243/1‑246 (MQ=255)
                                                                                                                      ctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGc                                                                                                                           >  2:396847/1‑251 (MQ=255)
                                                                                                                      ctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGGGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGgaa                                                                                                                                <  2:173243/246‑1 (MQ=255)
                                                                                                                         cccATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAAGCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGAt                                                                                                                        <  1:132525/251‑1 (MQ=255)
                                                                                                                          ccATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGCTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCgctc                                                                                                                       >  2:94086/1‑248 (MQ=255)
                                                                                                                                  ggAAAGGCGATTAATACCATGACACGGGCCGTGTAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCACCCATGCGGATGAGCGTGTCCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACACCACACGCAACAACAACGCGCAAGCCTGCTCGTACAATTGAAACAGCAGGACGAAACGAtctatggagc                                                                                                               >  2:92776/1‑244 (MQ=255)
                                                                                                                                                     tGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTCGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGGGACCAACTCAAATCTATTAAGGCCGATATCGCCGCGAACGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAGTTGAAAAAGCAGGAAGAAGCGATCTCTGAAGGCACCCGAAAGcagagcga                                                                                             >  2:113508/1‑243 (MQ=255)
                                                                                                                                                        cacGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGAtc                                                                                                                       >  2:9174/1‑221 (MQ=255)
                                                                                                                                                        cacGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGAtc                                                                                                                       <  1:9174/221‑1 (MQ=255)
                                                                                                                                                           gggCCGTGAAACCGCGCAGGTTTGCCATCAGGACCATCATCTACGCCAGCGTTCTTAACGCTGGCGTATTGTTGTGCGCCATTTCCGCCCATGCGGATGAGAGTGACCCACTCAAATCTATTCAGGCCGATATCGCAGCGAACGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGACTGCTCGCACAATTGAAACCGCAGGAAGAAGCGATCTCTGAAGCCAACCGTAAGCTGCGCAAAACGCa                                                                                       >  2:407864/1‑250 (MQ=255)
                                                                                                                                                            ggCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGATCGCACAATTGAAAAAGCAGGAAGAAGCGCTCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCaa                                                                                      >  1:193804/1‑250 (MQ=255)
                                                                                                                                                              ccGTGAAACCTCGCAGGTGTGCAAGCAGGCCCATCATCTACGCCAGCGTTCCTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGTCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAac                                                                                   <  1:183375/251‑1 (MQ=255)
                                                                                                                                                                gTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAacac                                                                                 >  1:41286/1‑251 (MQ=255)
                                                                                                                                                                      ccGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCtcaa                                                                           <  1:14737/251‑1 (MQ=255)
                                                                                                                                                                       cgcgCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCtcaat                                                                          <  1:170884/251‑1 (MQ=255)
                                                                                                                                                                           cAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAAc                                                                      <  1:324395/251‑1 (MQ=255)
                                                                                                                                                                                   cAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGCTATCGCCGCGAAAGAGCGCGCGGTACGCCAACAGCAACAACAACGCGCAAGCCTGCTCGCACCATTGAAAAAGCAGGACGAAGc                                                                                                                           >  2:72917/1‑190 (MQ=255)
                                                                                                                                                                                   cAATCAGGCCCATCATCTACGCCAGCGTTCTGAGCGCTGGCGTAGTGTTGTGCGCCTTTGCCGCCCATGCGGATGAGCGTGACCAACGCAAACCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGc                                                                                                                           <  1:72917/190‑1 (MQ=255)
                                                                                                                                                                                                 atcTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGCAGCCACCCTTAAGCTGCGCGCAACGCAAACCACGCTCAAGACACTCAATAAACAGATtgatg                                                    >  2:416484/1‑247 (MQ=255)
                                                                                                                                                                                                     aCGCCAGCGTTCTTAGCGCTGGCGTATTGGTGTGCGCCTTTTCCTCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTgatgagatga                                              <  1:215227/249‑1 (MQ=255)
                                                                                                                                                                                                       gCCAGCGTTCTTAGCGCTGGCGTATTGGTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAAcgc                                          <  2:41286/251‑1 (MQ=255)
                                                                                                                                                                                                           gCGTTCTTAGCGCTGTCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAATACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCg                                      <  1:213902/251‑1 (MQ=255)
                                                                                                                                                                                                             gTTCTTAGCGCTGGCGTATTGTTGTGAGCCGTTTCCGCCCATGCGGATGAGCGGGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAGCACGCTCAATCAACTCAATAATCAGATTGATGAGATGAACGCGTCGAt                                    <  2:43035/251‑1 (MQ=255)
                                                                                                                                                                                                              ttCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGAtt                                   >  2:354771/1‑251 (MQ=255)
                                                                                                                                                                                                                  tAGCGCTGGCTTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGc                                 <  1:203756/249‑1 (MQ=255)
                                                                                                                                                                                                                     cgcTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGCACGCGTCGATTGCCAAAc                            >  1:99533/1‑251 (MQ=255)
                                                                                                                                                                                                                     cgcTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAAc                            >  2:13138/1‑251 (MQ=255)
                                                                                                                                                                                                                       cTGGCGAATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAAGCCGATATCGCCGCGACAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACACAACACACGCTCAATCAACTCAATAAACAGATTAATGAGATGAAcgcagcgaatgcaagcctg                          >  1:357750/1‑235 (MQ=255)
                                                                                                                                                                                                                         ggCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGTTGAGATGAACGCGTCGATTGCCAAAcggg                         >  2:272163/1‑247 (MQ=255)
                                                                                                                                                                                                                           cGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACACCAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGagc                      >  1:68920/1‑251 (MQ=255)
                                                                                                                                                                                                                            gTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGagc                      >  1:408784/1‑250 (MQ=255)
                                                                                                                                                                                                                              aTTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGagcagc                   >  1:146004/1‑251 (MQ=255)
                                                                                                                                                                                                                                            tttCCGCCCATGCGGATAAGCGTGACCAACACAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGAGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCATGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACTCTCAATCAACTCAATAAACAAATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAgccgcct       <  2:14253/249‑2 (MQ=255)
                                                                                                                                                                                                                                               ccGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGATCAGCAAAAAGCCGCCCAGGa   <  2:416327/250‑1 (MQ=255)
                                                                                                                                                                                                                                               ccGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAg  <  2:146004/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
CGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAG  >  NC_000913/3786690‑3787179

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: