Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,275,6830TG50.0% 42.2 / 15.4 36T1108P (ACC→CCC) unknownChromosome partition protein MukB
Reads supporting (aligned to +/- strand):  ref base T (7/11);  new base G (18/0);  total (25/11)
Fisher's exact test for biased strand distribution p-value = 1.06e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTACGCAACGCACCAAGGGCTTTGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTGAAAGAACTTTGCAATTGAATCAACACTTGATTGTATTGTGCATATTGCGCTTGAACCTGACGCAATTGCTCGCGCTGTGTTTCGCGT  >  1/1275454‑1275933
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
gTACGCAACGCACCAAGGGCTTTGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑ATTTTTCGCTTTCCCTTCAATTGGGGGGGGGTTGTTTTTCATTGTAAg                                                                                                                                                                                                                                         >  1:28876/1‑251 (MQ=255)
gTACGCAACGCACCAAGGGCTTTGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑ATTTTTCGCTTCGCTTTCAA‑TGAGGGGGGGGTGTTTTTTAATGTAAGa                                                                                                                                                                                                                                        >  1:524729/1‑251 (MQ=255)
gTACGCAACGCACCAAGGGCTTTGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑AATTTTCGCTTCGCTTTTTA‑TTAGGGGGAGTTGTTTTTCTATGTaaaa                                                                                                                                                                                                                                        >  1:529489/1‑249 (MQ=255)
 tACGCACCCCACCAAGGGCTTTGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCCCCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGaa                                                                                                                                                                                                                                       <  1:398264/251‑1 (MQ=255)
 tACGCAACGCACCAAGGGCTTTGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGaa                                                                                                                                                                                                                                       <  1:335454/251‑1 (MQ=255)
 tACGCAACGCACCAAGGGCTTTGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGaa                                                                                                                                                                                                                                       <  1:94424/251‑1 (MQ=255)
       acgcaCCAAGGGCTTTGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑ATTTTTCGCTTCGCTTTCTA‑TGAGGGGGGGGTGTTTTTTAATGTAAAAAAGGCa                                                                                                                                                                                                                                  >  1:373038/1‑250 (MQ=255)
          caCCAAGGGCTTTGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGGAGAATTTTCC‑CCTCGCTTTTCA‑TGAGGGGGGGGTGTTTTTTAATGTAAAAAAGGCACTgg                                                                                                                                                                                                                              >  1:596137/1‑250 (MQ=255)
            ccAAGGGCTTTGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACg                                                                                                                                                                                                                            <  1:90286/251‑1 (MQ=255)
                      tGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCCCCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGaa                                                                                                                                                                                                                  <  1:321529/251‑1 (MQ=255)
                                      tAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAAc                                                                                                                                                                                                   <  1:89390/250‑1 (MQ=255)
                                         cTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTc                                                                                                                                                                                                 <  1:409917/249‑1 (MQ=255)
                                                  ggATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTc                                                                                                                                                                                                 <  1:256919/240‑1 (MQ=255)
                                                          atGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑ATTTTTCGCTTCGCTTTCAA‑TGAGGGGGAGTTGTTTTTCAATGTAAGAAAGGCACTGa                                                                                                                                                                                                                              >  1:642689/1‑203 (MQ=255)
                                                                       acgaTTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCCCCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCgcgc                                                                                                                                                                 <  1:85185/251‑1 (MQ=255)
                                                                                 cGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑ATTTTCCCCTTCGCTTTCCA‑TTGGGGGGAGTTGTTTTTCAATGTAAAAACGGCACCGACGATTAGTGGAAAGTTGTTGATAAAACTCATCCCCACGGATTTTTGCACGTTTTTCCGCGCCTTTGTCgg                                                                                                                                                        >  1:628760/1‑250 (MQ=255)
                                                                                             tCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGGGGGGTGGTTTTTAATGTAAGAACGGCCCTGACGATTAGGGGAAAGTTGTTGTTACAACTCATCCCGACCGATTTTTTCACGTTCTTTCGCGCCTTCGTCGGCACGTACAcccca                                                                                                                                           >  1:410297/1‑249 (MQ=255)
                                                                                             tCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑AATTTCCGCTTCGCTTTTAA‑TGAGGGGGGGTTGTTTTTTTATTTAAGAAAGGCACCGACGATTTGTGGGAAATTGTTTATACAAATCATCCCCACGGATTTTTTCACGGTCTTTCGCGCCTTCGTCGGCACGTACAcccca                                                                                                                                           >  1:135500/1‑249 (MQ=255)
                                                                                                            aaaCGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCCTCGCTTTCAA‑TGAGGGGGAGGTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGGTGTTGATAAAACTCATCCCGACGGATTTTTGCACGGTCTTCCGCGCCTTCGTCCGCACGTACACCCAATTCGCCGAtttt                                                                                                                               >  1:538732/1‑247 (MQ=255)
                                                                                                                                gCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑AATTTTCGCCTCGCTTTCAA‑TGAGGGGGAGTTGTTTTTCAATGTAAGGACGGCACTGACGATTAGTGGGAAGTTGTTGCTACAAATCATTCCGACGGAATTTTTCACGTTCTTCCTCGCCTTCGTCGGGACGTACACCCCATTCGCCGATTTTCTGCATTAATTCCTGCCACA‑ttt                                                                                                       >  1:442510/1‑251 (MQ=255)
                                                                                                                                   cATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGGGAGTTGTTTTTTAATGTAAGAACGGCACTGACGATTAGTGGGAAGTTGTTGATACAACTCATCCCGACGGATTTTTTCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGt                                                                                                    >  1:453211/1‑251 (MQ=255)
                                                                                                                                   cATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑ATTTTTCCCCTCGCTTTCAA‑TTCGGGGGGGTTGTTTTTTAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTTTTACAAATCATTCCCACGGGTTTTTTCACGTTTTTTCGCCCCTTCGTCGGCACGTAACCCCCATTCGCCGAGTTTCTGcctttatttcttctaccttttgt                                                                                                     >  1:629425/1‑229 (MQ=255)
                                                                                                                                                                gggTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCTCTTTCAA‑TGAGGGTGAGTTGTTTTTTAATGAAAGAACGGCAACGGCGCTTAGTGGAAAGTTGTTGAGACAACTCCACCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCTTCGGCACGTACAACCAATTTGCCGCTTTTCTGCCTTAAATTCTGCAAAA‑TTTGGGTTTTACTGTTGTAAGAACTTTGCAATTGa                                                                       >  1:484545/1‑251 (MQ=255)
                                                                                                                                                                gggTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑ATTTTTTGCTTCGCTTTCAA‑TGAGGGGGAGGTGTTTTTTAATGTAAGAACGGCACTGACGATTAGTGGGAAGGTGTTGATACAACTCATCCCGACGGGTTTTTTTACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTTGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTTAAAGAACTTTGCAATTGa                                                                       >  1:333078/1‑251 (MQ=255)
                                                                                                                                                                     ttGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGGGAGGTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGGAAATTGTTGATACAACTCATCCCGACGGATTTTTTCACGTTCTTCCGCGCCCTCGTCGGCACGTACACCCCATTCGCCGATTTCCTGCATTAATTCCTGCAACATTTTGGTTTTT‑CTTTTGAAAGAACTTTTCAATTTAATCaa                                                                  >  1:626601/1‑251 (MQ=255)
                                                                                                                                                                     ttGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑ATTTTCCGCTTCGCTTTTAA‑TGAGGGGGAGGTGTTTTTTAATGTAAGAACGGCACTGACGATTAGTGGAAAGGTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTTATTTCTGCAACA‑TTTGGGTTTTACTTTTGAAAGAACTTTTCAATTTAATCaa                                                                  >  1:510685/1‑251 (MQ=255)
                                                                                                                                                                     ttGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑ATTTTCCGCTTCGCTTTTAA‑TGAGGGGGAGGTGTTTTTTAATGTAAGAACGGCACTGACGATTAGTGGAAAATTGTTGATACAACTCATCCCCACGGATTTTTTCACGTTCTTCCCCGCCTTCGTCGGCACGTACACCCAAATCGCCCATTTCCCGCATTAATTCCTGCAACA‑TTTGGGTTTTACTGTTGAAAGAACTTTGCAATTGAATCaa                                                                  >  1:399049/1‑251 (MQ=255)
                                                                                                                                                                       gTAGTCCCGTTCCGCTTTACGAATGCGGCGGGGAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGGGGGGTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGGAAGGTGTTGATACAACTCATCCCGACGGATTTTTTCACGTTCTTCCGCGCCCTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTTCTGTTGAAAGAACTTTGCAATTGAATCAaca                                                                >  1:368457/1‑251 (MQ=255)
                                                                                                                                                                                ttCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTGAAAGAACTTTGCAATTGAATCAACACTTGattgt                                                       <  1:59022/251‑1 (MQ=255)
                                                                                                                                                                                       ttACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTGAAAGAACTTTGCAATTGAATCAACACTTGATTGTATTGTGc                                                <  1:200348/251‑1 (MQ=255)
                                                                                                                                                                                                            aTTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTGAAAGAACTTTGCAATTGAATCAACACTTGATTGTATTGTGCATATTGCGCTTGAACCTGa                             >  1:277082/1‑250 (MQ=255)
                                                                                                                                                                                                                         cttTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTGAAAGAACTTTGCAATTGAATCAACACTTGATTGTATTGTGCATATTg                                          >  1:624104/1‑224 (MQ=255)
                                                                                                                                                                                                                             cAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTGAAAGAACTTTGCAATTGAATCAACACTTGATTGTATTGTGCATATTg                                          >  1:69226/1‑220 (MQ=255)
                                                                                                                                                                                                                             cAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTGAAAGAACTTTGCAATTGAATCAACACTTGATTGTATTGTGCATATTGCGCTTGAACCTGACGCAATTGCTcgcg               >  1:167740/1‑247 (MQ=255)
                                                                                                                                                                                                                                  gAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTGAAAGAACTTTGCAATTGAATCAACACTTGATTGTATTGTGCATATTg                                          >  1:287316/1‑216 (MQ=255)
                                                                                                                                                                                                                                   aGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTGAAAGAACTTTGCAATTGAATCAACACTTGATTGTATTGTGCATATTGCGCTTGAACCTGACGCAATTGCTCGCGCTGTGttt       >  1:426454/1‑250 (MQ=255)
                                                                                                                                                                                                                                       tGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTGAAAGAACTTTGCAATTGAATCAACACTTGATTGTATTGTGCATATTGCGCTTGAACCTGACGCAATTGCTCGCGCTGTGTTTCGCgt  <  1:98887/251‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
GTACGCAACGCACCAAGGGCTTTGTCCGACATGGAACGTAGCTCATCGGCGGATAAATATGCCAATTCACGACGATTTAAGCGTTTTTCCACATCAGAATTGCGAGATAAACGCAGCACCACGCACCAGCTCATTTTAGCCGCCGTTACCAATTCGCGCTGGGTTTTGTAGTCCCGTTCCGCTTTACGAATGCGGCGGGTAAG‑ATTTTCCGCTTCGCTTTCAA‑TGAGGGTGAGTTGTTTTTCAATGTAAGAACGGCACTGACGATTAGTGGAAAGTTGTTGATACAACTCATCCCGACGGATTTTTGCACGTTCTTCCGCGCCTTCGTCGGCACGTACACCCAATTCGCCGATTTCCTGCATTAATTCCTGCAACA‑TTTGGTTTTTACTGTTGAAAGAACTTTGCAATTGAATCAACACTTGATTGTATTGTGCATATTGCGCTTGAACCTGACGCAATTGCTCGCGCTGTGTTTCGCGT  >  1/1275454‑1275933

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