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.

CGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATT  >  NC_000913/2018012‑2018495
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
cgcCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGTATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCAGTGCTCAGCAT‑CGCCAGCAACTT‑AATCAGTTGTCTCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAaaca                                                                                                                                                                                                                                            <  1:47190/251‑1 (MQ=255)
  ccggctgGGTGAAATGCGTCGCGGATGTCAGCAGTCGGAAGAACAGCTCAAAATGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCGGGAAAAA‑GCCATTACTCAGCAG‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTGGACATTGCCCTGAACGGTTTGCGAGAAAAAAaacaac                                                                                                                                                                                                                                          <  2:205696/248‑1 (MQ=255)
     gcttggTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGt                                                                                                                                                                                                                                        <  2:180796/246‑1 (MQ=255)
                                                   aaaaTGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACt                                                                                                                                                                                        >  1:170569/1‑251 (MQ=255)
                                                           gattgattATCAGATTGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AAGCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAAGCCACGGCCGCACTTCTTGCAg                                                                                                                                                                                <  2:445715/251‑1 (MQ=255)
                                                                        aatgaatATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCg                                                                                                                                                                    >  1:118368/1‑250 (MQ=255)
                                                                          tgaatATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc                                                                                                                                                                 >  1:102365/1‑251 (MQ=255)
                                                                           aaatATCGCAATACCCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCACCTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTA‑AATCAGGGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa                                                                                                                                                                <  1:155185/250‑1 (MQ=255)
                                                                                    aaTAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCACCTATCAGCAGTTTATCCAGAAGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCACCTT‑AATCAGTGGACGCAGACAGTTGACATTGCCCTGAACAGTTGGCGAGAAACAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAACACCGCCTCGATCAGaaaa                                                                                                                                                           >  1:30149/1‑247 (MQ=255)
                                                                                     aTAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGACAAAgatg                                                                                                                                                      >  2:163668/1‑251 (MQ=255)
                                                                                      tAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatgg                                                                                                                                                     >  1:382139/1‑251 (MQ=255)
                                                                                           ccAACAGCGATAGGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGGGGACGCAGAAAGTTGACATGGCCCGGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGCCACTGCAGGAACGGAAATCCAGGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAg                                                                                                                                                <  2:28086/250‑1 (MQ=255)
                                                                                               cggCGATATGAGTGCCGGGAAAACCAGCAACCGCTGGATCAACTATCAGCGGTTTATCCAGACGCTGGAAAAA‑GCCATTCCCCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATGGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAACCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCg                                                                                                                                            <  2:232090/249‑1 (MQ=255)
                                                                                                aGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGc                                                                                                                                           >  1:392754/1‑251 (MQ=255)
                                                                                                 gCGATATGAGTGCCGGGATACCCAGCAACCGCTGGATCAACTACCAGCAGTTTATCCAGACGATGAAAAAA‑TCCATTACTCAGCAG‑CGCCAGCAACTA‑AACCAGTGGCCGCAGAAAGTTGACTTTGCCCTGAACAGTTGGCGAAAAAAAAAACATCGTTTCCAGGCCTGG‑CAGCCACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGcc                                                                                                                                          <  1:194350/251‑1 (MQ=255)
                                                                                                   gATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGCGGACGCAGAAAGTTGACATTACCCTGAACAGTTGGCGAGCAACAAAAAAACGTTTGCCGACCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGACAAAGATGGCTGCGTTCGccc                                                                                                                                         >  1:299640/1‑250 (MQ=255)
                                                                                                            gCCGGGATAACCAGCAACCGCTGGACCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTCCTCAGCAT‑CGCCAGAAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTGGGCGAAAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCATTCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGcgccg                                                                                                                                <  1:49112/250‑1 (MQ=255)
                                                                                                               gggAAAACCAACAACCGCGGGTTAAACTACCAGCAGTTTACCCAGCCGCTGGAAAAA‑GCCATTCCCCCGAAT‑CGCCAGCACCTT‑AATCAGGGGACGCAAAAAGTTGACATTGCCATAAACATTTGGCGAGAAAAAAAACAAGTTTTGCAGGCCTGG‑CAGACACTCCAGGAACGGCAATCCACGGGGGCCCGTCTTGCAGAAACCCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCAt                                                                                                                            <  2:341356/251‑1 (MQ=255)
                                                                                                                                 gggATCAACTATCAGCAGTTTATCCAGACG‑TGGAAAAAAGCCATTACCCACCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCGGG‑CAGACACTGCAGGAACGGCAACCCACGGGGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGAtt                                                                                                          <  1:214335/250‑1 (MQ=255)
                                                                                                                                  ggATCAACTATCGGAAGTTTACCCAGACGCGGAAAAAA‑GCCATTCCCCAG‑ATACGCCAGCAACTA‑AAGAAGTGGACGCAGAAAGTTGACATTGCCCGGACAAGTTGGGGAGAAAAAAAACAACGTTTGCGGCCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACGTCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGAt                                                                                                           <  2:209577/249‑1 (MQ=255)
                                                                                                                                   gatcaacaaccagaggTTTATCCAACCGCTGGAAAAA‑GCAATTCCCCAGCAT‑CGCCAGAAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAAAAGTTGGCGAAAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACGGCAGGAACGGCAACCCACGGCGGCACTGCTTGCAGAAACCCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCg                                                                                                        <  2:392326/236‑1 (MQ=255)
                                                                                                                                               gcagTTTACCCATACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCGGAACAGTTGGCGAAAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGCCATCCCACGTCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCtt                                                                                            <  1:287307/251‑1 (MQ=255)
                                                                                                                                                    ttATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATGGCCCGGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGCCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGAtt                                                                                        <  1:289819/250‑1 (MQ=255)
                                                                                                                                                       cccAGACGCGGGAAAAA‑GCCATTACTCAGCAG‑CGCCAG‑AACTTGAACCAGTGGACGCAGAAAGTTGACATTGCCCGGAACAGTGGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACGGCAGGACCGGCAATCCACGGGGGCCCGGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTAccg                                                                                    <  2:110904/250‑1 (MQ=255)
                                                                                                                                                            aCGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAAAAAGTTAACATTGCCCTGAACAGTTGGCGAAAAAAAAAACAACTTTTGCAGGCCTGG‑CAGACCCTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAAcc                                                                                                                   <  1:289507/215‑1 (MQ=255)
                                                                                                                                                                      aa‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCCCTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGAcacc                                                                      >  1:297013/1‑250 (MQ=255)
                                                                                                                                                                           cATTACTGAGCTT‑CGCCAGAAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGTCCTGC‑CAGACACTGCAGGACCGGCAATCCACGGCGGCACTGCTTGCAGAACACCGCATCGATCAGAAAAAGATGGATGAGTTCGCCAAGCACGCCGCCCTGAGGCAACCTGAATGAATCGCTTAGCGCCCTGGATTACCTCCGACGTTGAccccaccac                                                                 >  2:50395/1‑251 (MQ=255)
                                                                                                                                                                              tACTCAGAAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACAtt                                                              >  2:157649/1‑251 (MQ=255)
                                                                                                                                                                                                      aaccagtgcaccaagaaattgaaCATTCCCCTGAACAGTGGGCAAAAAAAAAAACAACTTTTGCAGGCCGGG‑CGGCCTCCGCAGGCACGGCAACCCACGCCGCCATTGTTTGCAAAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAACCCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGAt                                        <  1:235929/229‑1 (MQ=255)
                                                                                                                                                                                                          aGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGATCGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGAtgctg                                    >  1:105493/1‑251 (MQ=255)
                                                                                                                                                                                                            tGGACGCAGAACGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAACACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGAtgctgc                                   >  2:67777/1‑250 (MQ=255)
                                                                                                                                                                                                             ggACGCAAAAAGTTGACATTGCCAGAAACAGTGGGAGAAAAAAAAAACAACGTGGGCAGGCCTGG‑CAGACACGGCGGGAACGGAAACCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGcac                                 <  2:170569/251‑1 (MQ=255)
                                                                                                                                                                                                                cGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGACGACGTTg                                                                           >  2:242525/1‑206 (MQ=255)
                                                                                                                                                                                                                cGCAAAAAGTTGACATTCCCCGAAACAGTGGGCGAAAAAAAAACCAACGTTTGCAGGCCGGG‑CAGACACTGCAGGAACGGCAACCCACGGCGGCACTGCTTGCAAAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTg                                                                           <  1:242525/206‑1 (MQ=255)
                                                                                                                                                                                                                   aGAAAGTTGACATTGCCCTGAACAGTTGGCAAAAAAAAAAACAACTTTTGCAGCCCGGG‑CAGACCCTGCAGGACCGGAAACCCACGCGGGCACTGTTTCCAGAAACCCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGAtt                           <  1:347028/251‑1 (MQ=255)
                                                                                                                                                                                                                          tGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTcgc                    >  2:440477/1‑251 (MQ=255)
                                                                                                                                                                                                                           gactttGCCCTGAACATTGGGCAAAAAAAAAAACAAC‑TTTGCCAGCCTGGGCACACACGGCAGGACCTCCAACCCACGGCGGCACTTCTTCCAAAAACCCCCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTcgcg                   <  1:373895/247‑1 (MQ=255)
                                                                                                                                                                                                                                         cAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCCCTGCTTGCAGAAAACCGCCTCGATCACAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTAACACCACCACATTGCCTAGCGGCAAAACCAGCGATGCTGCACAAGATTTTCTCGCGCTGTTGAGCGatgc     >  1:15071/1‑248 (MQ=255)
                                                                                                                                                                                                                                           gTTGGCGAAAAAAAAAACAACTTTTGCGGCCCTGG‑CAGACACTGCAGGAACGGAAATCCACGGCGGCCCTGCTTGCAAAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCAt   <  1:293766/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ttGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCCGAAAACCGCCTCGATCCGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGACAAAGACAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCAtt  >  1:159968/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
CGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCTGATTGATTATCAGAATGAATATCGCAATAACCTCAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAA‑GCCATTACTCAGCAT‑CGCCAGCAACTT‑AATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGG‑CAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATT  >  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: ‑