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,3310TC98.4% 230.3 / NA 63intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (34/28);  minor base A (1/0);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652090‑1652575
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
gtgtgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCCCTCCGGAGGCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCagagagg                                                                                                                                                                                                                                             >  1:12301/1‑250 (MQ=255)
   tgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAATGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGAATg                                                                                                                                                                                                                                          >  1:299613/1‑251 (MQ=255)
        cTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AAACAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATAATTGATGTTCACCCCAGAGAGTATGGc                                                                                                                                                                                                                                        >  2:731927/1‑248 (MQ=255)
          gTTTATACTGGAGTTCTCCCTGGTGAGCTGTGTGCACTTGGATGTGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcaga                                                                                                                                                                                                                                     <  1:686123/249‑1 (MQ=255)
            ttATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AAACAATCGGGGCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCCGAGAGTATGGcagaac                                                                                                                                                                                                                                   >  1:260780/1‑249 (MQ=255)
             tataCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATAGGGGCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaac                                                                                                                                                                                                                                   >  1:499923/1‑248 (MQ=255)
               taCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGAACCAA‑CGGGGCATTTATCTTATTAAGCCAGCAATCC‑ACGCCCTCCCGAGTCAGATGACATTAAAGAGACTGAGCAAAGAGCATAACATTTATTGTCACCTCCGAGAGTAAGGccgaacagaa                                                                                                                                                                                                                               >  2:835437/1‑250 (MQ=255)
                aCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaa                                                                                                                                                                                                                             <  1:731927/251‑1 (MQ=255)
                            ccTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCAc                                                                                                                                                                                                                 <  2:84221/251‑1 (MQ=255)
                                gTGAGCTTTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGAGCAGAAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGTCGGG‑AACCATTCGGGTCATTCATCTTATTAAGCCAGCAAGTG‑TCGCTCTCCTGAGTCAGATGACATTAACGAGACTTATCAAAGGGCATATCATTGATTTTCACCTCAGATATTATGGCGGAACAGAAAAACAAAAATGCAccggt                                                                                                                                                                                                             <  2:668532/251‑4 (MQ=255)
                                  gAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTg                                                                                                                                                                                                            <  2:549951/250‑1 (MQ=255)
                                  gAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTg                                                                                                                                                                                                            <  1:785557/250‑1 (MQ=255)
                                   aGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaaa                                                                                                                                                                                                                       <  2:486833/238‑1 (MQ=255)
                                   aGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGGATGGCAGAACAGAAAAACaaaa                                                                                                                                                                                                                       >  1:486833/1‑238 (MQ=255)
                                         gtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCaa                                                                                                                                                                                                    <  1:14949/251‑1 (MQ=255)
                                                  gATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgtg                                                                                                                                                                                           <  1:473716/251‑1 (MQ=255)
                                                     gggAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCa                                                                                                                                                                                        <  2:27978/251‑1 (MQ=255)
                                                            cATAGCTCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCCATCATTTTATTAAGACAGCAAACGCAC‑CTCTCCCGAGTCAGGAGACCTTAACGAGACTGAGCAAAGAGCATATCATTTGTGTTCACCCCAGAGAGGAAGGGAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCATCCTTCAGTGTCAGCgag                                                                                                                                                                                   >  2:284650/1‑249 (MQ=255)
                                                             aTAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTTATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGagagg                                                                                                                                                                                >  1:654150/1‑250 (MQ=255)
                                                                  tCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAAACGGGTCATTTATCTTTTTAAGCCAGCAATCG‑ACGCTCTCCGGAATCAAATGACATTAACGAGACTGAGAAAAGAGCATATCATTGATTTTCAACTCCGAGGATATGGCAGAACAGAAAAAAAAAAATGTACCTTTTTTTTTCAAACTTAAATGTCAACGCGAGTacaa                                                                                                                                                                            >  2:502859/1‑247 (MQ=255)
                                                                       aaaaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTAt                                                                                                                                                                      >  1:92466/1‑251 (MQ=255)
                                                                            ggAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTg                                                                                                                                                                  <  1:416130/250‑1 (MQ=255)
                                                                             gAACAATGTTGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCCGGTCATTCCTCTTATTAAGCCAGCAATCG‑AAGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc                                                                                                                                                                >  1:86055/1‑251 (MQ=255)
                                                                             gAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAAATGACATTAACGAGACTGAGCAAAGAGCAAATCATTGATGTTCACCTCAGAGAGTATGGCAAAACAGAAAAAAAAAAAAGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc                                                                                                                                                                >  1:576516/1‑251 (MQ=255)
                                                                              aaCAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATCACATTAAAGAGACTGAGCAAAAAGCATATCCTTGATGTTTAACTCCGAGAGTATGGCAGAACAGAAAAACCAAAATGCACCTTTGGTTTTCAACCTGCAATGCCAGCGAGAGTAAAAAATTATgggtgac                                                                                                                                                               >  2:30979/1‑246 (MQ=255)
                                                                               aCAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCCGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTTAACCTGAAGGGTCAGCGAGAGTAAAAAATTATGGTGAccc                                                                                                                                                              >  1:509853/1‑250 (MQ=255)
                                                                                cAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAACCAGCAATCG‑CCGCCCTCCGGAGGCAGATTGCCTTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCCGAGAGGATGGCAGAACAGAAAAAAAAAAATGCAACTTTGTTTTTCAAACTGAAGAGTCAGCGAGAGGAAAAAATTATGGTGAACAt                                                                                                                                                             >  2:760226/1‑251 (MQ=255)
                                                                                    gatgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAAACGGGTCATTCATCTTATTAAGCCAGCAATCG‑GCCCTCTCCGGAGGCAGGTGAAATTAAAGAGACTGAGCCAAAAGCATATCATTGATGTTTACCTCCGAGAGGATGGCCGAAAAGAAACAAAAAAAAGCCACTTTGTTTTTCAAACTTAAGTGTCATCGTAAGTAAAAAAATATTGTGACCAtttt                                                                                                                                                          >  2:444258/1‑250 (MQ=255)
                                                                                           aGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAACACACAGGCGGG‑TACCAAACGGGTCATTCATCTTCTTAAGGCAGCAAACG‑ACGCCCTCCGGATTCAGATGACATTATCGAGAATCAGCCAAAAGCAAATCCTTGATGTTCAACACAAAGAGGCTTGCAGAAACGAAAAACAAAAATGCACCTTTGTTTTTTCACACGAAGTGTCACAGAGAGTAAAAAAATTTGGTGACCCttttttcgttct                                                                                                                                                  >  2:380285/1‑243 (MQ=255)
                                                                                           aGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTg                                                                                                                                                   >  1:705344/1‑250 (MQ=255)
                                                                                             gagaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AAACAATCGGGTTATTCCTCTTTTTAAGCAAGAAAACG‑ACCCTCTCCGGAGGCAGAAGACATTAACGGGACTGAACAAAAAGCATATCATTGATGGTCACCTCCGAGAGGATGGCAGAACAGAAAAACAAAAAAGCCCCTTTGTTTTTTAACCTGAAGGGTTAGCGAGGGGAAAAAATTATGGTGAACATTTTACCGTTGACt                                                                                                                                                >  1:402501/1‑251 (MQ=255)
                                                                                              agaAATTTAGAAAAAGACCGTTTCACGGTACCAAAAACACAGGAGGG‑AAACAAATGGGGCATTTAAATTATTAAAGAAGCAATAG‑AAGGTATCCGGAGTGAGATGATATTAAAGAGAATGAGAAAAGAGCATATCATTGATGTTTACATTAGAGAGTATGGCAGAAAAGAAAAAAAAAAATTAAACTTTGTTTTTTAAACTGAAGAGGTAGCGAGAGAAAAAAATTATGGTGAACATTTTATAGTTGACtt                                                                                                                                               >  2:128680/1‑250 (MQ=255)
                                                                                                 aaTTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCTCCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAAt                                                                                                                                            <  1:203246/251‑1 (MQ=255)
                                                                                                                 gTTTCACGGTACCAAAAACACAGGCGGGCAACCAA‑CGGGTCATTCATCTTTTTAAACCCGCAATCG‑ACGCTCTCCCGGGTCAGATTACATTAACGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCCGAGAGTATAGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGAAAAAAATTTTGGGGACCATTTTACCGTTGTCTCAATAAGGCAGATGTgggg                                                                                                                             >  1:693186/1‑249 (MQ=255)
                                                                                                                   ttCACGGTACCAAAACCACAGGCGGG‑AACCAATCCGGTCTTTTATCTTATTAAGCCAGCCATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGGATGAGCAAAGAGCATATCATTGCTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTTAACCTGAAATGTCAGCGGGGGAAAAAAATTATGTGGACCCTTTTTACGTTGACTCAATAAAGCAGATGTggggagc                                                                                                                          >  2:225918/1‑247 (MQ=255)
                                                                                                                     cACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAg                                                                                                                        <  2:92466/251‑1 (MQ=255)
                                                                                                                                    aacggccgg‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCg                                                                                                         <  1:855129/244‑1 (MQ=255)
                                                                                                                                        gCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAAGCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCAGTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTACAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGGCCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCTGACt                                                                                                     <  2:585613/251‑1 (MQ=255)
                                                                                                                                          ggg‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTcc                                                                                                   >  1:292760/1‑251 (MQ=255)
                                                                                                                                                   aTCGGGTCATTCATCTTAGTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTGCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGCTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCTTCGa                                                                                           <  2:871997/251‑1 (MQ=255)
                                                                                                                                                     cGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaa                                                                                          >  2:854978/1‑250 (MQ=255)
                                                                                                                                                     cGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGATCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaa                                                                                          >  1:667960/1‑250 (MQ=255)
                                                                                                                                                          cattcatCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCat                                                                                    <  1:79776/251‑1 (MQ=255)
                                                                                                                                                                          ccAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaaa                                                                                         >  1:501896/1‑230 (MQ=255)
                                                                                                                                                                          ccAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaaa                                                                                         <  2:501896/230‑1 (MQ=255)
                                                                                                                                                                          ccAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATAc                                                                    <  1:444258/251‑1 (MQ=255)
                                                                                                                                                                               aaTCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATg                                                               >  2:318843/1‑251 (MQ=255)
                                                                                                                                                                                 tCG‑ACGATCTCCGGAGTCAGATGACATTCACGACACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGCGTATGGCAGAACAGAAAAACAAAAATGCCCCTTTGTTTTTCAACCTGAAGTGACAGCGAGAGAAAAACATTATGGTGACCATTTTACCGTTGACTCAATAAGGACGATGTGGGACGCAGCGATAAAACGTGCCGGACTCAGCCATCGAAAATCATATCAGTCGAGACATACTTATGcc                                                             >  1:617530/1‑251 (MQ=255)
                                                                                                                                                                                  cG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCct                                                            <  1:423158/251‑1 (MQ=255)
                                                                                                                                                                                          ctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCTCCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGACGAGACATACTTATGCCTGCTgg                                                       <  2:86055/249‑1 (MQ=255)
                                                                                                                                                                                                    aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGt                                                                                                             <  2:44080/185‑1 (MQ=255)
                                                                                                                                                                                                    aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGt                                                                                                             >  1:44080/1‑185 (MQ=255)
                                                                                                                                                                                                      aTGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTgg                                         >  1:108452/1‑251 (MQ=255)
                                                                                                                                                                                                           aTTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTa                                    <  2:576516/251‑1 (MQ=255)
                                                                                                                                                                                                             tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAc                                  >  1:486112/1‑251 (MQ=255)
                                                                                                                                                                                                              aaCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAc                                  <  2:499923/250‑1 (MQ=255)
                                                                                                                                                                                                               aCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAccc                                <  1:723743/251‑1 (MQ=255)
                                                                                                                                                                                                                   gaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCa                            <  2:509853/251‑1 (MQ=255)
                                                                                                                                                                                                                            aaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAACATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCATGCTGGTCGCTGACAGCTGGAGCTAACCCGGCATTAATAGc                    >  2:437739/1‑250 (MQ=255)
                                                                                                                                                                                                                              gagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGTAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCttttttatcaac                 >  2:108423/1‑243 (MQ=255)
                                                                                                                                                                                                                                         tGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTTGTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCATATGTGGGACGCAGCGATAAAACGTTCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCTGCATTTATAGCAAACCAGATGGGc       <  2:593730/250‑1 (MQ=255)
                                                                                                                                                                                                                                          gATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGcc      <  1:284650/250‑1 (MQ=255)
                                                                                                                                                                                                                                              ttCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGTCTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACACGGCATTTATAGCAAACCAGATGGGCCCTGc  >  2:381776/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
GTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652090‑1652575

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