Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,634 G→T G276G (GGG→GGT)  yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,6340GT95.1% 189.0 / ‑3.3 63G276G (GGG→GGT) yaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (28/30);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCT  >  NC_000913/46389‑46878
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
gAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTgg                                                                                                                                                                                                                                                   >  1:90185/1‑251 (MQ=255)
           gtgggtattaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACGGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCgg                                                                                                                                                                                                                                        <  2:198851/246‑1 (MQ=255)
           cTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGGTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTTCCAGGTGATCCCAATGTTCGCCATTTACGCCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCgg                                                                                                                                                                                                                                        >  1:545222/1‑251 (MQ=255)
             ggTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGcg                                                                                                                                                                                                                                      <  1:27274/251‑2 (MQ=255)
                 attaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTTGGTGTTGGCAAAAACGCGGcgctg                                                                                                                                                                                                                                   <  1:243356/250‑5 (MQ=255)
                        aaggggCGAGTAAAAGAGGGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGaacgt                                                                                                                                                                                                                           <  2:484239/248‑4 (MQ=255)
                        aaaGGGCGAGTAAAAGAGTGCGAAGAGATCATGATAAAACTGTTTGGCGAACAGGAGGCATACGATGAAGAGCAGACGCAGCAAGCCCGTGTTCGGGATCTGTTTAATCGCCGCCCTTCACCTGTTGTTCTGTTTGTAGCCGCAATCTGGCCATGTCAGGTGATCCCATAGTTCGCAATTTACACTTGTGGCCCGCGACTCGTCTGTATGTTGGGCTTGGGTGTTGGGAACAACGCggtgcgggggaacgg                                                                                                                                                                                                                           >  2:23149/1‑238 (MQ=255)
                             gctaGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTCCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGa                                                                                                                                                                                                                      <  1:224068/248‑1 (MQ=255)
                                            cGAAGAGATGATGATAAAACTGTTTGGCGAACCGTTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCACTGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCACAAACGCGGCGCTGGGGAACGTGATGATTAGCCTGTTCttt                                                                                                                                                                                                        >  1:196033/1‑250 (MQ=255)
                                            cGAAAAGATGATGAAAAAACCGTTTTGCGAACCGGTGGCTTTCGATGAAAAGCGGCAGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCCTTTTCCTTTTGTTCTGTTTGTTGACGCCATCTGGACCTGCCATGTGATCCCAATGTTCTCCATTTACACCTTTGGCCCGCAAATCGTTAGTTTGTTGGGATTGAGTGTTGGCAAAAACGCGGCGCTGGGGTACGTGGTGAGTAGCCTGTTGttt                                                                                                                                                                                                        >  2:632653/1‑250 (MQ=255)
                                             gAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAATCGCGGCGCTGGGGCACGTGGTGATTAGCCTGTTCTTTAt                                                                                                                                                                                                      >  2:341949/1‑251 (MQ=255)
                                                   atgatgatAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCAATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGACTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCgg                                                                                                                                                                                                <  1:4070/251‑1 (MQ=255)
                                                      atgatAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTg                                                                                                                                                                                             >  1:90872/1‑251 (MQ=255)
                                                       tgataaaccTGTTTGGCGAACCGGTGGCTATCGATCAAGAGCAGCCTCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTGTCCTTTTGATCCGCTTGTTGCCGACATCTGGACCTGCCAGGCGATCCTACTGCTCGCCATTTACACCCTTGGACCGAAAATCGTTGGTTTGTTGGGATTGGGCATTAGCAGAAAAGCGGAGTTGGGGTACGTGGGGTTTAGCCTGTTATTTATGCTCGGCTGt                                                                                                                                                                                            >  2:515448/9‑251 (MQ=255)
                                                            aaaCTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTcc                                                                                                                                                                                       >  2:118849/1‑251 (MQ=255)
                                                             aaCTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTccg                                                                                                                                                                                      <  1:329949/251‑1 (MQ=255)
                                                                 gTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTccgccg                                                                                                                                                                                   <  1:218096/250‑1 (MQ=255)
                                                                                ggtcttctatgaatagacgcCGCAGCAAACCCTTTTGCGCGATCTCTTTATTCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCAGTTACACCTTTGGCCCGCAAATCGTTGGGTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTTGGGACCGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAaca                                                                                                                                                                   <  2:622451/233‑1 (MQ=255)
                                                                                             gagCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCt                                                                                                                                                                                                               >  1:3976/1‑194 (MQ=255)
                                                                                                 agcCGCCGCAAACCCGTTTTCGCGATCTGTTTAATAGCCGCCCTTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCCTTTCCACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGACCGTGGTGATTAg                                                                                                                                                                                                                  >  2:402800/1‑187 (MQ=255)
                                                                                                 agcCGCAGCAAACCCGTTTTTGCGATCTGTTTAATCGCCGCCAGTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCCGCCAGGTGATCCCAATGTTCGCCATTTACCCCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAg                                                                                                                                                                                                                  <  1:402800/187‑1 (MQ=255)
                                                                                                     gcagcaAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCAttgt                                                                                                                                              >  1:237571/1‑251 (MQ=255)
                                                                                                           aaCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTg                                                                                                                                        >  1:243067/1‑251 (MQ=255)
                                                                                                           aaCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGCGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAATAACTGCCGGACGGCGTCCATTGTTGATTg                                                                                                                                        >  1:635782/1‑251 (MQ=255)
                                                                                                                  tttCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTTTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCAGTGTTGATTGGCAGCt                                                                                                                                  <  1:580026/250‑1 (MQ=255)
                                                                                                                            tGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGACATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGGCCATTGTTGATTGGCAGCTTTGCCatgatg                                                                                                                       >  2:584196/1‑251 (MQ=255)
                                                                                                                             gTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCat                                                                                                                           <  1:178797/246‑1 (MQ=255)
                                                                                                                                 aaTCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGCCCTGCCAGGTGAGCCCAATGTTCGCCATTTACACCTTGGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgct                                                                                                                  <  1:616268/251‑1 (MQ=255)
                                                                                                                                    cgccgcCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAacac                                                                                                                                                                  >  2:155807/1‑200 (MQ=255)
                                                                                                                                    cgccgcCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAacac                                                                                                                                                                  <  1:155807/200‑1 (MQ=255)
                                                                                                                                    cgccgcCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgctggc                                                                                                               >  1:377468/1‑251 (MQ=255)
                                                                                                                                        gccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAACTCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGCCGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgctggctctg                                                                                                           >  2:565165/1‑251 (MQ=255)
                                                                                                                                           aTTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGcgctggg                                                                                                                                                                                                                                 >  2:214403/1‑124 (MQ=255)
                                                                                                                                           aTTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGcgctggg                                                                                                                                                                                                                                 <  1:214403/130‑7 (MQ=255)
                                                                                                                                             tttCCTTTTGTTCTGTTTGTTGCCGCCATCTGGGCCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAAGCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCgg                                                                                                       <  1:142829/250‑1 (MQ=255)
                                                                                                                                              ttCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCgggt                                                                                                     >  1:460681/1‑249 (MQ=255)
                                                                                                                                                 cTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTTCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGTCGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTg                                                                                                  >  1:182377/1‑251 (MQ=255)
                                                                                                                                                 cTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTGTTg                                                                                                  >  1:244881/1‑251 (MQ=255)
                                                                                                                                                      gTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTTGGGCt                                                                                             >  1:209545/1‑251 (MQ=255)
                                                                                                                                                           gtttgttGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATcc                                                                                        >  1:436444/1‑251 (MQ=255)
                                                                                                                                                           gtttgtgGCCCCCATCTGGACCTGCCAGTGGATCCCAATGTTCGCCATTTACTCCTGTGGCCCGCAAATCGTTGTGTTGTTGGGATTGGTTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTGCTTTAGGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGTCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATcc                                                                                        <  2:237571/251‑1 (MQ=255)
                                                                                                                                                                     gccATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGTTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCCAGATGACGCTGGCGATGACGGTTTTtggggtggttcacgtgatggg                                                                              >  2:367092/1‑234 (MQ=255)
                                                                                                                                                                      ccATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGGTCCCGGGTATgggg                                                                             >  1:404272/1‑251 (MQ=255)
                                                                                                                                                                          cTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCACATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTACACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTgggggttgtgccgggtttggggatct                                                                         >  2:555232/1‑229 (MQ=255)
                                                                                                                                                                           tggcccTGCCAGGTGATCCCAATGTTTGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTGGGTATTGGTTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGAt                                                                                                                                                                                 <  2:22434/142‑1 (MQ=255)
                                                                                                                                                                           tGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGAt                                                                                                                                                                                 >  1:22434/1‑146 (MQ=255)
                                                                                                                                                                           gggACCTGCCAGGGGATCCCAAGGTTCGCCATTTACGCCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATctg                                                                        <  2:90185/250‑1 (MQ=255)
                                                                                                                                                                            ggACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATct                                                                         <  1:118849/249‑1 (MQ=255)
                                                                                                                                                                             gACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGGATTCCGCCGATGCTGTGGTTCAACACTGCCGGAAGGCGTCCATTGTTGATTGGCAGCTTTGCCATGGTGAATCCGGAGCTGGCGGTTTTGGGGCtgggccacgacattgggggttggt                                                                      >  2:22015/1‑228 (MQ=255)
                                                                                                                                                                               ccTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATctggctg                                                                    <  1:493978/251‑1 (MQ=255)
                                                                                                                                                                                      ggTGATCCGAATGTGCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATTCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGAGTGGCAGCTTTGCCATGATTACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGa                                                             <  1:646487/251‑1 (MQ=255)
                                                                                                                                                                                             ccAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCtt                                                      <  1:112850/251‑1 (MQ=255)
                                                                                                                                                                                                 tGTTCGCCATTTACCCCTTTGGCCCGCAAATCGTTGGTTTGGTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgctggc                                                                                                               <  2:557159/190‑1 (MQ=255)
                                                                                                                                                                                                 tGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgctggc                                                                                                               >  1:557159/1‑190 (MQ=255)
                                                                                                                                                                                                  gTTCTCCAGTTACCCCTTTGGCCCGCAAATCGGTGGTTTTGGGGGATTGGGGGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCCCGGCTGTATTCCGCCGATGCTGTGGTTAAACCCTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATTACGCGGGCGCTGGCTGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGTGCTTTGct                                                  <  2:218386/250‑2 (MQ=255)
                                                                                                                                                                                                        cATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTTCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATg                                          <  1:640553/251‑1 (MQ=255)
                                                                                                                                                                                                         aTTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑c                                        <  1:555232/251‑1 (MQ=255)
                                                                                                                                                                                                             acacCTTTGGCCCGCAAATCGTTGGGTTGTTTGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCtt                                     <  2:171072/250‑1 (MQ=255)
                                                                                                                                                                                                                       cccGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGTCGCTGGGGAACGTGGTGATGAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGTTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCTTTTTCTCTGGCg                          <  1:410519/251‑1 (MQ=255)
                                                                                                                                                                                                                           aaaaTCGTGGGGTTGATGGGATGGGGTGTTGGCAAAAACGCGGCGCTGGGAAACGTGGTGATTAGCCTGTTCTTTATGCGCGGCTGTATGCCGCCGATGCTGTGGTTAAACACTGCCTGACGTCGTCCATTGTTGATTGGCAGCTGTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCGTATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCTTTGTCTCTGGCGGGcc                      <  2:644146/250‑1 (MQ=255)
                                                                                                                                                                                                                                        tgttgGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCTTTTTCTCTGGCGGGCCGGGTAATTTGCa          <  1:565165/250‑1 (MQ=255)
                                                                                                                                                                                                                                         gttgGGATTGGGAGTTGGCAAAAACGCGGCGCTGGGGAACGCGGTGATTAGCCTGTTCTTTATGCTCGGATGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCTATGATGACGCTGGCGCTGGCGGCTTTGGGGCTGATACCGGGTATGGGGATCTCGCTGGGTGGGGTGGCCTTTGCGGAG‑TATGGCCTTTTT‑TCTCGCGGGACGGGGAATTTGCAgg        >  1:43415/1‑250 (MQ=255)
                                                                                                                                                                                                                                               aTTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGGATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGGTGATCCCGGGTTTGGGGGACTGGCTGGGAGGGATGGCGCTTGCGGGGGTTTG‑C‑TTTTTTTCTGGCGGGCCgggggattttgcatgggctc  >  2:19161/1‑234 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
GAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCT  >  NC_000913/46389‑46878

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