Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,786,972 T→C A45A (GCT→GCC envC → murein hydrolase activator EnvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,786,9720TC98.5% 214.5 / NA 65A45A (GCT→GCCenvCmurein hydrolase activator EnvC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (30/34);  minor base A (1/0);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAGCGCAGCCTCGCCGCACAACTGGATGCCGCATTC  >  NC_000913/3786729‑3787212
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
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                                                                                                                                                                                                                                        <  2:94382/210‑1 (MQ=255)
                                            tGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCCGATATcgccg                                                                                                                                                                                                                                        >  1:94382/1‑210 (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)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
GCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAACGCAAAACACGCTCAATCAACTCAATAAACAGATTGATGAGATGAACGCGTCGATTGCCAAACTGGAGCAGCAAAAAGCCGCCCAGGAGCGCAGCCTCGCCGCACAACTGGATGCCGCATTC  >  NC_000913/3786729‑3787212

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: