Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,963,203 G→A 8.2% G616G (GGC→GGT)  flhA ← putative flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,963,2030GA 8.2% 148.7 / 4.2 49G616G (GGC→GGT) flhAputative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base G (21/24);  new base A (2/2);  total (23/26)
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.68e-01

CGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTGCAGCACTTTATGCAGTGTGGTGAGCGTGACGACGCCAGGAACGAGATCTTCCGT  >  NC_000913/1962921‑1963481
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                      
cggcggCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGc                                                                                                                                                                                                                                                                                  <  2:208281/289‑1 (MQ=255)
      cTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAgcgc                                                                                                                                                                                                                                                                             <  2:24917/288‑1 (MQ=255)
                       aaCAGTAAATTTGCGAATAAGGTTCCAATTATTTGCCGCCAATTGTGGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTc                                                                                                                                                                                                                                                           <  2:500813/289‑1 (MQ=255)
                                      aataaGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGTATAGCGCTTCCTGAGTTTGCGCCAGTACACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCCCCTGTAGTCACCAACGTTTCAGCCGTGTATCgcg                                                                                                                                                                                                                                            >  2:122511/1‑287 (MQ=255)
                                        taaGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAgg                                                                                                                                                                                                                                           >  1:88642/1‑288 (MQ=255)
                                              tctcATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAATCTGCGGCGCAGGAAGCGAGTCAATAATGTTCGCAGCGCGCGGTTCACCAACAATACTGGCGGCGCACCCAGCACCTTCTGTCGGGCTAGCGCTTCCTGAGTTTGCTCCAAGTAACGAACCGCCAGCCACTGACTCCGGT‑CTTCCCCGCCCTGCAGCTCCTGTACTACACAACGTTCCCGGGGGCGTTCGAGGCCCATac                                                                                                                                                                                                                                    >  2:433535/1‑288 (MQ=255)
                                              tctcATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATa                                                                                                                                                                                                                                     >  1:384640/1‑288 (MQ=255)
                                               ctcATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCCGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAAc                                                                                                                                                                                                                                   >  2:32683/1‑289 (MQ=255)
                                                 cATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTc                                                                                                                                                                                                                                                           >  2:430109/1‑263 (MQ=255)
                                                 cATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTc                                                                                                                                                                                                                                                           <  1:430109/263‑1 (MQ=255)
                                                 cATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACAt                                                                                                                                                                                                                                 >  1:204694/1‑289 (MQ=255)
                                                              cAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCttt                                                                                                                                                                                                                    <  2:137659/289‑1 (MQ=255)
                                                                        cGTCATGCGGTTATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAAc                                                                                                                                                                                                          <  1:166650/289‑1 (MQ=255)
                                                                        cGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAAc                                                                                                                                                                                                          <  2:292818/289‑1 (MQ=255)
                                                                        cGTAATGCGGGTATGTCGGTTATCGGACAGTTCCAGATTCGACAGGACCACTAACGGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAAc                                                                                                                                                                                                          <  2:383021/289‑1 (MQ=255)
                                                                                    atGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTa                                                                                                                                                                                              >  2:477739/1‑289 (MQ=255)
                                                                                           ttacaaGACAGTTCCAGATTCGACAGGACCACTAACGGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGcc                                                                                                                                                                                        <  2:204694/283‑1 (MQ=255)
                                                                                                                gACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCCcggccccg                                                                                                                                                                   >  2:423951/1‑282 (MQ=255)
                                                                                                                gACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCacg                                                                                                                                                                   >  1:284896/1‑288 (MQ=255)
                                                                                                                gACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCACCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCacg                                                                                                                                                                   >  1:107053/1‑288 (MQ=255)
                                                                                                                       ccACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGt                                                                                                                                                           <  1:293632/289‑1 (MQ=255)
                                                                                                                          cTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTaa                                                                                                                                                        >  1:409190/1‑289 (MQ=255)
                                                                                                                          cTAACGGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTaa                                                                                                                                                        <  1:105200/289‑1 (MQ=255)
                                                                                                                               tGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCAt                                                                                                                                                   >  1:197643/1‑289 (MQ=255)
                                                                                                                                           cggcgCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCttt                                                                                                                                       >  1:295764/1‑289 (MQ=255)
                                                                                                                                                                 atGGTCGCAGCGCGGGGTTCACCATCAATACTAGCGGCGACCCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAg                                                                                                                 <  2:384640/289‑1 (MQ=255)
                                                                                                                                                                    gTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGt                                                                                                              >  1:256414/1‑289 (MQ=255)
                                                                                                                                                                     tagCAGCGCGTGGTTAACCAACAATACATGCGGCGCACCCAGCAGCTCCTGACGGGATAGCGCTTCCTGATTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGt                                                                                                              <  1:268729/286‑1 (MQ=255)
                                                                                                                                                                              gTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATg                                                                                                    >  1:248210/1‑289 (MQ=255)
                                                                                                                                                                                         acaaTACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCatat                                                                                          >  1:9231/1‑288 (MQ=255)
                                                                                                                                                                                          caaTACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATcg                                                                                        <  2:284896/289‑1 (MQ=255)
                                                                                                                                                                                            aTACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATcgcg                                                                                      >  1:245152/1‑289 (MQ=255)
                                                                                                                                                                                                 ggcggcGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCg                                                                                 >  2:288546/1‑289 (MQ=255)
                                                                                                                                                                                                   cggaggaCCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGGGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCgg                                                                                <  2:503761/282‑1 (MQ=255)
                                                                                                                                                                                                                         cGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGAtttt                                                         >  1:358431/1‑289 (MQ=255)
                                                                                                                                                                                                                             aTAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTgca                                                      <  1:32683/288‑1 (MQ=255)
                                                                                                                                                                                                                             aTAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCACCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTgca                                                      <  2:107053/288‑1 (MQ=255)
                                                                                                                                                                                                                             aTAGCGCTTCCTGAGGTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTgca                                                      <  1:71612/288‑1 (MQ=255)
                                                                                                                                                                                                                                  gcTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTGCAGCACt                                                 <  1:306409/288‑1 (MQ=255)
                                                                                                                                                                                                                                        tGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTGCAGCACTTTATGCa                                          <  1:288546/289‑1 (MQ=255)
                                                                                                                                                                                                                                             ttGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTGCAGCACTTTATGCAGtgtg                                     <  2:358431/289‑1 (MQ=255)
                                                                                                                                                                                                                                                cgcCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTTATCGAGGAGATTTTGCAGCCCCTTATGCAGTgggggt                                  >  2:352508/1‑287 (MQ=255)
                                                                                                                                                                                                                                                               cGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTGCAGCACTTTATGCAGTGTGGTGAGCGTGACGACGCCa                   <  1:423951/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ggCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTGCAGCACTTTATGCAGTGTGGTGAGCGTGACGACGCCAGGAACGAGAt         <  2:409190/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   cctccCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGGGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTGCAGCACTTTATGCAGTGTGGTGAGCGTGACGACGCCAGGAACGAGATCTTCCGt  <  1:348221/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      cccccACCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTGCAGCACTTTATGCAGTGTGGTGAGCGTGACGACGcc                    >  2:475700/1‑266 (MQ=255)
                                                                                                                                                                                                                                                                                      cccccACCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTGCAGCACTTTATGCAGTGTGGTGAGCGTGACGACGcc                    <  1:475700/266‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         ccGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCtt                                                                                                                                        <  2:447184/147‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         ccGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCtt                                                                                                                                        >  1:447184/1‑147 (MQ=255)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                      
CGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCC‑CTGGCTCCAGTCCTCCCCCGCCCTGCAGCGCCTGTAGTAACAAACGTTCCAGCGGTGTATCGAGGCCAATAACATGGACTTCATCTTTGCCAGGAAACCACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCATCGAGGAGATTTTGCAGCACTTTATGCAGTGTGGTGAGCGTGACGACGCCAGGAACGAGATCTTCCGT  >  NC_000913/1962921‑1963481

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑