Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,243,825 C→T A235T (GCC→ACC)  ycgR ← flagellar brake protein YcgR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,243,8250CT98.2% 156.4 / NA 54A235T (GCC→ACC) ycgRflagellar brake protein YcgR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/29);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATG  >  NC_000913/1243578‑1244070
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
aGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTtc                                                                                                                                                                                                                                                    <  2:360065/251‑1 (MQ=255)
  aTCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTtctc                                                                                                                                                                                                                                                  >  1:158952/1‑251 (MQ=255)
               gCAATATGCGCTGGTGGTGTCATACGCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCgagagaga                                                                                                                                                                                                                                      <  1:126608/250‑1 (MQ=255)
                           ggtggtGTCATACGCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTg                                                                                                                                                                                                                         <  1:131114/251‑1 (MQ=255)
                               gtgtCATACGCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTTCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATATTCCGCTGTAAt                                                                                                                                                                                                                     <  2:51132/251‑1 (MQ=255)
                                   cATACGCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGcc                                                                                                                                                                                                                 >  1:473782/1‑251 (MQ=255)
                                        gCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGCAAAGGCAGTCAGCAAAATCAACGATTTGGATGCTGACGAGATCCTCGAACACGTAAGGGGAAATTACCCCTCGCCGCAGGCTCCCCGCTAAATCCACAAACTTGAGCCGGCCATGGACCGCAAGGAAATCAGTCGCCCTATTTGTCCGCTTTTTTCCCGGCTTCTCGCTCGGGGGGGTAAAAAATTCGCTGGAATTGCCGCTccc                                                                                                                                                                                                           >  2:500672/1‑250 (MQ=255)
                                          tAACGGGGCAGGTGCGCTGCTACGGAATTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGCTGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCCGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGGAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCAc                                                                                                                                                                                                          >  2:417896/1‑251 (MQ=255)
                                           aaCGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGCAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAACGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACACACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCAcc                                                                                                                                                                                                         >  2:258217/1‑251 (MQ=255)
                                           aaCGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCAcc                                                                                                                                                                                                         >  1:4219/1‑251 (MQ=255)
                                                 gCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCggg                                                                                                                                                                                                   <  1:531681/251‑1 (MQ=255)
                                                   aGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCGCCGTCGGGCt                                                                                                                                                                                                 <  2:57801/251‑1 (MQ=255)
                                                   aGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCt                                                                                                                                                                                                 >  2:516862/1‑251 (MQ=255)
                                                            gctACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCTATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTaa                                                                                                                                                                                         <  2:27065/250‑1 (MQ=255)
                                                             ctACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTa                                                                                                                                                                                          >  2:392873/1‑248 (MQ=255)
                                                                           gTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGCTGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGAGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTACTTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTc                                                                                                                                                                           >  2:103906/1‑249 (MQ=255)
                                                                                    gAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAgt                                                                                                                                                                >  1:131610/1‑251 (MQ=255)
                                                                                     aaaaaGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAgtg                                                                                                                                                               >  2:375094/1‑251 (MQ=255)
                                                                                           gCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATg                                                                                                                                                         <  2:310627/251‑1 (MQ=255)
                                                                                            cGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATgg                                                                                                                                                        >  2:150360/1‑251 (MQ=255)
                                                                                                cAGCAAAATCAACGAGTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTGGAGCAGGCACTTGACGCGATGTAAATCAGTCGCGCTCTTTGTCCGCTTTTTCCCGGGTGTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACTTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTCATGGTTTc                                                                                                                                                    <  2:283368/251‑1 (MQ=255)
                                                                                                   cAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCAtt                                                                                                                                                 >  1:550525/1‑251 (MQ=255)
                                                                                                    aaaaTCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGGGGGGATGGATGCATTc                                                                                                                                                >  2:358746/1‑251 (MQ=255)
                                                                                                          aaCGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACAGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTACCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGATCCACCGTCGGGCTGACGTTAAGAAAACGGAAGATCAGACGGGGAGGGGAGATGGATTCATTCTTGccc                                                                                                                                          >  1:427516/1‑250 (MQ=255)
                                                                                                                ttGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCatcaat                                                                                                                                     <  1:245919/250‑1 (MQ=255)
                                                                                                                 tGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCTCCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCatcaatc                                                                                                                                    <  1:258217/250‑1 (MQ=255)
                                                                                                                      gcagACGAGTTCCTCGAACACGGAGCGCGAAATCACCCTGCGCCGCAGGCACCGCGCTATATCGTCAAACTTTAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATGGCCGCTCGACCTTCGGGCGGACGTTAAGAATACGTAAGCTCAGACGGGGAGTGGTGATGGTGGCATTCTTTCCATCAATCACTTTg                                                                                                                              <  2:67409/248‑1 (MQ=255)
                                                                                                                      gCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACttt                                                                                                                               <  1:177943/250‑1 (MQ=255)
                                                                                                                                    cGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATgg                                                                                                                 <  1:620167/250‑1 (MQ=255)
                                                                                                                                    cGAACACGCAACGCGAAATCACCCTGAGCAGCAGGCACCGAGCTATATCCACAACCTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGAGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGCGAGCGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGGGTGGTGATGGTTTCATTCTTGCCCTCCATCACTTTGCGCTCGCTGATga                                                                                                                 >  2:598202/1‑249 (MQ=255)
                                                                                                                                     gAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATgg                                                                                                                 >  1:144267/1‑249 (MQ=255)
                                                                                                                                      aaCACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATgg                                                                                                                 >  1:121987/1‑248 (MQ=255)
                                                                                                                                             aGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTg                                                                                                       <  1:84031/251‑1 (MQ=255)
                                                                                                                                                            ggcgcCGCAGGCGCCGCGCTATATCTGCAAACTTGAGCAGGCACTGGACGCGATGTAAAGCAGGCGCGCACGTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGGTGTAATTGCCGCTCCACCGTCGTGCTGACGTTAAGAAAACGGAAGCTGAGACGGGGAGTGGTGATGGTTTCATTCTTGGCATCGATCACTTTGCGCGCGCTGATGGAGATTAACTGGGCGTCAAAGTGaaa                                                                                        <  1:323773/250‑1 (MQ=255)
                                                                                                                                                                ccgcAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAAGCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGGTTCATTCGTGCCATCAATGACTTTGCGCTCGCTGATGGAGATTAACGGGGCGGCAAAGGGaaaatca                                                                                    <  2:385941/251‑4 (MQ=255)
                                                                                                                                                                             gcTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCa                                                                       <  1:150360/251‑1 (MQ=255)
                                                                                                                                                                                atatCTACAAACTTGAGCAGGCACTGGACGCGATGGAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCAt                                                                      <  1:417896/249‑1 (MQ=255)
                                                                                                                                                                                 tatCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGtt                                                                   >  2:461413/1‑251 (MQ=255)
                                                                                                                                                                                  atCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTg                                                                  >  1:344057/1‑251 (MQ=255)
                                                                                                                                                                                    gtACAAACTTGAGCAGGCACGGGACGGGATGTAAAGCAGTCGCGCACTTTGTCAGCTTTTTCCCGGGTTTCTCTCTCGAGAGAGAAAATAATCCGCTGTAAGGGCCTCTGCACCGTCTGGCTAACGTTAAGAAAACGGAAGCTCAGACGATGCGTGGTGATGGTTTCATTCTTGCCATCAATCGCTTTGCTCTCGCTGATGGAGAGTAAGGGGGCGTCAAAGTAAAAAATTCCGCTTTGCCCTATGTTGAc                                                                <  2:374705/250‑1 (MQ=255)
                                                                                                                                                                                      aCAAACTTGAGCTGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGa                                                                 <  2:339442/248‑1 (MQ=255)
                                                                                                                                                                                      aCAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACt                                                               <  1:521334/250‑1 (MQ=255)
                                                                                                                                                                                       cAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTc                                                             >  2:475180/1‑251 (MQ=255)
                                                                                                                                                                                             tGAGCTGGCACTGGACGCGATGTAACTCAGCAGCGAACTATGTCCGCTTTTTCACGGGTTTCTCGATAGAGAGAGAAAATAATCCGATGTAATTGCCGCTCCCCCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGGGGTGATGGTTTCATTCTTGCCATAAATAACTTTGCGCTGGCTGATGGGGATTAATTGTGCGTAAAAGTGACAAACACCCAATTGCCCCCTGTTGACTTCAAtatg                                                       >  2:356656/1‑248 (MQ=255)
                                                                                                                                                                                               aGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTg                                                       >  2:278399/1‑249 (MQ=255)
                                                                                                                                                                                                         ggACGCGATGGAAAGCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCGCCACCGTCGTGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTCGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCa                                            <  2:246678/250‑1 (MQ=255)
                                                                                                                                                                                                               gATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACGTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCcttc                                     <  2:223409/251‑1 (MQ=255)
                                                                                                                                                                                                                     aaaTCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGt                                 <  1:540758/250‑1 (MQ=255)
                                                                                                                                                                                                                     aaaTCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGt                                 >  1:229542/1‑250 (MQ=255)
                                                                                                                                                                                                                             cgcgCACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCa                        <  1:597641/251‑1 (MQ=255)
                                                                                                                                                                                                                              gcgcACTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCTCCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAg                       <  2:467985/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ccTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCtt                   <  1:23588/250‑1 (MQ=255)
                                                                                                                                                                                                                                  aCTTTGTCCGCTTTTTCCCGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCtt                   <  1:530641/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccGGGTTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATg  >  1:419515/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
AGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATG  >  NC_000913/1243578‑1244070

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