Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6610TG40.0% 49.2 / 16.8 56pseudogene (943/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base T (10/23);  new base G (20/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 7.78e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCT  >  NC_000913/3862513‑3862807
                                                                                                                                                       |                                                                                                                                                      
ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCCGTTCCGTAATTCCCACCAGCATGGC‑GGTTAAGGtt                                                                                                                                                       <  2:515570/150‑1 (MQ=255)
  gCCCATGGTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGGACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCAACAGCATGGCGGGTTAAGGtgcg                                                                                                                                                     >  2:330836/1‑148 (MQ=255)
     cATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGg‑gggg                                                                                                                                                >  1:524431/1‑150 (MQ=255)
         ttGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGC‑ggg                                                                                                                                                 >  1:292656/1‑146 (MQ=255)
         ttGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGGAATTCCCACCAGCATGGC‑GGATCAAGTTTG‑GGGaaac                                                                                                                                             >  2:181987/1‑148 (MQ=255)
           gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGAAT‑CAGc                                                                                                                                         <  2:327315/151‑1 (MQ=255)
            cccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAAGGTGC‑ggggaaca                                                                                                                                            >  2:397941/1‑143 (MQ=255)
             ccaccaCACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGAAT‑CAGCaa                                                                                                                                       <  2:1191532/151‑1 (MQ=255)
                ccaccccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGAAT‑CAGCAAGcc                                                                                                                                    <  2:962054/151‑1 (MQ=255)
                 cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGAAT‑CAGCAAGccc                                                                                                                                   <  1:458097/151‑1 (MQ=255)
                            gaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGG‑GGGGAA‑CAGCAAAGCCG‑CCA‑Ccttt                                                                                                                        >  1:908671/1‑149 (MQ=255)
                               aCATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAGGGTGC‑GGGGAA‑CAGCAAACCCC‑CCC‑Ctttttacc                                                                                                                    >  2:1061053/1‑146 (MQ=255)
                                            cATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAGGGTGC‑GGGGAT‑CAACCAGGCCG‑CCACCTTTTACCC‑ATTTTTCGGGGCg                                                                                                      >  2:458466/1‑151 (MQ=255)
                                             atggAGGCCGCCAGCCGGGCGTGACCCGCAAAAAGCAACGGTAAAATAAAAAGGAAGGTAAATC‑CCA‑GCGGTTCGTTAATCCCCCCCAGCATGGC‑GTTTAGGGTGGC‑GGGAAT‑AAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCg                                                                                                           <  2:292813/142‑1 (MQ=255)
                                              ttGAGGCCGCCAGCACGGGGTGTACCGCAAAAAGCAAGGGAAAAATGAAAAGGAA‑GTAAATTCCCG‑GCGGTTCGTTAATCCCCACCAGCAGGGC‑GTTTAAGGTTGC‑GGGAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGt                                                                                                         <  1:853587/146‑1 (MQ=255)
                                               tGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTTC‑GGGAAT‑CAGCAAACCCC‑CCC‑CCTTTTCCCCCTTTTTCGGGGCg                                                                                                      >  2:373652/1‑148 (MQ=255)
                                               tGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTGAGGGTGC‑GGGGATAAAGCAA‑CCCG‑CCA‑CCTTTAACAGCTTTTTTGGTGCgg                                                                                                     >  2:279650/1‑149 (MQ=255)
                                                 aGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGATTAGGGTGCGGGGAAC‑CCCCAACCCCC‑CCC‑CTTTTCCCCGTTTTTCgggggggca                                                                                                    >  2:1245872/1‑143 (MQ=255)
                                                    ccgccAGCGGGGGGTTTACCGAAAACAGCAAGGGTAAAATAAACAGGAAGTAAATTT‑CCG‑GCGGTTGGATAATCCCCACCAGCAGGGG‑GTTTAGGGTGGG‑GGGAAT‑CAGCAGCCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGa                                                                                               <  2:310957/150‑1 (MQ=255)
                                                    ccgccAGCACGGCGTGTCCCGAAAACAGCAACGGTGAAATGAACAGAAAGGTAAATT‑CCG‑GCGGTTCGGTAATCCCCACCAGCATGGC‑GGTTAAGGTGGG‑GGAAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTg                                                                                                <  1:752279/149‑1 (MQ=255)
                                                     cgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAGGGTGC‑GGGGAT‑CAGCAAACCCGCCCC‑CTTTTACCC‑ATTTTTCGGGGCGGGAGTg                                                                                                >  1:1090874/1‑148 (MQ=255)
                                                       ccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTGC‑GGGGAA‑CAGCAAACCCC‑CCC‑CTTTTTCCCGATTTTTCGGGGCGGCAGGGaa                                                                                              >  1:373718/1‑148 (MQ=255)
                                                              ggCGTGTACCCAAAACAACACGGGTAAAAAAAGCAGAAAGGAAAATC‑CCG‑GCGTTTCGTTAATCCCCACCAGCAGGGC‑GTTTAGGGTTGC‑GGAAAT‑AAGCAACCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt                                                                                        <  2:339164/147‑1 (MQ=255)
                                                                  tgtACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGCAATTCCCACCAGCATGGC‑GGTTAAGATATC‑CGGAAT‑CAGCAAACCCG‑CCC‑CCCTTACCCGCGTTTTCGTTGCgccggcgaa                                                                                              >  2:1109231/1‑129 (MQ=255)
                                                                      cgccAAACAGAAACGGTAAAATGAACAGG‑AGGTAAAAT‑CCAGGGGGTTCGGTATTTCCCACGAGCATGGC‑GTTTAGGGTGGG‑GGAAAT‑CACCAGGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAt                                                                            <  1:180604/148‑1 (MQ=255)
                                                                       cGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGAAT‑CAGCAAACCCG‑CCC‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAaaaaaa                                                                             >  1:756192/1‑149 (MQ=255)
                                                                        gCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGG‑GGGAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGGGGGGCAGGGAAGTAAATCGCTAAAGaaa                                                                             >  2:525861/1‑148 (MQ=255)
                                                                          aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGC‑GGGGAT‑CAGCAAGCCCG‑CCA‑CCTTTACCCGATTTTTCGGTGCGGCAGGGAAGGACATCGCCAAAGAAATAccc                                                                        >  2:152670/1‑151 (MQ=255)
                                                                          aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGA‑GGGAAA‑CAGCAAACCCG‑CCA‑CCTTTACCCGATTTTCCGGTGCGGGAGTGAAGTTAATCGCCAaaaaaaaaccc                                                                        >  2:1239704/1‑146 (MQ=255)
                                                                            aCAGCAAGGGAAAAATAAAAAGGAAGGTAAATT‑CCG‑GGGTTGCGGTAATCCCCGCCAGCATGGC‑GTTTAGGGTTGG‑GGAAAA‑CAGAAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCAc                                                                      <  1:1250134/151‑1 (MQ=255)
                                                                              aGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGC‑GGGAAT‑CAGCAAGCCCG‑CCC‑CCTTTACCCGAATTTTCGGGTCGGGAGTGAAGGAAATCCCTAAAAAAATAccccc                                                                      >  2:577019/1‑147 (MQ=255)
                                                                                 aagggTGAAATGAACGGGAAGGTAAATC‑CCA‑GCGGTTCGGTAATCCCCCCCAGCATGGC‑GGTTAAGTTTGC‑GGGAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCAc                                                                      <  1:1168770/143‑1 (MQ=255)
                                                                                 aaCGGTAAAATGAACAGGAAGGTAAATT‑CCG‑GCGGTTCGGTAATCCCCCCCAGCATGGC‑GGTTAGGGTTGG‑GGGAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCAc                                                                 <  1:1060477/151‑1 (MQ=255)
                                                                                    ggTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGC‑GGGAAT‑CAACAAGCCCG‑CCA‑CCTTTTCCCCATTTTCCCGTGCGGCAGGGGAGGACATCGCTAAAGAAATGCCCACGGCAcc                                                                >  1:684871/1‑149 (MQ=255)
                                                                                    ggTAAAAAGAAAAGAAAGGTAAATT‑CCA‑GCGTTTGGGTAATCCCCCCCAGCATGGC‑GGTTAAGGTGGG‑GGGAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCAc                                                                 <  1:297867/148‑1 (MQ=255)
                                                                                          ataaaCAGGAAGGAAAATT‑CCA‑GCGTTCCGGTAATCCCCACCAGCATGGC‑GGTAAAGGTTGG‑GGAAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt                                                        <  1:208735/148‑1 (MQ=255)
                                                                                              aCAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGAAT‑CAGCAAACCCC‑CCA‑CTTTTACCCGATTTTTCGGGGCGGCCGTGAAGTAAATCGCTAAAAAAAAGCCCACGGCCCCCAAGAttttt                                                      >  1:1105495/1‑149 (MQ=255)
                                                                                              aCAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGC‑GGGAAT‑CAGCAAGCCCG‑CCC‑CTTTTACCCGATTTTCCGGGGGGGCAGTGAAGTAGATCGCTAAAGAAAAGCCCACGGCACCAAAGAttttt                                                      >  1:451020/1‑149 (MQ=255)
                                                                                                     gatAAATC‑CCA‑GCGGTTGGGTAATCCCCACCAGCGTGGC‑GTTTAGGGGGGG‑GGAAAT‑AAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTAcc                                              <  2:463821/148‑1 (MQ=255)
                                                                                                         aaTT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAAGATGG‑GGGGAT‑CTGCCAGCCCG‑CCA‑CATTTACCCGATTTTTCGGGGCGGCAGGGAAGGACCTCGCTCAAAAAATGCCCACGGCAACAAAGATTTTTTAGTTAcc                                              >  2:551345/1‑146 (MQ=255)
                                                                                                         aaTT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGGAG‑CCGCAAGCCCG‑CCC‑CCTTTACCCGATTTTCCGGTGCGGCGGTGAAGTACATCGCACAAAAAATACCCAAGGCACCAAAGATTTTTGAGTTACCCTg                                           >  1:1180869/1‑149 (MQ=255)
                                                                                                         aaTT‑CCA‑GCGGTTCGGTAATCCCCCCCAGCTTGGC‑GTTTAAGGTTGG‑GGGAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCa                                         <  1:340323/151‑1 (MQ=255)
                                                                                                         aaTT‑CCA‑CCGGTGGGTTAATCCCCCCGAGCATGGC‑GGTAAAGGGGGC‑GGGAAA‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCa                                         <  2:738340/151‑1 (MQ=255)
                                                                                                           tt‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTTC‑GGGGAT‑CAGCAAGCCCG‑CCA‑CCTTTTCCCGATTTTCCGGTGCGGGAGTGAAGGACAACGCTCAAAAAATGCCCACGGCACCAAAGAATTTTTAGGTTACGTGCAgg                                       >  2:684781/1‑151 (MQ=255)
                                                                                                            t‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGGGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg                                      >  2:712136/1‑151 (MQ=255)
                                                                                                            t‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGC‑GGGAAT‑CAGCAAGCCCG‑CCC‑CTTTTTCCCCCTTTTTCGGGGCGGCGGTGAAGTACATCGCTAAAAAAATACCCACCGCCCCAAAGATTTTTTAGTTACCCTGCAggg                                      >  1:592698/1‑151 (MQ=255)
                                                                                                                 accGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGGAT‑CAGCAAACCCG‑CCC‑CTTTTACCCCATTTTTCGGTGCGGCAGTGGAGTACATCGCCAAAGAAAAGCCCACGGCCCCCAAGATTTTTGAGTTACCCCGCAggg                                      >  1:354752/3‑148 (MQ=255)
                                                                                                                    ggTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGC‑GGGGAT‑CAGCAAGCCCG‑CCC‑CTTTTACCCGATTTTTCGGGGCGGCAGTGAAGGACATCGCTAAAAAAATGCCCACGGCACCCAAGATTTTTGGGTTACCGTGCAGGGCaaa                                  >  1:750269/1‑149 (MQ=255)
                                                                                                                      ttCGGTAATTCCCACCAGCATGGC‑GGTTAAAGTTGC‑GGGAAT‑CAGCAAACCCG‑CCC‑CCTTTACCCCATTTTCCGGTGCGGCCGTGAAGTACATCGCTAAAGAAATACCCACGGCACCAAAGATTTTTGAGTTACCGTGCCGGGCaaaa                                 >  2:1239143/1‑148 (MQ=255)
                                                                                                                       tCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGGAT‑CAGCAAACCCG‑CCA‑CTTTTTACCGAATTTTCGGTGGGGCAGTGGAGTACATCCGTAAAAAAATGCCCACGGGCCCCAAGATTTTTTAGTTACCCTGCAGGGCAAAAACGc                             >  1:1062945/1‑151 (MQ=255)
                                                                                                                       tCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGAAT‑CAGCAAACCCG‑CCA‑CTTTTACCCCATTTTCCGGTGCGGCAGTGAAGTACATCGCCCAAAAAATACCCACCGCACCAAAGATTTTTGAGGTACCGTGCAGGGCAAAACCGc                             >  2:684356/1‑151 (MQ=255)
                                                                                                                                ccccccAGCTTGGC‑GTTTAAGGTTGC‑GGGAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCtt                          <  1:581950/145‑1 (MQ=255)
                                                                                                                                ccccccAGCAGGGC‑GTTTAAGGTGGC‑GGAAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTc                         <  1:1245181/146‑1 (MQ=255)
                                                                                                                                ccccccAGCAGGGC‑GGTTAAGGTGGC‑GGGAAA‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCgg                      <  1:98228/149‑1 (MQ=255)
                                                                                                                                ccccccAGAATGGA‑GATTAGGGTTGG‑GGGAAG‑CAACAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTc                         <  1:330667/146‑1 (MQ=255)
                                                                                                                                 cccccaGCTTGGG‑GGTTAAGGTTGG‑GGGAAC‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaa                   <  1:488696/148‑1 (MQ=255)
                                                                                                                                           ggC‑GGTTAGGGTGGC‑GGAAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAttt          <  1:1185677/150‑1 (MQ=255)
                                                                                                                                                   aaGTTTGC‑GGAAAT‑CAGAAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCt  <  1:181886/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGC‑GGGAAT‑CAGCAAGCCCG‑CCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCT  >  NC_000913/3862513‑3862807

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: ‑