Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,518,356 A→G intergenic (‑121/‑16) fecI ← / → ykiE RNA polymerase sigma factor FecI/protein YkiE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,518,3560AG98.0% 175.1 / ‑6.2 51intergenic (‑121/‑16)fecI/ykiERNA polymerase sigma factor FecI/protein YkiE
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (27/23);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTG  >  NC_000913/4518111‑4518599
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
gCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAagg                                                                                                                                                                                                                                                >  1:82053/1‑251 (MQ=255)
gCAATGTCATCTGCATCAAAAGCAGACTGGAGGTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAagg                                                                                                                                                                                                                                                <  1:105900/251‑1 (MQ=255)
  aaTGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAaggta                                                                                                                                                                                                                                              <  1:45713/251‑1 (MQ=255)
       cATCTGCTTCAAAAGCTGACTGGAGTTTGCGCGTCAGCCAACTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCTGTAAGACTTGTGGAGTGGACTAAAAGTTAATTTTACTGTAGTCATATTAATTTGGAAATAGTTATCATTACAAATGGAATTAAAATTTGTTCCAATTATCTATTTTAACATGGTTTTTTTTCTAGTGTTGTAaggtaggttt                                                                                                                                                                                                                                         >  2:29770/1‑249 (MQ=255)
       cATCGGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAAGTGTTTACAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGtat                                                                                                                                                                                                                                         <  2:498930/251‑1 (MQ=255)
          cTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAaggtaggt                                                                                                                                                                                                                                           >  1:29519/1‑246 (MQ=255)
          cTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAaggtaggt                                                                                                                                                                                                                                           <  2:29519/246‑1 (MQ=255)
                cAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATggg                                                                                                                                                                                                                                 <  1:6974/250‑1 (MQ=255)
                            ggAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCctct                                                                                                                                                                                                                      <  2:328115/249‑1 (MQ=255)
                                                    tCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTCGCATGTGTCGGTCAGAAt                                                                                                                                                                                             <  2:418094/250‑1 (MQ=255)
                                                    tCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATg                                                                                                                                                                                            <  2:206266/250‑1 (MQ=255)
                                                     cAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCCGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATg                                                                                                                                                                                            >  2:196708/1‑250 (MQ=255)
                                                                   gtgtgCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTTGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTTGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGTAATGTGTCGGTCAGAATGACTCACGAtgtgtt                                                                                                                                                                              <  2:289977/250‑3 (MQ=255)
                                                                       gCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGc                                                                                                                                                                          <  1:196708/250‑1 (MQ=255)
                                                                         cATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTAc                                                                                                                                                                       <  1:95234‑M1/251‑2 (MQ=255)
                                                                                cGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGACTGACTCATGATGGGGTCTGCTActcaaggg                                                                                                                                                                >  2:389431‑M1/1‑243 (MQ=255)
                                                                                 gACTCGAACGTTAAGGAAGCTGTGGCAGTGGCGCGGTCAGACATGCGAAGTGCATCAAAAGGTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTATGCTTCTGGCATGTGTCGGTCAGAATGACTCATGTTGTGGTCTGCTActcaaggtg                                                                                                                                                               >  1:447306‑M1/1‑242 (MQ=255)
                                                                                 gACTCGAACGTTAAGGAAGCTGTGGAAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACCCTTGTTTTGGTTTGCTActccagggg                                                                                                                                                               >  2:22583‑M1/1‑242 (MQ=255)
                                                                                      gAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaatt                                                                                                                                                          >  1:247255‑M1/1‑237 (MQ=255)
                                                                                                 aaGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagt                                                                                                                                                >  2:15316‑M1/1‑226 (MQ=255)
                                                                                                               gcgcGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacg                                                                                                                                 <  2:201735‑M1/251‑40 (MQ=255)
                                                                                                                         aaatgCGGAGGGCATCAAAAGATAATTATCCCGTAGTCATATAAATATGAGAATAGTTATCATTACAATAGGAAATAAAATTGTTTCCAATAGACATTTTTACCATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTATGCTActcaaggtggaattgtgtctcagtataccgcctgaaacggggattaatc                                                                                                                       >  2:327466‑M1/3‑202 (MQ=255)
                                                                                                                                          aaaGTTAATTATCACGTCGTCATATTAATTTGAGAATGGTTATCATTACAATTAGAAATAAAATTGTTTACAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGTGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGGGCTCTGCTActcaaggtggcattgtgtcgcagtataccgccagaaacggggcttaatactctttaccagttaccc                                                                                                      >  2:214829‑M1/1‑185 (MQ=255)
                                                                                                                                               tAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggct                                                                                                  <  1:15316‑M1/250‑71 (MQ=255)
                                                                                                                                                        cGTATTCAGATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTAGAAAATGGGATGGAGCCTCTGCTGCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGATActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatcctctttacctgttacccggctacggtagcta                                                                                        <  1:38418‑M1/251‑81 (MQ=255)
                                                                                                                                                          gTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctat                                                                                       >  1:17000‑M1/1‑170 (MQ=255)
                                                                                                                                                             gTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatccc                                                                                    <  1:307335‑M1/251‑85 (MQ=255)
                                                                                                                                                                  atTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggcattgtgtcgcagtatacagccagaaacggggattaatactctttacctgttacccggctacggtagctagcccagttt                                                                               >  2:357363‑M1/1‑162 (MQ=255)
                                                                                                                                                                    tAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttga                                                                              >  1:110992‑M1/1‑160 (MQ=255)
                                                                                                                                                                        atGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttga                                                                              >  1:164224‑M1/1‑156 (MQ=255)
                                                                                                                                                                        atGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggccacggtagctatcccagttgaatatc                                                                         >  2:76352‑M1/1‑156 (MQ=255)
                                                                                                                                                                          gagaATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggttacggtagctatcccagttgaagatca                                                                        <  1:67376‑M1/250‑97 (MQ=255)
                                                                                                                                                                           agaATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATCGCCATTTTTATCATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGTAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttaccctgctacggtagctatcccagttgaagatcacc                                                                      >  2:42872‑M1/1‑153 (MQ=255)
                                                                                                                                                                            gaATGGTTATCATTACAATTGGAAATAAAATNGTTTCCAATAGACATTTTTAGCATTTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGACAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtatgccgccagaagcggggattaatactctttacctgttacccggctacggtagctagcccaggtgaagatcaccc                                                                     <  2:233928‑M1/251‑100 (MQ=255)
                                                                                                                                                                                    aTCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactg                                                             >  2:432382‑M1/1‑144 (MQ=255)
                                                                                                                                                                                                            gTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccgcaaactgttacccaggcgcatacggtccatc                                     <  2:66106‑M1/251‑132 (MQ=255)
                                                                                                                                                                                                            gTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatc                                     <  2:82053‑M1/251‑132 (MQ=255)
                                                                                                                                                                                                              ttCCAATAGACATTTTTAACATGTTGTTTTTCTAATTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatca                                    <  1:294518‑M1/250‑133 (MQ=255)
                                                                                                                                                                                                               tCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTATGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGGCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcac                                   <  2:17000‑M1/250‑134 (MQ=255)
                                                                                                                                                                                                                ccAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgc                                 >  2:332222‑M1/1‑116 (MQ=255)
                                                                                                                                                                                                                     aGACATTTTTAACATGTTGTTTTTCTAAGTGTGGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCTGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttc                             <  1:22583‑M1/250‑140 (MQ=255)
                                                                                                                                                                                                                           ttttAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagc                      >  1:166906‑M1/1‑105 (MQ=255)
                                                                                                                                                                                                                               aaCATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATggg                                                                                                                                                                                                                                 >  2:242728/1‑44 (MQ=255)
                                                                                                                                                                                                                               aaCATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATggg                                                                                                                                                                                                                                 <  1:242728/44‑1 (MQ=255)
                                                                                                                                                                                                                                aCATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaactt                 >  2:86373‑M1/1‑100 (MQ=255)
                                                                                                                                                                                                                                   tgttgtTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggt              >  2:162622‑M1/1‑97 (MQ=255)
                                                                                                                                                                                                                                    gttgttTTTCTAAGTGTTGTAAGGTAGGTATAAAATTGGATGGAGCCTCTGCATCTGGCATGTGGCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccatttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggt              <  1:29770‑M1/250‑155 (MQ=255)
                                                                                                                                                                                                                                      ttgttTTTCTAAGAGTTGTAAGGAAGGGATAAAATGGGATGGTGCCTCTGCTTCTGGCCTGTGTCGGGCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagcaacggggattaatactctttacctgttacccggctacggtagctatcctagttgaagatcacccacaaactgttacccaggcgcatacggtccatcaagcgttcgcccagccacttggtcc            >  2:453364‑M1/1‑95 (MQ=255)
                                                                                                                                                                                                                                       tgttTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcat           >  1:113169‑M1/1‑94 (MQ=255)
                                                                                                                                                                                                                                         tttttCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggcttaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcctct         >  1:199561‑M1/1‑92 (MQ=255)
                                                                                                                                                                                                                                                 aGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccat  >  1:3944‑M1/1‑84 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
GCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTG  >  NC_000913/4518111‑4518599

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