Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,652,331 T→C intergenic (+2346/+397) ydfD → / ← ynfP lysis protein/protein YnfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,652,3310TC95.7% 156.5 / NA 46intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (20/24);  minor base A (2/0);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGG  >  NC_000913/1652086‑1652568
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
tATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCag                                                                                                                                                                                                                                               <  1:105548/250‑1 (MQ=255)
 aTGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagag                                                                                                                                                                                                                                             <  1:553493/251‑1 (MQ=255)
   ggTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagagag                                                                                                                                                                                                                                           <  1:36363/251‑1 (MQ=255)
      gtgtgtCTGTGTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATATGAAAAATGGAACAATGATGATCAGGAGAAATGTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGACAGCAATCGACGCTCTCCGGATTCAGATGACATTAACGAGACTGATCAAAGAGCATATCATTGATGTTCACCTCTGAGAGTa                                                                                                                                                                                                                                         <  2:497580/250‑1 (MQ=255)
              gTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaac                                                                                                                                                                                                                               <  1:580776/251‑1 (MQ=255)
              gTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcag                                                                                                                                                                                                                                   <  2:174841/248‑1 (MQ=255)
              gTTTAGACTGGAGTTCGCCCTGGAGAGTTGTGTGCACTTGGAGGGGAGGACATAGATCTGAAAAATGGATCAATGATGATCAGGAGAAATTTAGCAAAAGACCGTGTCACGGTACTAAAAACACAGGCAGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATAGCATTGATGTTCACCTCAGAGAGTATGGcagaa                                                                                                                                                                                                                                 <  1:365565/250‑1 (MQ=255)
                    aCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaa                                                                                                                                                                                                                          <  1:278245/250‑1 (MQ=255)
                       ggAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCCACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTTATGTTCACCTCAGAGGGTATGGCAGAACAGGAAAAc                                                                                                                                                                                                                      >  1:263431/1‑250 (MQ=255)
                                   ggTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCtt                                                                                                                                                                                                         <  2:128894/251‑1 (MQ=255)
                                           gtgtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTc                                                                                                                                                                                                <  1:428921/251‑1 (MQ=255)
                                                    tGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCAACGGGGCAATCTTCTTATTAAGACAGAAATCGACGCCCTCCGGAGTCAGATGACATTTACGAGACTGAGCAAAGAGCATATCATTGCTGTTTACCTCAAAGAGTATGGCAGAACAGCAAAACAAAAATATACATTTGTTTTTCAACCttaa                                                                                                                                                                                        >  2:201900/1‑247 (MQ=255)
                                                     ggATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTTTTAAACCAGCAATCGGCGCTCTCCGGAGTCAGATGACATTAACGAAAATGAGCACAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGCTGAAAAGAAAAAAAAAAAAGGCCCTTTGTTTTTCAACCTGAAgt                                                                                                                                                                                      >  2:115058/1‑251 (MQ=255)
                                                        tGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGGCATTCATTTTATTAAGCCAGCAATCGACGCATTCCGGAGTAAGAAGAAATTAACGAGACTGAGGAAAAAGCATAATATTGGTGTTCAACTCAGAGAGTATGGCAGAAGAGGAAAAAAAAAAAGGACCTTTTTTTTTTAAACTGGAGTGTc                                                                                                                                                                                   >  1:181145/1‑251 (MQ=255)
                                                           taggaCATAGATCTGAAAACTGGAACAATGCTGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAAAGGGTCATTCATCTTATTAAACCAACAAACGACGCTCTACGGGGTCAGATGAAATTTACGGGACTGAGCAAAAAGCATATTATTTTTTTTCACCTCAGTGAGGATGGCAAACAAAAAAAACAACAATGCACCTTTGTTTTTTAACCTGaaaagtca                                                                                                                                                                                  >  2:332264/2‑243 (MQ=255)
                                                                      tCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTTATGTTCACCTCAGAGAGGATGGCAGAAAAGAAAAAAAAAAATGCACCTTTGTTTTTCAACCTGAAGAGTCCGCGAGAGGaaaa                                                                                                                                                                     >  1:305916/1‑250 (MQ=255)
                                                                                 gAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaa                                                                                                                                                                                                                          <  1:469893/189‑1 (MQ=255)
                                                                                 gAACAATGATGAACAGGAGAAATTTAGCAAANGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACAGTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc                                                                                                                                                         <  2:277681/251‑1 (MQ=255)
                                                                                 gAACAAAGATGAACAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCCCCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGCCAGCGAGAGTAAACAATTATGGTGAc                                                                                                                                                         <  2:314615/251‑1 (MQ=255)
                                                                                          tgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGCAACCAATCGGGGCAATCATCTTATTAAGGCAACAATCGACGCAATCCGGAGTCAGATGACATTAACGAGACTGAGAAAAAAGCATATAATTTATGTTTACATCAGGGAGTATGGCAGAAAAGCAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGCGAGGAAAAAATTATGGTGACCAttttcac                                                                                                                                                >  1:242776/1‑248 (MQ=255)
                                                                                           gatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACATAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCATCCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTAcc                                                                                                                                                <  2:32246/250‑1 (MQ=255)
                                                                                                                 gACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAACCGGGCCAATCTATTTATTAAAACAGAAATAGGCGCCCTCCGGAGTCAGATGAACTTAACGAGAACTAGCAAAGCACATATCATTGAAGTTCAACTCCGAGAGGATGGCCGAAACAAAAAACAAAAATGCACCTTTTGTTTTTCACCTGAACTGTAAGCGAGGAGAAAAAATTATGGTGCACATTTTGCCGGTGCATCAATAAGGCAGAtgtg                                                                                                                         >  2:392646/1‑251 (MQ=255)
                                                                                                                    cGTTTCACGGTACCAAAACCACAGGCGGGAACCAACCGGGTCCTTCATCTTATTAAGCCAGCAATCCACGCCCTCCGGAGTCAGATGACATTAAAGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGGGAGTATGGCAGAACCGAAAAACAAAAAAGCACCTTTTTTTTTCAAACAGAAGTGTCAGCGAGAGAAAAAAATTATGGGGAACATTTTATCCTTGACTCAATAAAGCAGATGTgggg                                                                                                                      >  2:452756/1‑250 (MQ=255)
                                                                                                                    cGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGTCCTTTATCTTATTAAACCCACAAACGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAAAAAAAAAGCACCCTTGTTTTTCAAACTTAAATGTCAGGGAGAGAAAAAAATTATGGTGACCCTTTTACCGTTGAATCCATAAGGCAGATGggggg                                                                                                                      >  2:152802/1‑250 (MQ=255)
                                                                                                                    cGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGACATTCATCTTCTTAAACCAGCAAACGAAGCACTCCGGGGTCAGATGACCTTAACAGGACCGAGAAAACAGCAAATCATTGATGTTCACCTCAGAAAGTATGGCAGAACAGAAAAACAAAAATGCCCCTTTGTTTTTCAACCTGAAATGTCCGCGAGAGTAAAAAATTTTTGGGAACCTTTTAACGTTTACTCCATAAGGCAGATGggggg                                                                                                                      >  1:582099/1‑250 (MQ=255)
                                                                                                                     cgtgcACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGAAGCTCTCCTGAGTCAGATGACATTAACCAGACTGAGCAAAGAGGATATCAATTATGGTCACCTCAGAGAGGATGGCAGAACAGAAAAACAAAAATGCCACTTTGTTTTTCAACCTGAAGTGTCAGCGAAAGTAAAAAATTATTGTGACCCTTTTACCGTTGACTCAATAAGGCAGATGTGGGAc                                                                                                                     <  1:339920/247‑1 (MQ=255)
                                                                                                                                ccAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaaa                                                                                                          <  1:332264/251‑1 (MQ=255)
                                                                                                                                  aaaaaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGGGGGACGCAGCGATaaaa                                                                                                         >  1:208052/1‑250 (MQ=255)
                                                                                                                                    aatcacaGGCGGGATCCAATTGTGGCATTCATCTTAATAAGCCAGCATTCGGCGCTCTCCGGAGTCAGATGACACGAGCGAGACTGAGCAAAGAGCATATCATTGATTTTCCCCTCAGAGAGTATGGCATAACAGAAAAACAAAAATGCACCTGTGTTTTTCAACCTGAAGTGTCAGCGAGAGGAAAAAATTATGCTGACCATGTTACCGTTGACTCGATAGGGCAGATGTGGGACGCAGCGCGAAAACgg                                                                                                      <  2:42473/248‑2 (MQ=255)
                                                                                                                                          aGGCGGGAACCCATCGGGTCATTCATCTTATTAAGCCAGAAATCGACGCTCTCCGGAGGCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAACAATGCACCTTTGTTTTTCAACCTGAAGCGTCAGCGAGAGTAAAACATTATGGTGACCCTTTTACCGTTGACTCAATAAGGCAGCTGTGGGACGCAGCGATAAAACGTGCCggc                                                                                                >  2:533626/1‑250 (MQ=255)
                                                                                                                                              gggAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTcc                                                                                            >  2:152822/1‑251 (MQ=255)
                                                                                                                                               ggAACCAATCGGGAAGTTCATCTTATTAACACAGCAATCGTCGATATCCGGAGTCAGATGACATTAAAGAGGCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAAAAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAACGAGAGTAAAAAATTATGGTGACAATTTTACCGTTGACTCAATAAGGCAGAGGTGGGACGCAGCGCTCAAAAGAGCCGGAcaccc                                                                                           >  2:376712/1‑247 (MQ=255)
                                                                                                                                                    cAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATc                                                                                      <  2:224884/251‑1 (MQ=255)
                                                                                                                                                    cAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTACGCCATc                                                                                      >  2:73935/1‑251 (MQ=255)
                                                                                                                                                           gTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATc                                                                               <  1:152802/251‑1 (MQ=255)
                                                                                                                                                                       tattaAGCCAGCAATCGACGCTGTCGGGAGTCAGATGACATTAACGAGACTGTGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGAATGGCAGAACAGACAAACAAACATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCCGATGTGGTACGCAGAGATAAACCGTGCCGGACTCCGCCATCGAACATAATCTCAGt                                                                       >  2:174257/1‑248 (MQ=255)
                                                                                                                                                                            aGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACAt                                                               <  1:15517/251‑1 (MQ=255)
                                                                                                                                                                                         cGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCCTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGGCCATTTTACCGTTGACTCAATAAGTCAGATGTGTGACGCAGCGATAAAACGTGCCGGACTCCGCCAGCGAAAATCATATCAGTCGAGACATACTTATGCctgct                                                  <  2:415910/251‑1 (MQ=255)
                                                                                                                                                                                             ctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTc                                              >  1:13617/1‑251 (MQ=255)
                                                                                                                                                                                                             cATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACATCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGCTGTGGGACGCAGCGATAAAACGTGCCGGCCTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATTCCTGCTGGTCCCTGACCGCTGGTGCt                              >  2:86625/1‑251 (MQ=255)
                                                                                                                                                                                                                   cGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCg                        <  1:389025/251‑1 (MQ=255)
                                                                                                                                                                                                                     agaCTGAGCAAAGAGCTTAGCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAGGCACCTTTGGTGTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACAATTTGACCGTTGTCTCAATAAAGCAGATGTGGGACGCAGCGATAAAACGTGCCGGTCTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGGCGCTGTCAGCTGGTGCTAAGCCGGc                      <  2:444509/251‑1 (MQ=255)
                                                                                                                                                                                                                     agaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGc                      >  1:265186/1‑251 (MQ=255)
                                                                                                                                                                                                                       aCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCCCCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAt                    <  2:98179/251‑1 (MQ=255)
                                                                                                                                                                                                                                         tCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCCTCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATg   >  1:522348/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATgg  >  1:510274/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
TATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGG  >  NC_000913/1652086‑1652568

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