Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,786,984 A→G A49A (GCA→GCG envC → murein hydrolase activator EnvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,786,9840AG98.4% 207.2 / ‑6.1 63A49A (GCA→GCGenvCmurein hydrolase activator EnvC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (29/33);  total (29/34)
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.68e-01

CGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAGCGCAGCCTCGCCGCACAACTGGATGCCGCATTCCGTCAGGGCGAGCATA  >  NC_000913/3786743‑3787228
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
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)
                                       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)
                                                                                                                                                                                              ccatgcgGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAgcgcc                                               >  1:74732/5‑250 (MQ=255)
                                                                                                                                                                                              ccatgcgGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAGCGCa                                               <  1:298598/247‑1 (MQ=255)
                                                                                                                                                                                              ccatgcgGATGAGCGTGAACAACTCACATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACACTTGCAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAgcgc                                                >  2:185236/5‑250 (MQ=255)
                                                                                                                                                                                                      aTGAGCGTGACCAACTCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGACGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTAAATCAACTCCATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAGCGCAGCCTcgcc                                       >  1:67615/1‑251 (MQ=255)
                                                                                                                                                                                                          gCGTGCCCAACTCAAATCTATTCCGGCAGATATCGCCGCGAAAAAGCGCGCGGTACGCAAAAAGCAACAACAACGCGAAAGCCTGCTCGCACAATTGAAAACGCAGGAAGAAGCCATCTCTGAAGCCCCCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACACAAAAAACAGATTGCTGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCACAGGAGCGCAGCCTCGCCGcac                                   >  2:185295/1‑251 (MQ=255)
                                                                                                                                                                                                                     tCAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGAGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAGCGCAGCCTCGCCGCACAACTGGATGcc                        >  2:344913/1‑251 (MQ=255)
                                                                                                                                                                                                                      cAAATCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCACAAAGCCGCCCAGGAGCGCAGCCTCGCCGCACAACTGGATgtc                        <  2:386754/250‑3 (MQ=255)
                                                                                                                                                                                                                          tCTATTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAGCGCAGCCTCGCCGCACAACTGGATGCCGCAtt                   >  1:91298/1‑251 (MQ=255)
                                                                                                                                                                                                                             aTTCAGGCCGATATCGCCGCGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAGCGCAGCCTCGCCGCACAACTGGATGCCGCATTc                  <  1:332185/249‑1 (MQ=255)
                                                                                                                                                                                                                                           cgccgcGAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAGCGCAGCCTCGCCGCACATCTGGATGCCGCATTCCGTCAGGGCGAACATa  <  1:370714/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
CGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAGCGCAGCCTCGCCGCACAACTGGATGCCGCATTCCGTCAGGGCGAGCATA  >  NC_000913/3786743‑3787228

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