Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6610TG40.0% 38.9 / 27.8 67pseudogene (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 (15/24);  new base G (24/2);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 9.34e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑AATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGC  >  NC_000913/3862513‑3862803
                                                                                                                                                     |                                                                                                                                                     
ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTg                                                                                                                                                     <  2:65432/151‑1 (MQ=255)
      aTGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAGGGTGggg‑g                                                                                                                                                >  1:304337/1‑149 (MQ=255)
      aTGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑aa                                                                                                                                             <  1:904400/151‑1 (MQ=255)
      aTGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑aa                                                                                                                                             <  2:804427/151‑1 (MQ=255)
           gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑AATC‑AGc                                                                                                                                       <  2:223634/151‑1 (MQ=255)
                  acaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGTGG‑G‑AATC‑AGCAAGgcc                                                                                                                                 >  2:69452/1‑150 (MQ=255)
                     ccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAAGGTTCGG‑G‑GAGC‑AGAAAGAGCG‑cc                                                                                                                             >  2:85790/1‑150 (MQ=255)
                     ccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGATCGGTAATTCCCACCAGCATGGC‑GGTTTAGGGTGCGG‑G‑GATC‑AACAAGGCCT‑cc                                                                                                                             >  2:12024/1‑150 (MQ=255)
                       aaaGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAAGGTGCGG‑G‑AAAC‑AGCAAAgccg‑cccc                                                                                                                           >  2:626041/1‑148 (MQ=255)
                          gagaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAGGGTGCGG‑G‑GATC‑AGCAAACCCC‑CCA‑Ccttt                                                                                                                      >  1:776671/1‑151 (MQ=255)
                             aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGTTGGGG‑G‑GATC‑AGCAAACCCC‑CCC‑Ccttt                                                                                                                      >  1:45783/1‑148 (MQ=255)
                                   tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGGTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTGCGG‑G‑GAGG‑AACAAACCCC‑CCC‑Ctttttctccagtt                                                                                                             >  2:1037453/1‑142 (MQ=255)
                                      ggacGACATGGGGCCCCCCAGCCGGGCGTGCCCCGAAAACAGCAACGGAAAAAAAAAAAGGAAGGAAAATCCCGGGGGTTCGGTACTCCCCCCCAGC‑TGGG‑GGGTTAGGTTGGGG‑GAAATA‑AGAAAGCCCG‑CCA‑CTTTTA‑CCCGAtt                                                                                                             <  2:976772/144‑1 (MQ=255)
                                      ggTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTGCGG‑G‑GATC‑AGCAAACCCC‑CCC‑CTTTTT‑CCCGCtttt                                                                                                           >  1:863402/1‑149 (MQ=255)
                                       gTCGACATGGAGGCCGCCAGCGCGGCGTGTCCCACAAACAGCAACGGTAAAATGAACAGGAAGGTAATTCCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GTTTAGGTTGGCGG‑G‑AAAC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCg                                                                                                        <  1:1171068/151‑1 (MQ=255)
                                        tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTTCGG‑G‑GATT‑AGGAAAACCG‑CCACCTTTTA‑CCCG‑TTTTcc                                                                                                         >  1:707100/1‑149 (MQ=255)
                                         cGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAAGGTGCGG‑G‑GAAC‑AACAAACCCC‑CCC‑CTTTTT‑CCCCAtttt                                                                                                           >  1:566746/1‑146 (MQ=255)
                                          gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAAGATTGGG‑G‑GAGC‑AGCAAACCCC‑CCC‑CCTTTT‑cccccgtttttggt                                                                                                      >  2:926714/1‑139 (MQ=255)
                                           acttggggcccgccAACACGGCGTTCCCCCCAAACAACAACGTGGAAAGGAACAGGAAGGTAAACCCCAGGGTTCCGGTAATCCCCCCCAGCTTGGC‑GGTT‑AGGTTGGGG‑GAAATA‑AGAACGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGt                                                                                                      <  2:90360/140‑1 (MQ=255)
                                            cATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTGCGG‑G‑GATC‑ACCAAACCCG‑CCA‑CTTTTA‑CCCCCTTTTCCggcg                                                                                                     >  1:1058488/1‑147 (MQ=255)
                                                   gccgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTTCGG‑G‑GATC‑AGCAAACCCG‑CCC‑CTTTTT‑CCCGATTTTCCCGGGCGGGAggg                                                                                             >  1:226635/1‑148 (MQ=255)
                                                     cgccAGCAGGGGGTGTACCGCAAACAAAAACGATAAAAGGAACAGGAAGAAAAATTCCAGGGTTTCGATATTCCCCACCAGCAGGGC‑GGTTAAGGTTGCGG‑A‑ATTC‑AGAAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGaa                                                                                           <  1:1071373/150‑1 (MQ=255)
                                                         aGCACGGCGTTAACCGAAAACAGAAGCGGTGAAATAAACAGGAAGGAAAATTCCAGCGGTTGGGAAATCCCCCCCGTCAGGGC‑GTTTAAGGTGGCGG‑G‑AATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                                        <  1:1279273/149‑1 (MQ=255)
                                                         aGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTTCGG‑G‑GATC‑AGCAAACCCC‑CCC‑CTTTTT‑CCCCATTTTTCCGTGCgg                                                                                                  >  1:608603/1‑139 (MQ=255)
                                                            aCGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTGGGGG‑G‑AAAC‑A‑CAAAACCC‑CCC‑CTTTTT‑CCCCCTTTTTCGGGGGGGccggggagaac                                                                                       >  1:151502/1‑137 (MQ=255)
                                                              ggCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGATAAGGTTGTGG‑G‑AATT‑AGCCAGCCCG‑CCC‑CTTTTA‑TCCGATTTTTCCGTTCGGCAGGGAAgga                                                                                        >  1:416442/1‑142 (MQ=255)
                                                              ggCGTGTACCGAAACCAGCAAGGGTGAAATGAACGGGAAGGAAAATCCCAGGGGTTCGTTATTCCCCCCCAGCGGGGC‑GGTTAGGGTGGCGG‑G‑AATC‑AGAAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCg                                                                                   <  2:622405/149‑1 (MQ=255)
                                                               gCGTGTACCGAAAACAGCAACGGAGAAATAAACAGAAAGGTAAATTCCGGGGGTTCGGTAATCCCCACCAGCATGGC‑GGTTAGGTTGGCGG‑A‑AATC‑AGAAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTa                                                                                <  1:255709/151‑1 (MQ=255)
                                                                cGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGGGG‑G‑AAAT‑AGCAAACCCC‑CCA‑CTTTTT‑TCCGCTTTTTCGGGGCGGCGGGGGAGTACCACGCTa                                                                                >  1:529298/1‑150 (MQ=255)
                                                                 gtgacccGAAAACAGCACCGGTAAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCCCCAGTGGGGC‑GGTTAAGGTTGGGG‑G‑AACC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTaaa                                                                              <  1:384441/146‑1 (MQ=255)
                                                                    tACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGCGG‑G‑GATC‑AGCAAGACCG‑CCC‑CTTTTT‑CCCGATTTTTCGGGGCGGGAGTGGAGGAAATCGCTaaa                                                                              >  2:343285/1‑148 (MQ=255)
                                                                      ccccAAACAGCAACGGGGAAGAAAGAAGGAAGTAAAATCCCAGCGGTT‑GTTAAACCCCCCCCGGAGGGT‑CGTTTAGGGGGGGG‑G‑AAACAAGCACGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaa                                                                           <  2:1060169/146‑1 (MQ=255)
                                                                         cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTGCGG‑G‑GATC‑AGCAAACCCC‑CCC‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAATACATCGCTAAAGAAAAGcc                                                                      >  2:696331/1‑151 (MQ=255)
                                                                           accagcagGGGAGAAAGGAAAGGGAAGAAAAGTCCCAGGGGTTCGTTATTCCCCCCCAGCA‑GGG‑GGTTGGGGGTGGGGAA‑AACC‑ACCACGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGcc                                                                      <  2:331321/147‑1 (MQ=255)
                                                                           aaCAGCAACGGAAAAAGAAAGAGGAAGGAAAATTCCAGGGGTTCGTTAATCCCCACCAGCAGGGC‑GTTTAGGGTGGGGG‑A‑AATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCa                                                                    <  2:776643/151‑1 (MQ=255)
                                                                            aCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTTCGG‑G‑GATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATACCCAc                                                                   >  1:1251632/1‑151 (MQ=255)
                                                                                  aCGGTGAAATAAACGGGAAGGTAAATTCCGGCGTTTCGTTAATCCCCCCCGGCATGGC‑GGTTAGGGTTGGGG‑G‑AATC‑AGAAGGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGc                                                                <  2:137656/148‑1 (MQ=255)
                                                                                   cGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑AAAC‑AGCAAACCCG‑CCA‑CTTTTA‑ACCCATTTTCCCGGGCGGCAGGGGAGTAAAACGCTAAAAAAAAACCCACCGCA‑cc                                                            >  1:48815/1‑150 (MQ=255)
                                                                                    ggaaaaaTGAAGGGAAAGGTAACTCCCGGCGGTGCGGTAATCCCCACCAGCGGGGC‑GGTGAGGGTGGGGG‑A‑AACA‑AGCAACCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCa                                                           <  2:1209162/146‑1 (MQ=255)
                                                                                    ggTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑GATC‑AGCAAACCCG‑CCC‑CTTTTA‑CCCGATTTTCCGGGGCGGCAGGGAAGGACAACGCTTAAGAAAAGCCCACGGCA‑cc                                                            >  1:269675/1‑149 (MQ=255)
                                                                                    ggTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGTTGCGG‑G‑GATC‑AGCCAGCCCGCCCC‑CTTTTA‑CCTAATTTTCC‑GTGCGGGCGGGAAGGACATCGCaaaaaaaaatccccagg                                                                 >  2:1200109/1‑136 (MQ=255)
                                                                                        aaaggaaaaggaAGTTAATCTCCAGCGGTTCGTAAATCCCCCCCAGATTGGG‑GGTAAAGGTGGGGG‑G‑AAAA‑AGAAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGa                                                       <  2:1006857/150‑1 (MQ=255)
                                                                                         aataaaagggAAGATAAATCCCCGGGGTTCGGAAATCCCCCCCAGCGTGGC‑GGTTAGGTTTGCGG‑A‑AATC‑AGCAACCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGAtt                                                     <  2:451738/143‑1 (MQ=255)
                                                                                         aaggacagggAGGGTAAATTCAGGCGGTTGGGTAACCCCCCCCGGCATGGC‑GGTAGAGGTGGGGG‑A‑AATC‑AGCAACCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGAtt                                                     <  2:637260/143‑1 (MQ=255)
                                                                                         aaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAAGTTGGGG‑G‑GATC‑AGGAAGCCCC‑CCC‑CTTTTT‑CCCGATTTTTCGGGGGGGGCGGGGAGTACATCGCTAAAGAAAAACCCACCGCA‑CCAAAGaat                                                     >  2:420945/1‑149 (MQ=255)
                                                                                         aaTGAAAGGGAAGTTAAATCCCAGCGGTTCGGTATTCCCCCCCAGCAGGGC‑GTTAAAGGTGGCGG‑A‑AATA‑AGAAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGAtt                                                     <  1:800034/151‑1 (MQ=255)
                                                                                          aTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGATTGGG‑G‑GAGC‑AGGAAGCCCG‑CCC‑CTTTTT‑CCCCGTTTTTCGGGGCGGCAGGGGAGAAAATCGCCAAAGAAAAACCCACCGCC‑CCCAAGAatt                                                    >  2:1063286/1‑151 (MQ=255)
                                                                                           tGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTGCGG‑G‑AATC‑AGCAAGCCCG‑CCA‑CTTTTA‑TCCGATTTTTCCGTGCGGCCGAGAAGTACATCGCTCAAGAAATACCCACCGCA‑CCCAAAAtttt                                                   >  1:1151637/1‑151 (MQ=255)
                                                                                                aGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGCGG‑G‑GATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCCATTTTTCCGGGCGGCAGGGAAGAACATCGCTAAAGAAAAACCCACCGCA‑CCAAAAAttttt                                                  >  1:544760/1‑147 (MQ=255)
                                                                                                aGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGCGG‑G‑AAGC‑AGCAAACCCC‑CCC‑CTTTTT‑ACCCATTTTCCCGTGCGGCAGGGAAAGACAACGCTAAAGAAATGCCCACGGCA‑CCAAAAAttttttggt                                              >  1:62370/1‑147 (MQ=255)
                                                                                                        aaaTCCCAGGGGTTCGGTAATCCCCCCCAGCAGGGG‑GTTTAAGGTTGGGG‑G‑AATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGt                                        <  1:518150/149‑1 (MQ=255)
                                                                                                         aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAAGGTGCGG‑G‑GATC‑AGCAAGGCCC‑CCA‑CTTTTACCCCGTTTTTCCGG‑GGGGGAGTGAAGGACATCGCTTAAAAAAAACCCACCGCACCCAAAAAttttttgttaacctgc                                       >  2:659775/1‑139 (MQ=255)
                                                                                                         aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑GATC‑AGCAAGCCCG‑CCA‑CTTTTT‑CCCGATTTTCCGGTGGGGCAGTGAAGGACATCGCTAAAGAAATGCCCACCGCA‑CCAAAGATTTTTTAGTTACCGt                                        >  1:765428/1‑148 (MQ=255)
                                                                                                         aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGCGG‑G‑AATA‑AGCAAGCCCC‑CCA‑CTTTTA‑TCCCATTTTCCCGGGCGGGCGTGGAGTAAATCGCTAAAGAAATGCCCACCGCA‑CCCAAAATTTTTTAGTTACCgg                                        >  2:492924/1‑147 (MQ=255)
                                                                                                           ttCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGCGG‑G‑AATC‑AGCAAACCCG‑CCA‑CCTTTT‑CCCCATTTTCCGGTGCGGCCGGGAAGTAAATCGCTAAAGAAAAGCCCACGGCA‑CCAAAGATTTCTTAGATACCGTGCCgg                                   >  2:1034315/1‑151 (MQ=255)
                                                                                                              cAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑AATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCCATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAAGCCCACGGCA‑CCCAAGATTTTTGAGTTACCGTGCAggg                                  >  1:834886/1‑149 (MQ=255)
                                                                                                                gCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGTTGCGG‑G‑GATC‑AGCAAACCCC‑CCC‑CTTTTT‑CCCGATTTTTCGGTGCGGCAGGGAAGGACATCGCCAAAGAAAAGCCCACGGCC‑CCCAAGATTTTTTAGTTACCGTGCAggg                                  >  2:617927/1‑147 (MQ=255)
                                                                                                                 cGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTGCGG‑G‑GATC‑AGCCAGCCCC‑CCA‑CCTTTA‑CCCGCTTTTCCGGGGGGGGAGTGGAGTACAACGCTAAAGAAAAGCCCACCGCA‑CCAAAGAATTTTTAGTTTCCGTGCAgggggaa                              >  2:909522/1‑146 (MQ=255)
                                                                                                                  ggTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑AATC‑AACAAACCCG‑CCC‑CCTTTT‑CCCCATTTTTCGGGGCGGGAGGGGAGGACCTCGCTAAAAAAAAGCCCACGGCA‑CCAAAAAATTTTTAGGTAACGTGCAGGGGaaaa                             >  2:582215/1‑150 (MQ=255)
                                                                                                                    ttCGGTAATCCCCACCAGCAGGGG‑GTTTAGGGTTGGGG‑G‑ATAC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAAcc                           <  1:926752/150‑1 (MQ=255)
                                                                                                                       ggTAATCCCCCCCAG‑AAGGG‑GGTTTAGGTTGGGG‑G‑AAACAAGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCg                          <  2:812286/148‑1 (MQ=255)
                                                                                                                        gTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑AATA‑AGCAAGCCCG‑CCC‑CTTTTA‑CCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCCAAAGAAAAACCCACGGGC‑CCAAAGATTTTTTAGTTACCGTGCAGGGCAAAAACGc                         >  1:684348/1‑148 (MQ=255)
                                                                                                                           aTTCCCACCAGCATGGC‑GGTTTAGGTTGCGG‑G‑GATC‑AGCAAGCCCC‑CCA‑CTTTT‑‑CCCGATTTTCCGGTGCGGGAGTGAAGGACATCGCTTAAGAAATGCCCACCGCA‑CCCAAGATTTTTTAGGTAACGGGCAGGGCAAAAACGCCCTcc                    >  1:720082/1‑149 (MQ=255)
                                                                                                                           aTTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑AATC‑AACAAGCCCC‑CCA‑CTTTTA‑CCCGCATTTTCGGTGGGGCAGTGAAGTAAATCGCTAAAGAAATGCCCACGGCA‑CCCAAAATTTTTGTGTTACCGCGCAGGGCAAAAACGCCtttc                    >  2:709151/1‑148 (MQ=255)
                                                                                                                                  ccaGCATGGC‑GGTTAAGGGTGCGG‑G‑GATC‑AGCAAGCCCG‑CCA‑CTTTTT‑CCCGATTTTCCGGGGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCCAAGATTTTTTAGTTACCGTGCCGGGCAAAACCGCCTTCCGGGaa               >  2:1046934/1‑148 (MQ=255)
                                                                                                                                         ggC‑GGTAAGGGTGGAGG‑G‑AATG‑GGGAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTc     <  1:645759/151‑1 (MQ=255)
                                                                                                                                             ggTTAGGGTTGCGG‑G‑AATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGc  <  1:423088/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑AATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGC  >  NC_000913/3862513‑3862803

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