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,3310TC97.8% 158.2 / ‑4.2 46intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (20/25);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

TGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTA‑TGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652088‑1652573
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
aatcggcgctcttccgatagtGGAGTTCGCCTTGGGGAGCTGTGGGCATGTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTG‑CACCTCagag                                                                                                                                                                                                                                                <  2:483910/230‑1 (MQ=255)
 ggTGTGTTTCTGGTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGAACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAGCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCagagag                                                                                                                                                                                                                                              <  1:95988/251‑1 (MQ=255)
  gtgtgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGt                                                                                                                                                                                                                                             <  1:589860/251‑1 (MQ=255)
     tgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCCATCGACGCTCTCCGGAGTCAAATGACATTAACTAGACTGAGAAAAAATCATATCATTGATGAT‑CACCTCAGAGAGGATg                                                                                                                                                                                                                                          >  1:144018/1‑251 (MQ=255)
     tgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATAGGGGCATTTATCTTTTTAAGCCAGCAATCGACGCTCTCCGGAGACACATGACATTAACGGGACTGAACAAAGAGCATATCATTGATGTTTCACT‑CAGAGTGTATg                                                                                                                                                                                                                                          >  1:834063/1‑251 (MQ=255)
      gtgtCTGTTTATACTGGAGTTCGCCCATGTGAGCTGTGTGCACTTGGATGTGAAGACATAGAACTGAAAAATGGAACAATGATTATCATGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCGCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGGT‑CACCTCAGAGAGTATgg                                                                                                                                                                                                                                         <  2:549069/248‑1 (MQ=255)
      gtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGAAGCTCTCCGGAGTCAGATGACATTAACGAGACTGGGCAAAGAGCATATCATTGTTGTT‑CACCTCCGAGAGTaagg                                                                                                                                                                                                                                         >  1:642813/1‑248 (MQ=255)
      gtTTCTGGGTATACTGGAGGTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGAGCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCATAGAGTATgg                                                                                                                                                                                                                                         <  2:851047/251‑1 (MQ=255)
            gTTTATACGGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGcagaa                                                                                                                                                                                                                                    <  2:385144/250‑1 (MQ=255)
              ttATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGCAACAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGTAGACAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑TACCTCCGAGAGGATGGcagaac                                                                                                                                                                                                                                   >  1:716071/1‑249 (MQ=255)
                ataCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGTGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGcagaacaga                                                                                                                                                                                                                                <  2:485378/250‑1 (MQ=255)
                  aCTGGAGTTCGCCCTGGTGAGCTTTGGGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCATGAGAAAGTTAGCAAAAGACCGTTTCACGGTACCAAAACCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGaaaa                                                                                                                                                                                                                             <  2:595286/251‑1 (MQ=255)
                  aCTGGAGGTCGCCCTGGTGAGCTGTGTGCACGTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAAGTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGAATGAGCAAAGAGCATATCATCGATGTT‑CAACTCAGAGAGTATGGCAGAACAGaaaa                                                                                                                                                                                                                             <  2:97831/251‑1 (MQ=255)
                     ggAGTTCGCCCTGGTTAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAAc                                                                                                                                                                                                                           <  2:523949/250‑1 (MQ=255)
                                  gTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATTGCAGAACAGAAAAACAAAAATGCACCttt                                                                                                                                                                                                             <  1:700547/251‑1 (MQ=255)
                                  gTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGTCATTCATCTTATTAAGCCCGCAAACGACGCTCTCCGGGGTCAAATTACATTAACGAGACTGCGCAAAAAGCATATCATTGATTTT‑CAACTCCGAGGGAATGGCAGAACAGAAAAAAAAAAAAGAACCttt                                                                                                                                                                                                             >  2:700547/1‑251 (MQ=255)
                                    gAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTg                                                                                                                                                                                                            <  2:495857/250‑1 (MQ=255)
                                           gtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCaa                                                                                                                                                                                                    <  1:634808/251‑1 (MQ=255)
                                                        ggaggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACAGTTTCACGGTACCAAAAACACAGGCGGGAAACAAAAGGGTCATTTATCTTATTAAACCCGCAAAAAGCGCTCTCCGGGGTCAGATGACATTAACGAAACTGAGCAAAAAACATATTATTGATGGT‑AACCTCCGAGGGTATGGCAAAACAACAAAAAAAAAATACAACTtttttttttcaaatgaaatgtaa                                                                                                                                                                                        >  2:646672/1‑235 (MQ=255)
                                                                     cTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGGGGGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGGaaaaa                                                                                                                                                                           >  1:648486/1‑250 (MQ=255)
                                                                       gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTGTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                                                                                                                                                                             <  2:681456/246‑1 (MQ=255)
                                                                                        tgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCCTTCATCTTATTAAGCAAGCAATCGTCACTCTCCGGAGTCAGATTACCTTAACGAGACTTAGCAAAAAGCATATCATTGATGTT‑CACCTCACGGGGGATGGCACAACAGAAAAACAAAAAAGAACCTTTTTTTTTTAACCTGAAGGGTACGCGAGAGTAAAAAAATTTTGTTGACC‑TTTTAcc                                                                                                                                                      >  2:112645/1‑251 (MQ=255)
                                                                                                         gCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCCTCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTAACGTTGACTCAATAAGGc                                                                                                                                      >  1:866587/1‑250 (MQ=255)
                                                                                                           aaaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCCTTCATCTTATTAAACCAGCAATCGACGCTCTCCCGAGTCAGATGACATTTACAAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCCAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTTACTCAATAAGGc                                                                                                                                      >  1:225567/1‑248 (MQ=255)
                                                                                                                     ttCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACg                                                                                                                          <  1:20545/250‑1 (MQ=255)
                                                                                                                                aaaaaCACAGGCGGGAACCAATCGGGTCATTCATCTTTTTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTTACGAGACTGGGCAAAAAGCATATCAATGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAACTGAAGTGTCAGCGAGAGGAAAAAA‑TTATGGTGACCATTTTTCCGTTGAATCAATAAAGCCGATGGTGGACGCcgcgcaaatac                                                                                                              >  1:15780/1‑244 (MQ=255)
                                                                                                                                 aaaaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGATCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACg                                                                                                             <  1:641987/251‑1 (MQ=255)
                                                                                                                                          gCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTTTTTCTCAACCTTAAGTTTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTGACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACt                                                                                                    <  1:566801/251‑1 (MQ=255)
                                                                                                                                          gCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACt                                                                                                    <  1:112645/251‑1 (MQ=255)
                                                                                                                                           cGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTc                                                                                                   >  1:824473/1‑251 (MQ=255)
                                                                                                                                                           cattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCat                                                                                   >  2:176355/1‑251 (MQ=255)
                                                                                                                                                                             aGCAAACGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACtt                                                                 <  2:225567/251‑1 (MQ=255)
                                                                                                                                                                                  tCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTA‑TGcc                                                           >  1:848951/1‑251 (MQ=255)
                                                                                                                                                                                    tacactctCCGGAGACAGATGCCATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAACTGCACCTTTGTTTTTCAACCTGATGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACAGTAGACTCAATAAGGCAGATGTGTGACGCAGCGATAAAACGTGACGGACTCCCCCATCGAACATCATATCAGTCCAGAAATACTTA‑TGCctg                                                         >  2:111579/5‑251 (MQ=255)
                                                                                                                                                                                    gACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTA‑TGCctg                                                         <  1:142366/251‑1 (MQ=255)
                                                                                                                                                                                       gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGACGAGACATACTTA‑TGCctgctg                                                      >  2:108360/1‑251 (MQ=255)
                                                                                                                                                                                        ctctCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGCACAGAAAAACACAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTA‑TGCCTGCTgg                                                     >  1:51870/1‑251 (MQ=255)
                                                                                                                                                                                        ctctCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTA‑TGCCTGCTgg                                                     >  1:533199/1‑251 (MQ=255)
                                                                                                                                                                                            ccGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaaaa                                                                                                                                                                          <  2:602000/132‑1 (MQ=255)
                                                                                                                                                                                            ccGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaaaa                                                                                                                                                                          >  1:602000/1‑132 (MQ=255)
                                                                                                                                                                                                      aTGACATTAACGAGACTGCGCAAAGAACATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCGGAAGAGTCAGCGAGAGTAAAAAA‑TTATGGTGACAATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCAGACTCCGCCATCGAAAATCATATAAGTCGAGACAAATTTA‑TGCTTGCTGATCGCTGACAACTgg                                       >  2:598931/1‑251 (MQ=255)
                                                                                                                                                                                                        gACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTA‑TGCCTGCTGGTCGCTGACAGCTGGt                                      <  1:241575/250‑1 (MQ=255)
                                                                                                                                                                                                                       gAGCAAAGAGCATATCATTGATGTT‑CACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCGTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGAAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTA‑TGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAtt                        <  1:367966/249‑1 (MQ=255)
                                                                                                                                                                                                                               agCATATCATTGATGTT‑CACCTAAGAGAGTATGGCAGAACAGCAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGCGAGTAAAAAA‑TTATGGAGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAAAGAGCCGGACTTCGCCATCGAAAATCATATCAGTCGTGACAT‑CTTAATGCCAGCTGGTTGCTGACAGCTTGTGCTAAACAGGCATTTATAgacaa               >  2:251736/1‑246 (MQ=255)
                                                                                                                                                                                                                                      cATTGATGTT‑CACCTCAGAGAGTATGGCAGAAAAGAAAAACAAAAATGCACCTTGGTTTTTCAACCTGAATTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAAGAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTA‑TGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATgg       <  2:848951/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aTGTT‑CACCTCAGTGATTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGGGACCATTTTACCGTTGACTCAATAATGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCTGCCATCGAAAATCATATCAGTCGAGACATACTTA‑TGCCTGCGGGGCGCGGACAGCTGGTGCTAGCCCGGCAGTTATAGCAACCCAGATGGGCCAt  <  2:99998/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
TGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTT‑CACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAA‑TTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTA‑TGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652088‑1652573

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