Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,018,2520GA27.5% 62.6 / 10.0 40E103K (GAA→AAA) fliJflagellar biosynthesis protein FliJ
Reads supporting (aligned to +/- strand):  ref base G (17/12);  new base A (0/11);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 7.84e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATT  >  NC_000913/2018012‑2018495
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
cgcCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGTATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCAGTGCTCAGCATCGCCAGCAACTTAATCAGTTGTCTCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAaaca                                                                                                                                                                                                                                            <  1:47190/251‑1 (MQ=255)
  ccggctgGGTGAAATGCGTCGCGGATGTCAGCAGTCGGAAGAACAGCTCAAAATGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCGGGAAAAAGCCATTACTCAGCAGCGCCAGCAACTTAATCAGTGGACGCAGAAAGTGGACATTGCCCTGAACGGTTTGCGAGAAAAAAaacaac                                                                                                                                                                                                                                          <  2:205696/248‑1 (MQ=255)
     gcttggTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGt                                                                                                                                                                                                                                        <  2:180796/246‑1 (MQ=255)
                                                   aaaaTGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACt                                                                                                                                                                                        >  1:170569/1‑251 (MQ=255)
                                                           gattgattATCAGATTGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAAGCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAAGCCACGGCCGCACTTCTTGCAg                                                                                                                                                                                <  2:445715/251‑1 (MQ=255)
                                                                        aatgaatATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCg                                                                                                                                                                    >  1:118368/1‑250 (MQ=255)
                                                                          tgaatATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc                                                                                                                                                                 >  1:102365/1‑251 (MQ=255)
                                                                           aaatATCGCAATACCCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCACCTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTAAATCAGGGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa                                                                                                                                                                <  1:155185/250‑1 (MQ=255)
                                                                                    aaTAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCACCTATCAGCAGTTTATCCAGAAGCTGGAAAAAGCCATTACTCAGCATCGCCAGCACCTTAATCAGTGGACGCAGACAGTTGACATTGCCCTGAACAGTTGGCGAGAAACAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAACACCGCCTCGATCAGaaaa                                                                                                                                                           >  1:30149/1‑247 (MQ=255)
                                                                                     aTAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGACAAAgatg                                                                                                                                                      >  2:163668/1‑251 (MQ=255)
                                                                                      tAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatgg                                                                                                                                                     >  1:382139/1‑251 (MQ=255)
                                                                                           ccAACAGCGATAGGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGGGGACGCAGAAAGTTGACATGGCCCGGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGCCACTGCAGGAACGGAAATCCAGGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAg                                                                                                                                                <  2:28086/250‑1 (MQ=255)
                                                                                               cggCGATATGAGTGCCGGGAAAACCAGCAACCGCTGGATCAACTATCAGCGGTTTATCCAGACGCTGGAAAAAGCCATTCCCCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATGGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAACCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCg                                                                                                                                            <  2:232090/249‑1 (MQ=255)
                                                                                                aGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGc                                                                                                                                           >  1:392754/1‑251 (MQ=255)
                                                                                                 gCGATATGAGTGCCGGGATACCCAGCAACCGCTGGATCAACTACCAGCAGTTTATCCAGACGATGAAAAAATCCATTACTCAGCAGCGCCAGCAACTAAACCAGTGGCCGCAGAAAGTTGACTTTGCCCTGAACAGTTGGCGAAAAAAAAAACATCGTTTCCAGGCCTGGCAGCCACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGcc                                                                                                                                          <  1:194350/251‑1 (MQ=255)
                                                                                                   gATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGCGGACGCAGAAAGTTGACATTACCCTGAACAGTTGGCGAGCAACAAAAAAACGTTTGCCGACCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGACAAAGATGGCTGCGTTCGccc                                                                                                                                         >  1:299640/1‑250 (MQ=255)
                                                                                                            gCCGGGATAACCAGCAACCGCTGGACCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTCCTCAGCATCGCCAGAAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTGGGCGAAAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCATTCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGcgccg                                                                                                                                <  1:49112/250‑1 (MQ=255)
                                                                                                               gggAAAACCAACAACCGCGGGTTAAACTACCAGCAGTTTACCCAGCCGCTGGAAAAAGCCATTCCCCCGAATCGCCAGCACCTTAATCAGGGGACGCAAAAAGTTGACATTGCCATAAACATTTGGCGAGAAAAAAAACAAGTTTTGCAGGCCTGGCAGACACTCCAGGAACGGCAATCCACGGGGGCCCGTCTTGCAGAAACCCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCAt                                                                                                                            <  2:341356/251‑1 (MQ=255)
                                                                                                                                 gggATCAACTATCAGCAGTTTATCCAGACGTGGAAAAAAGCCATTACCCACCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCGGGCAGACACTGCAGGAACGGCAACCCACGGGGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGAtt                                                                                                          <  1:214335/250‑1 (MQ=255)
                                                                                                                                  ggATCAACTATCGGAAGTTTACCCAGACGCGGAAAAAAGCCATTCCCCAGATACGCCAGCAACTAAAGAAGTGGACGCAGAAAGTTGACATTGCCCGGACAAGTTGGGGAGAAAAAAAACAACGTTTGCGGCCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACGTCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGAt                                                                                                           <  2:209577/249‑1 (MQ=255)
                                                                                                                                   gatcaacaaccagaggTTTATCCAACCGCTGGAAAAAGCAATTCCCCAGCATCGCCAGAAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAAAAGTTGGCGAAAAAAAAAACAACGTTTGCAGGCCTGGCAGACACGGCAGGAACGGCAACCCACGGCGGCACTGCTTGCAGAAACCCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCg                                                                                                        <  2:392326/236‑1 (MQ=255)
                                                                                                                                               gcagTTTACCCATACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCGGAACAGTTGGCGAAAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGCCATCCCACGTCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCtt                                                                                            <  1:287307/251‑1 (MQ=255)
                                                                                                                                                    ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATGGCCCGGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGCCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGAtt                                                                                        <  1:289819/250‑1 (MQ=255)
                                                                                                                                                       cccAGACGCGGGAAAAAGCCATTACTCAGCAGCGCCAGAACTTGAACCAGTGGACGCAGAAAGTTGACATTGCCCGGAACAGTGGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACGGCAGGACCGGCAATCCACGGGGGCCCGGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTAccg                                                                                    <  2:110904/250‑1 (MQ=255)
                                                                                                                                                            aCGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAAAAAGTTAACATTGCCCTGAACAGTTGGCGAAAAAAAAAACAACTTTTGCAGGCCTGGCAGACCCTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAAcc                                                                                                                   <  1:289507/215‑1 (MQ=255)
                                                                                                                                                                      aaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCCCTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGAcacc                                                                      >  1:297013/1‑250 (MQ=255)
                                                                                                                                                                           cATTACTGAGCTTCGCCAGAAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGTCCTGCCAGACACTGCAGGACCGGCAATCCACGGCGGCACTGCTTGCAGAACACCGCATCGATCAGAAAAAGATGGATGAGTTCGCCAAGCACGCCGCCCTGAGGCAACCTGAATGAATCGCTTAGCGCCCTGGATTACCTCCGACGTTGAccccaccac                                                                 >  2:50395/1‑251 (MQ=255)
                                                                                                                                                                              tACTCAGAATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACAtt                                                              >  2:157649/1‑251 (MQ=255)
                                                                                                                                                                                                      aaccagtgcaccaagaaattgaaCATTCCCCTGAACAGTGGGCAAAAAAAAAAACAACTTTTGCAGGCCGGGCGGCCTCCGCAGGCACGGCAACCCACGCCGCCATTGTTTGCAAAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAACCCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGAt                                        <  1:235929/229‑1 (MQ=255)
                                                                                                                                                                                                          aGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGATCGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGAtgctg                                    >  1:105493/1‑251 (MQ=255)
                                                                                                                                                                                                            tGGACGCAGAACGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAACACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGAtgctgc                                   >  2:67777/1‑250 (MQ=255)
                                                                                                                                                                                                             ggACGCAAAAAGTTGACATTGCCAGAAACAGTGGGAGAAAAAAAAAACAACGTGGGCAGGCCTGGCAGACACGGCGGGAACGGAAACCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGcac                                 <  2:170569/251‑1 (MQ=255)
                                                                                                                                                                                                                cGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGACGACGTTg                                                                           >  2:242525/1‑206 (MQ=255)
                                                                                                                                                                                                                cGCAAAAAGTTGACATTCCCCGAAACAGTGGGCGAAAAAAAAACCAACGTTTGCAGGCCGGGCAGACACTGCAGGAACGGCAACCCACGGCGGCACTGCTTGCAAAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTg                                                                           <  1:242525/206‑1 (MQ=255)
                                                                                                                                                                                                                   aGAAAGTTGACATTGCCCTGAACAGTTGGCAAAAAAAAAAACAACTTTTGCAGCCCGGGCAGACCCTGCAGGACCGGAAACCCACGCGGGCACTGTTTCCAGAAACCCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGAtt                           <  1:347028/251‑1 (MQ=255)
                                                                                                                                                                                                                          tGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTcgc                    >  2:440477/1‑251 (MQ=255)
                                                                                                                                                                                                                           gactttGCCCTGAACATTGGGCAAAAAAAAAAACAACTTTGCCAGCCTGGGCACACACGGCAGGACCTCCAACCCACGGCGGCACTTCTTCCAAAAACCCCCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTcgcg                   <  1:373895/247‑1 (MQ=255)
                                                                                                                                                                                                                                         cAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCCCTGCTTGCAGAAAACCGCCTCGATCACAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTAACACCACCACATTGCCTAGCGGCAAAACCAGCGATGCTGCACAAGATTTTCTCGCGCTGTTGAGCGatgc     >  1:15071/1‑248 (MQ=255)
                                                                                                                                                                                                                                           gTTGGCGAAAAAAAAAACAACTTTTGCGGCCCTGGCAGACACTGCAGGAACGGAAATCCACGGCGGCCCTGCTTGCAAAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCAt   <  1:293766/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ttGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCCGAAAACCGCCTCGATCCGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGACAAAGACAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCAtt  >  1:159968/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
CGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATT  >  NC_000913/2018012‑2018495

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