Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,977,914 A→C M95R (ATG→AGG)  flhD ← DNA‑binding transcriptional dual regulator FlhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,977,9140AC96.0% 268.0 / ‑0.9 75M95R (ATG→AGG) flhDDNA‑binding transcriptional dual regulator FlhD
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (27/45);  total (29/46)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

CATGCCTTTCGGCGGTGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTG  >  NC_000913/1977677‑1978159
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
caTGCCTTTCGGCGGTGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATg                                                                                                                                                                                                                                           <  1:52892/251‑1 (MQ=255)
 aTGCCTTTCGGCGGTGGGCTTCCGCGCAGTTCTTTATAAAGGTTTATCAGGCGTCCGCGTCTTAACTGTGATTAGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGGGATCCAATCCATTGCCCGCTGAATATCCCGCGCTTCCTGGACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATg                                                                                                                                                                                                                                           >  1:128968/1‑250 (MQ=255)
    ccTTTCGGCGGTGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAAtt                                                                                                                                                                                                                                       <  1:154160/251‑1 (MQ=255)
            cggTGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGgag                                                                                                                                                                                                                               <  2:298567/251‑1 (MQ=255)
                      cctcgcAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGGGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAAc                                                                                                                                                                                                                      <  2:802455/247‑1 (MQ=255)
                      ccGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAAcg                                                                                                                                                                                                                     <  1:422263/251‑1 (MQ=255)
                        gcgcAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAAcgcg                                                                                                                                                                                                                   <  1:872043/251‑1 (MQ=255)
                        gcgcAGTTCTTTAGAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCCTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTTAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAAcgcg                                                                                                                                                                                                                   <  2:356396/251‑1 (MQ=255)
                         cgcAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCgg                                                                                                                                                                                                                  <  1:281458/251‑1 (MQ=255)
                             gTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGAGCGTCAACGCGGGAAt                                                                                                                                                                                                              <  1:129288/251‑1 (MQ=255)
                                  ttATAAAGTTTTATCAGGCGGCCGCGACTTAACTTTGTTTCGCCTTCCAGCAACCGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCACTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTAGTAAATTGCTGGAGATCGTCAACGCGGGAATCTTg                                                                                                                                                                                                          <  2:128968/250‑1 (MQ=255)
                                   tataAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTGCTTCAGGCTGATTAACAGCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGc                                                                                                                                                                                                         <  1:103418/250‑1 (MQ=255)
                                         gTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACt                                                                                                                                                                                                  >  1:53133/1‑251 (MQ=255)
                                          ttttATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGAGTGTGGAGAGCCTGCTGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGGGTCAACTg                                                                                                                                                                                                 >  2:461520/1‑251 (MQ=255)
                                               tCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTa                                                                                                                                                                                             >  1:784385/1‑250 (MQ=255)
                                                        cgcgACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCTGGGTGATCAATTCCATTGCAAGATGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCACGCGTGTTGAGAGCATGATGCAGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATAGTCtggt                                                                                                                                                                                   >  2:550117/1‑251 (MQ=255)
                                                          cgACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCtggtgg                                                                                                                                                                                 >  2:10407/1‑251 (MQ=255)
                                                               tAACTGTGTTTCGCGTTCCCGCATCTGCAAACGAGCGCCCAGGGTGAGTAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAAGGCTTTTTTCACGCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAATAATCGTCTGGTGGCTg                                                                                                                                                                              <  2:801323/249‑1 (MQ=255)
                                                                 aCTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGCCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCaa                                                                                                                                                                          <  1:725774/251‑1 (MQ=255)
                                                                  cTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCaa                                                                                                                                                                          >  1:35990/1‑250 (MQ=255)
                                                                    gtgtTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAAc                                                                                                                                                                       <  1:386936/251‑1 (MQ=255)
                                                                      gtTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACgg                                                                                                                                                                     >  1:764491/1‑251 (MQ=255)
                                                                         tcgatttCCAGCATCGGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCAGTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAAGCGTGCGTCAAATGAGTAATCGTCTGGTGGCTGTCAAAACGGaa                                                                                                                                                                   <  2:183944/246‑1 (MQ=255)
                                                                         tCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCAGTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTGTCGTCATCTGAGGAATCGTCTGGTGGCTGTCAAAACTGaa                                                                                                                                                                   <  1:200719/250‑1 (MQ=255)
                                                                          cGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGt                                                                                                                                                                 <  1:272016/251‑1 (MQ=255)
                                                                                     atctgcgcaCGAGCGCCCAGGGTGATCTATTCCATTGCCAGCTGAATATCCAGCGCTTACTGAACAATGCTTTTTTCACTCATGATCAGTCCCTTTTCTTGCGCCGCGCTTCTTCAGGCTGCTTAACATCCTTCAGCAAGCGTGTTGCGAGCCTGATGCAGGTATGACTTTGCTGGAGATCGTCCACGAGGGACTCATGCGTCAACTGAGTACTCGACTGGTGTCTGTCACCACGACAGTAACAAACCAGt                                                                                                                                                      >  2:858896/9‑251 (MQ=255)
                                                                                              cGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGtt                                                                                                                                            <  1:46945/251‑1 (MQ=255)
                                                                                                   gcCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGc                                                                                                                                       <  1:310087/251‑1 (MQ=255)
                                                                                                       aGGGTGATCAATTCCATTTCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTGTCACTCATGATCAGTCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGATATCGTCAACGCGTGAATCTTGCGTCAATTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGc                                                                                                                                   <  1:572574/251‑1 (MQ=255)
                                                                                                       aGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAACCCAGTTGATTGGTTTCTGCCAGc                                                                                                                                   >  1:413976/1‑251 (MQ=255)
                                                                                                         ggagATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCtt                                                                                                                                 <  2:728449/248‑1 (MQ=255)
                                                                                                         ggTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCtt                                                                                                                                 <  1:322747/251‑1 (MQ=255)
                                                                                                            taTCAATTACATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGACTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGGCTGGTGGCTGTCAAAACGAAAGTGTAAAATCAGTTATGTGGTTTCTGCCatctaaac                                                                                                                              >  2:582797/2‑244 (MQ=255)
                                                                                                               cAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCa                                                                                                                            <  1:874264/250‑1 (MQ=255)
                                                                                                                          ccAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAag                                                                                                                  >  1:226820/1‑249 (MQ=255)
                                                                                                                              cTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAAc                                                                                                                                                                                                                      <  1:236621/146‑1 (MQ=255)
                                                                                                                              cTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAAc                                                                                                                                                                                                                      >  2:236621/1‑146 (MQ=255)
                                                                                                                                gAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAAAGTCAGTg                                                                                                          >  2:465764/1‑251 (MQ=255)
                                                                                                                                      cccGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTa                                                                                                    <  1:648630/251‑1 (MQ=255)
                                                                                                                                              tCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTg                                                                                              <  1:717129/249‑1 (MQ=255)
                                                                                                                                              tCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTc                                                                                            >  1:727260/1‑251 (MQ=255)
                                                                                                                                                cTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGc                                                                                          >  1:278229/1‑251 (MQ=255)
                                                                                                                                                  gAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGcc                                                                                         <  2:764491/250‑1 (MQ=255)
                                                                                                                                                  gAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGcc                                                                                         >  1:135323/1‑250 (MQ=255)
                                                                                                                                                        tGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCt                                                                                                                                                                                      <  2:171517/152‑1 (MQ=255)
                                                                                                                                                        tGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCt                                                                                                                                                                                      >  1:171517/1‑152 (MQ=255)
                                                                                                                                                            tttttCTCTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCa                                                                                <  1:436862/249‑1 (MQ=255)
                                                                                                                                                                tCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGGCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTAt                                                                           <  2:784385/250‑1 (MQ=255)
                                                                                                                                                                 cACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATg                                                                          <  1:134687/250‑1 (MQ=255)
                                                                                                                                                                 cACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGc                                                                         <  2:141363/251‑1 (MQ=255)
                                                                                                                                                                   cTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCg                                                                       >  1:324319/1‑251 (MQ=255)
                                                                                                                                                                        gATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGGGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGCCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCgaac                                                                    >  2:665925/1‑247 (MQ=255)
                                                                                                                                                                           cAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGACATTTCTTCATTTATGCCGCGACgcac                                                               >  2:603561/1‑248 (MQ=255)
                                                                                                                                                                            aGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAAc                                                              >  2:24914/1‑251 (MQ=255)
                                                                                                                                                                               cccTTGTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATa                                                           <  2:384545/251‑1 (MQ=255)
                                                                                                                                                                                ccTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGCAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATa                                                           >  2:319153/1‑250 (MQ=255)
                                                                                                                                                                                     tCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCACGCGTGTTGAGAGCCTGCTGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACAATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCCCCATTTCTTCATTTATGCCGAGACGAAACATAGCggcc                                                     >  2:31294/1‑249 (MQ=255)
                                                                                                                                                                                        tGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGAc                                                     <  2:752229/248‑1 (MQ=255)
                                                                                                                                                                                            cAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGt                                              <  1:298439/251‑1 (MQ=255)
                                                                                                                                                                                                    cttcAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAAc                                       >  1:468316/1‑250 (MQ=255)
                                                                                                                                                                                                      tcAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAAc                                       <  2:226820/248‑1 (MQ=255)
                                                                                                                                                                                                      tcAGGCTGATTAACATCATTCAGCAAGCGTGATGAGAGCCTGATTCCGGTATGAAGTTGCTGGAGTGCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAAc                                       <  2:312632/248‑1 (MQ=255)
                                                                                                                                                                                                       cAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATc                                   >  1:307107/1‑251 (MQ=255)
                                                                                                                                                                                                           cTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAAcc                                                                                                                                                         >  1:185615/1‑130 (MQ=255)
                                                                                                                                                                                                           cTGATTAACATCATTCAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAAcc                                                                                                                                                         <  2:185615/130‑1 (MQ=255)
                                                                                                                                                                                                                      catTCAGCAAGAGTGATGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGGCGCTTTGTCCTGAACAATCAACCGCTGTGCAagt                    >  2:500914/1‑251 (MQ=255)
                                                                                                                                                                                                                          cAGCAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAgtagta                <  2:324319/251‑1 (MQ=255)
                                                                                                                                                                                                                             cAAGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAat             <  1:675869/251‑1 (MQ=255)
                                                                                                                                                                                                                              aaGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTAAGTGCCGCAAACGTTGGCGCCATTTATTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTAAACAATAAAAAGATGGGCAAGTAGTaaca             >  2:463394/1‑248 (MQ=255)
                                                                                                                                                                                                                               aGCGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGAATGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAatat           <  2:53133/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTGTTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGa         <  1:533485/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gtTGAGAGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACaa      <  1:765557/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gagaGCCTGATGTCGTTATGAATTTGCTGGTGATCGTCAACGCGGGAATCTTGCGTCAACTTAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTTCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCGTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGt    <  2:413370/250‑1 (MQ=255)
                                                                                                                                                                                                                                       gagaGCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGtt   <  1:603561/251‑1 (MQ=255)
                                                                                                                                                                                                                                        agagCCTGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACAGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTAACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCAGTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTg  <  1:858896/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
CATGCCTTTCGGCGGTGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTG  >  NC_000913/1977677‑1978159

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