Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,668	0	A	G	40.6%	60.2 / 21.2	64	pseudogene (936/1617 nt)	glvC	pseudogene, 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 A (15/23);  new base G (26/0);  total (41/23)
Fisher's exact test for biased strand distribution p-value = 1.19e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑AATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGC  >  NC_000913/3862519‑3862813                                                                                                                                                         |                                                                                                                                                       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)                                                                                                                                                       g‑AATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGgcgc  >  2:66161/1‑147 (MQ=255)                                                                                                                                                         |                                                                                                                                                       ATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGG‑G‑AATC‑AGCAAGCCCG‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGC  >  NC_000913/3862519‑3862813

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