Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6380TC41.8% 78.2 / 12.0 55pseudogene (966/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 (20/12);  new base C (0/23);  total (20/35)
Fisher's exact test for biased strand distribution p-value = 4.81e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGG  >  NC_000913/3862488‑3862787
                                                                                                                                                          |                                                                                                                                                           
cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAAtt                                                                                                                                                             >  2:291456/1‑151 (MQ=255)
 ggTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAAtt                                                                                                                                                             <  1:279421/150‑1 (MQ=255)
        ccTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAg                                                                                                                                                    <  1:1113075/151‑1 (MQ=255)
           aGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCAt                                                                                                                                                 <  1:160850/151‑1 (MQ=255)
               aaTCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTAACGCAAACAGCAACGGTGAAA‑TGAAAAGG‑AAGGTAAATTTCAG‑CGGTTCGGTAATT‑CCCAACAGCATgg                                                                                                                                               >  2:881380/1‑149 (MQ=255)
                aTCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑gg                                                                                                                                           >  1:368078/1‑151 (MQ=255)
                         ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGtt                                                                                                                                   <  2:540913/150‑1 (MQ=255)
                         ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGtt                                                                                                                                   <  1:942274/150‑1 (MQ=255)
                         ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGTT‑g                                                                                                                                 <  2:1119206/151‑1 (MQ=255)
                          cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGTT‑Gc                                                                                                                                <  2:1134003/151‑1 (MQ=255)
                             ccATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAAGtttgggg                                                                                                                              >  1:320712/1‑146 (MQ=255)
                                      ccaccaCACCAAAGAGATACATTACGGTCGTCATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCAGCaa                                                                                                                   <  1:810976/151‑1 (MQ=255)
                                         ccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCAGCAAGcc                                                                                                                <  1:872499/151‑1 (MQ=255)
                                         ccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCAGCAAGc                                                                                                                 <  1:833205/150‑1 (MQ=255)
                                               caAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAAAAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑ACCACCAGCATGGC‑GGGTTAAGTT‑GGGGG‑GATccggcacaccgccac                                                                                                          >  2:598881/1‑137 (MQ=255)
                                                               ggccGACAGTGGGCCCGCCAGCAGGGGGT‑GAACCGAAAACAGCAAGGGAAAAA‑TGAAAGGG‑AAGGAAAACCCCGG‑CGTTTCGTTAACC‑CCCACCAGCAGGGT‑GGTTAGGGTG‑GGGGA‑AACCAGCAACACCG‑CCACTTTTA‑CCCGATTTTcc                                                                                        <  2:1198712/148‑1 (MQ=255)
                                                               ggTCGACATGGAGCCCCCCGGCACGGGGT‑GTACCGCAAACACAAATGTAGAAA‑T‑AACTGGAAATATAAATACGTG‑CGGTTCGGTAATT‑CCCCCCGGCATGGC‑‑GTTTAGGTT‑GGGGGAAACCAGAAGCCCCG‑CCACTTTTA‑CCCGAtttt                                                                                          <  2:599070/149‑1 (MQ=255)
                                                               ggTCGAAATTGAGGCCGCCACCACGGCGT‑GTCCCGCAAACAGCACGGGTAAAA‑TGAACAGG‑AAGGAAAATCCCAG‑GGTTTCGGTAATC‑CCCACCAGCATGGC‑GGTTAAGTTT‑GCGGA‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGAttt                                                                                           <  1:602206/148‑1 (MQ=255)
                                                                   gACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGGT‑GGGGG‑GATCAGCAAACCCGCCCCCTTTTA‑CCCG‑TTTTCTGGt                                                                                     >  1:1387773/1‑150 (MQ=255)
                                                                   aaatttgggccccccAGCACGGCTT‑TCACCAAAAACAGCAACGGTAAAA‑TAAACATG‑AATGTAAATTCCAG‑GGGTTCGGTACCC‑CCCCCCAGCAGGGC‑GTTAAAGGTG‑GGGGA‑AAACAGCATGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTg                                                                                    <  2:952953/138‑1 (MQ=255)
                                                                    agaTTGAGCCCCCCGGCGCGGCGT‑GCCCCAAAAACACCAGCGGAAAAA‑TAAAAAAA‑AAGATAAATCCCGG‑CGTTTC‑GTAATTCCCCCCCAGAGGGGC‑GGTTAGGGTG‑CCGGA‑AACCAGAAACCCCG‑CCACTTTTA‑CCCGAtttt                                                                                          <  2:1406955/142‑1 (MQ=255)
                                                                     cATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGCGGGTAAAGGTT‑GGGGG‑AAAAGCCAAACCCC‑CCCCTTTTC‑CCCCTTTTTCCgtg                                                                                     >  2:602477/1‑147 (MQ=255)
                                                                          aGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTTAAGGT‑GCGGG‑GAACAGCCAACCCG‑CCCCCTTTT‑CCCCCTTTTTCGGGGCgg                                                                                 >  2:373057/1‑147 (MQ=255)
                                                                           ggCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGGAATT‑CCCACCAGCATGGC‑GGGTAAGGTT‑GCGGG‑AATGAATAACACCG‑CCCCTTTTT‑CACGAGTTTTCCCTGCGGcg                                                                               >  2:856853/1‑147 (MQ=255)
                                                                            gccgccAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGGTAAGGTT‑GCGGG‑GAACAGGAAACCCC‑CCCCCTTTA‑CCCCAtttttccggggggccgtg                                                                            >  1:1093100/1‑136 (MQ=255)
                                                                             ccgccAGC‑CGGGGTGGTCCCGAAAACACCAA‑GGTGAAATTGAACAGA‑AGGGAAAATTCCGG‑GGGTCCGTTAATC‑CCCACAAGCGTGGG‑GGTTAGGGTG‑GGGGA‑ACTCGACAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGaa                                                                          <  1:494827/151‑1 (MQ=255)
                                                                             cccccggcccgccGT‑GTCCCCCAACCACAAAGGGTGAAA‑TGAGCAGG‑AAGAT‑AATCCCAGCCGGTTGTGTATCC‑CCCCCCAGGATGGG‑GGTGAAGGTG‑GGGGC‑AGCCAGAACGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCa                                                                               <  2:46914/134‑1 (MQ=255)
                                                                             cccccAGCACGGCGT‑GTCCCGAAAACAGCAAGGGTAAAA‑TGAAAAGA‑AAGGTAAATCCCAG‑GGGTGCGGTAATT‑CCCCCCAGCAGGGC‑GGTTAGGGTG‑GGGGG‑AATCAGAAGGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGaa                                                                          <  1:17556/148‑1 (MQ=255)
                                                                             ccccaaGCACGGGTT‑TTCCCACAAACAGAAACGGTGAAA‑TGAAAAGA‑AGGGTAAATTCAGG‑GGGTTCGGTAATC‑CCCCCCAGCAGGGC‑GGTAATGGTG‑GCGTG‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCa                                                                               <  1:1008380/141‑1 (MQ=255)
                                                                               gccAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTTATGTT‑TTGGG‑GGATAACCAGCCCG‑CCCCTTTTACCCCGATTTT‑CGGTGCGGGAgggg                                                                           >  2:615848/1‑145 (MQ=255)
                                                                                ccAGCACGGCGT‑TTACCGAAAACAGCAACGGTAAAA‑TGAACAGG‑AAGTTAAATTCCGG‑CGTTTCGGTAATC‑CCCCCCAGCATGGC‑GGTTAGGGTT‑GCGGA‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                       <  1:68834/151‑1 (MQ=255)
                                                                                      cGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGGTAAGGGT‑GCGGG‑GATCAGCAAACCCG‑CCCCCTTTT‑CCCGCTTTTCCGGGGCGGCAGGGGAGGACaacac                                                                 >  1:713173/1‑147 (MQ=255)
                                                                                       ggCGT‑GTCCCAAAAACAGCAACGGTAAAA‑TGAACAGG‑AAGGAAAATCCCGG‑CGGTTCGTTAATC‑CCCACCAGCTGGGG‑GTTTAGGGTT‑GGGGA‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGt                                                                        <  1:232084/143‑1 (MQ=255)
                                                                                       ggCGG‑GTCCCGCAAACAGAAGGGGTAAAG‑TGACGGGG‑AAGTTAAATGCCGG‑GGGTGCGGTAATC‑CCCCCCAGCATGGC‑GGTTAGGGTT‑GGGGG‑AATCAACACGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCg                                                                  <  1:234538/149‑1 (MQ=255)
                                                                                        ggctggtACCGAAAACAGCAAGGTTGAAA‑GAAAAAGG‑AAGTAAAATCCCGG‑CGTTTGGTAAATC‑CCCCCCAGCATGGC‑GGTAAGGGTT‑GCGGG‑AACAAACAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                       <  2:1087314/139‑1 (MQ=255)
                                                                                                ccGCAAACAGCAAGGGGAAAA‑TGAACAGG‑AAGGAAAATTCCGG‑GGGTGCGGTAATC‑CCCCCCAGCGGGGG‑GTTTAGGGTT‑GGGGA‑AATCAGAAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGa                                                           <  2:256378/148‑1 (MQ=255)
                                                                                                 cGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑GATCAGCAAACCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCCGTGAAGTACATCGCTAAAGaaa                                                         >  1:389900/1‑149 (MQ=255)
                                                                                                           aaCGGTAAAA‑TGAACAGG‑AAGGTAAATCCCAG‑CGGTTCGGTAATC‑CCCCCCAGCATGGC‑GGTTAGGGTT‑GCGGG‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACg                                                 <  1:1409147/147‑1 (MQ=255)
                                                                                                              gttaaaa‑TGAACAGG‑AAGTTAATTTCCAG‑GGGTTGGTTAATC‑CCCCCCGGCAGGGG‑GGTTAAGGTT‑GCGGA‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCa                                              <  1:624899/143‑1 (MQ=255)
                                                                                                                     aaagacgagaaaGGTAATTCCCGG‑GGTTTCGGTAATC‑CCCACCAGCAGGGC‑GTTTAGGGGG‑GGGGG‑AAGCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAtt                                     <  2:505416/141‑1 (MQ=255)
                                                                                                                            gg‑AGGGAAAATTCCAG‑CGGTTCGGTAATC‑CCCCCCAGCATGGG‑GGTTAGGGTG‑GGGGA‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGtt                             <  1:991093/151‑1 (MQ=255)
                                                                                                                               aagtaaaaTCCCGG‑CGGTTCGGTAATC‑CCCACCAGCGGGGC‑GG‑TAAGGTTGGCGGA‑AAGCAGAAACCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTa                            <  2:558935/145‑1 (MQ=255)
                                                                                                                               aaGTTAAATCCCAG‑CGGTTCGGTAATC‑CCCCCCAGCATGGC‑GTTTAAGGTT‑GGGGG‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt                                    <  1:1385983/142‑1 (MQ=255)
                                                                                                                                  gaaaaTCCCCG‑GGGTTCGGAAATC‑CCCACCGGCAGGGC‑GGTAAGGGTT‑GCGGG‑AACCACCAACCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCg                         <  1:1274139/148‑1 (MQ=255)
                                                                                                                                    aaaTTCCAG‑CGTTTCGGTAATC‑CCCCCCAGCAGGGC‑GTTTAAGGTT‑GGGGG‑AATCAGAAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCg                         <  1:1042089/148‑1 (MQ=255)
                                                                                                                                    aaaTTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑AAACAGCAAAACCG‑CCACTTTTA‑CCCGATTTTTCTGTGGGGCAGTGAAGGACATCGCCAAAGAAAAAGCCACGGCACCAAAAAATTTTTTGTTAcc                          >  2:518647/1‑147 (MQ=255)
                                                                                                                                    aaaTCCCAG‑GGTTCCGGTACTC‑CCCCCCGGAATGGG‑GGTTAGGGTT‑GGGGA‑AATAAGAAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGCGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCg                         <  2:1293892/148‑1 (MQ=255)
                                                                                                                                        tCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑AATCAGCAAACCCG‑CCCCTTTTA‑CCCGATTTTCCGGGGCGGCAGTGAAGGACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTTAGTTACCGTGCAggg                  >  1:325041/1‑151 (MQ=255)
                                                                                                                                               gtttCGGTATTC‑CCCCCCAGCTTGGC‑GTTTAGGGTT‑GCGGG‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCaa               <  1:1267900/146‑1 (MQ=255)
                                                                                                                                                  tCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCCATTTTCCGGGGCGGCAGTGAAGAACATCCCTAAAGAAATGCCCACGGGCCCAAAGATTTTTTAGTTACCGTGCAGGGCAAAACCGc         >  1:196052/1‑151 (MQ=255)
                                                                                                                                                  tCGGTAATC‑CCCCCCAGCAGGGC‑GTTTAAGGTG‑GCGGA‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCg          <  1:626587/150‑1 (MQ=255)
                                                                                                                                                   cGGTAATT‑CCCACCAGCATGGC‑GGTTTAGGGT‑GCGGG‑GATCAGCCAGCCCC‑CCCCTTTTT‑CCCCATTTTTCGGGGCGGGAGGGGAGGACATCGCTAAATAAAATCCCACCGCACCAAAGATTTTTTAGTGGCCGGGCAGGGGAAAACCGcc        >  2:1470113/1‑151 (MQ=255)
                                                                                                                                                   cGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑GATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTTCGGTGCGGCCGTGAAGGAAATCGCTAAAAAAAAGCCCACGGCACCCAAAATTTTTTAGTTACCGTGCAGGGCAAAACCGcc        >  1:870884/1‑151 (MQ=255)
                                                                                                                                                    ggTAATC‑CCCACCAGCTTGGC‑GGTTAGGGTT‑GCGGG‑AAACAGCAGGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAgg                   <  1:71447/139‑1 (MQ=255)
                                                                                                                                                      tAATT‑CCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑GATCAGCAAGCCCG‑CCCCTTTTA‑ACCGATTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTTAGTTACCGGGCAGGGCaaaaaccccct      >  1:1410854/1‑143 (MQ=255)
                                                                                                                                                      tAATT‑CCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑AATCAGCAAGCCCC‑CCCCCTTTA‑CCCGATTTTCCGGGGCGGGCGGGAAGTACATCGCTTAAGAAAAGCCCACCGCACCAAAGATTTTTTAGGTACCGTGGAGGGCAAAACCGccct      >  2:999891/1‑148 (MQ=255)
                                                                                                                                                       aaTT‑CCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGGGAAGGACATCGCTAAAGAAAAGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTcc    >  1:517393/1‑151 (MQ=255)
                                                                                                                                                         tt‑CCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AAACAGCAAGCCCG‑CCACCTTTA‑CCCGATTTTTCGGGGCGGGGGTGAAGTACATCGCTAAAGAAATGCCCACCGCACCCAAGATTTTTTAGTTACCGTGCAGGGCAAAAACGCCTTCCgg  >  2:27111/1‑151 (MQ=255)
                                                                                                                                                          |                                                                                                                                                           
CGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGT‑GTACCGCAAACAGCAACGGTGAAA‑TGAACAGG‑AAGGTAAATTCCAG‑CGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCAGCAAGCCCG‑CCACTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGG  >  NC_000913/3862488‑3862787

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