Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,285,853	0	G	T	41.0%	103.3 / 27.3	61	A126S (GCC→TCC)	agaD	N‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base G (6/30);  new base T (25/0);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 9.99e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGT  >  NC_000913/3285711‑3286002                                                                                                                                                   |                                                                                                                                                          aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTgccag‑cca                                                                                                                                                   <  1:335508/150‑1 (MQ=255)      cGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGaa                                                                                                                                             <  2:992200/150‑1 (MQ=255)          ggAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTTTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTG‑CT‑CCT‑CATTTTCCCG‑CCA‑GGGgaaac                                                                                                                                          >  1:1130609/1‑147 (MQ=255)           gaagaaAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTT‑CT‑CAT‑CATTTTTCCG‑CCA‑CGGaaa                                                                                                                                            >  1:502139/1‑146 (MQ=255)               aaaaaGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑ggg                                                                                                                                  <  2:420234/151‑1 (MQ=255)                aaaaGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTT‑TT‑CAT‑CATTTTCCCT‑GCA‑CGGGAACctgcctggg                                                                                                                                  >  1:508691/1‑146 (MQ=255)                  aaGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGcc                                                                                                                               <  2:884087/151‑1 (MQ=255)                     gagaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTT‑CT‑CAT‑CCTTTTCCAG‑CCC‑GGGGAACCCGCC‑gggg                                                                                                                                 >  1:902152/1‑146 (MQ=255)                     gagaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTT‑TT‑CAT‑CGTTTTCCGG‑CCA‑GGGgaagaagcctggggcggg                                                                                                                             >  1:271978/1‑136 (MQ=255)                       gaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTTTCTTTGCTTTTTTCTGGTTTCCTTTTTTCCCGATTTTGGC‑GGGACTTTT‑CT‑CCT‑CTTTTGCCAG‑CCA‑GGGAAACCTGCT‑TGGGCCGTTTc                                                                                                                          <  2:40245/151‑1 (MQ=255)                         aaaCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑Ta                                                                                                                       <  1:512221/151‑1 (MQ=255)                          aaCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTT‑TT‑CCT‑CATTTTCCAG‑CCA‑GGGGAAACCGCT‑GGGGGCGGTTC‑ttt                                                                                                                      >  1:621421/1‑149 (MQ=255)                            ccGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TAttt                                                                                                                    <  2:684198/151‑1 (MQ=255)                             cgcgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTT‑CT‑CAT‑TATTTTTCAC‑CCC‑GGGGAAACCCCC‑GGGGCCGGTTC‑TAtttt                                                                                                                   >  1:1144299/1‑151 (MQ=255)                             cgcgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TAtttt                                                                                                                   <  1:578475/151‑1 (MQ=255)                              gcgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTT‑TT‑CAT‑CATTTTTTAG‑CCC‑GGGGAAACTGCT‑GGGGGCGGTTT‑TCttttt                                                                                                                  >  2:671423/1‑151 (MQ=255)                               cgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTc                                                                                                                 <  1:891134/151‑1 (MQ=255)                                gACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTT‑TTCCTT‑CTTTTTCCCC‑CCG‑GGGAAA‑CTGCT‑GGGGGCGGttt‑tttttttt                                                                                                                 >  2:499439/1‑149 (MQ=255)                                gACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTT‑GT‑CATCCTTTTTCCAC‑CCG‑GGGAAACC‑GGT‑GGGGCCgggttttttttt                                                                                                                   >  1:509158/1‑138 (MQ=255)                                  cacCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGcc                                                                                                              <  1:138333/151‑1 (MQ=255)                                        tAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTG‑CT‑CAT‑CATTTGTCCG‑CCC‑GGGGAAACTGCT‑GGGGGCGGGTT‑TTTTTTTTGCCGtttt                                                                                                         >  2:728143/1‑149 (MQ=255)                                         aaaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTT‑TT‑CCT‑CCTTTTCCAG‑CCA‑CGGAAAACTGCT‑GGGGGCGGTTC‑TATTTTTTGcc                                                                                                              >  1:202828/1‑144 (MQ=255)                                           aGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTg                                                                                                     <  2:854040/151‑1 (MQ=255)                                            ggCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTGTTT‑CAA‑TTTTTTCC‑G‑CCA‑GGGGAAACAGCT‑GGGGGCGATTCTTTTTTTTCTCCGtttcctcc                                                                                                     >  2:1152793/1‑146 (MQ=255)                                             gCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTT‑TT‑TCT‑CATTTTTCAGCCCC‑GGGAAACCGCCG‑GGGGCGGTttt‑ttttttttcccttt                                                                                                           >  1:1123451/1‑137 (MQ=255)                                               cTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAATTG‑CT‑CAT‑CATTTTCCCG‑CCA‑GGGGAAACCGCT‑GGGGGCGGTTC‑TTTTTTTTTCCCTTTTCCTGCtt                                                                                                  >  1:276799/1‑150 (MQ=255)                                                     ggggCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCtttt                                                                                           <  2:712676/149‑1 (MQ=255)                                                       ggCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCttttt                                                                                          <  2:786980/150‑1 (MQ=255)                                                         cACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAATTG‑GT‑CAT‑CATTTTCCAC‑CCAGGGGAAACCT‑CT‑TGGGGCGCTTT‑TATTTTTTCCCCTTTATCCTCtttttttttt                                                                                          >  1:201385/1‑148 (MQ=255)                                                           cTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTg                                                                                      <  1:485301/150‑1 (MQ=255)                                                                     cAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTg                                                                           <  1:350255/151‑1 (MQ=255)                                                                     cAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTg                                                                           <  2:577624/151‑1 (MQ=255)                                                                      aaTCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTgg                                                                          <  1:863170/151‑1 (MQ=255)                                                                      aaTCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTgg                                                                          <  1:897276/151‑1 (MQ=255)                                                                       aTCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTTCCCC‑CCC‑GGGCAACCCCCC‑GCGGCCCCCTC‑TCTTTTTTTCCTCTTCCCTTCTTATCTTTTTTCTCCGCGGCGGCTggg                                                                         >  2:295531/1‑150 (MQ=255)                                                                       aTCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGa                                                                         <  2:47836/151‑1 (MQ=255)                                                                         cGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGAc                                                                        <  1:1153113/150‑1 (MQ=255)                                                                          gCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAATTT‑CT‑CCT‑CCTTTTCCCA‑CCC‑GGGGAACCCGCC‑GGGGGCGGTTC‑TTTTTTTTTCCGGTTTCCTGCTTTTTTTTTTTCCGCGCGGCGGCCGGGccc                                                                      >  1:91028/1‑151 (MQ=255)                                                                           ccGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTT‑TT‑CAT‑CCTTTTCCAG‑CCC‑GGGGAAACTGCT‑GGGGGCGGTTT‑TTTTTTTTGCCGGTTTCCTTCTTTTTTTTTTTCTTCCCGTCGGGTGGGcccc                                                                     >  2:1169234/1‑150 (MQ=255)                                                                             gggATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTG‑CT‑CTT‑CTTTTTTCAG‑CCC‑GGGGAAACTGGT‑GGGGGCGGGTT‑TTTTTTTTGCCGCTTTTCTGCTTTTTTTTTTTTTGCGCGTCGGCCGGGccccc                                                                    >  2:1116851/1‑148 (MQ=255)                                                                             gggATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAGATTT‑CT‑CAT‑CATTTTTCTT‑CCT‑GGGGAAACCGCT‑AGGGGCGGTTCTTTTTTTTGGCCGGTTTCCTGCTT‑TCTTTTTTCTTCCCGTggggttgatgcccc                                                                   >  1:1038992/1‑140 (MQ=255)                                                                                aTTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTTG‑TT‑CTT‑CATTTTTTCC‑CCT‑GGTGACCCTTC‑‑GGGGGCGCTTT‑TATTTTTTCTCGGTTACCTGCTTTTTTTTTTTCTGCTCGCCGGGTggg                                                                         >  2:601302/1‑141 (MQ=255)                                                                                   ggCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGt                                                              <  1:173580/150‑1 (MQ=255)                                                                                    gcgatgcgatTTTCTGGTTTACTTTGGTGCCGATTATGGC‑GGGAATTTG‑CC‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCG‑GGGGTCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGACGGTTa                                                            <  2:1071360/151‑1 (MQ=255)                                                                                     cgatgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTT‑TT‑CAT‑CCTTTTTCAG‑CCC‑GGGGAACCTCCT‑GGGGGCGGTTC‑TATTTTTTTCCGGTTTTCTGCTTATCTTTTTTCTTCGCGTCGGCTGGGCCCCCCCCGGtttt                                                           >  1:183196/1‑149 (MQ=255)                                                                                     cgatgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAGTTG‑TT‑CAT‑TACTTTTCCG‑CCA‑GGGAGAACTGCTGGGGGCGGGTTT‑TTTTTTTTGCGGTTT‑TCTGCTTTTTATTTTTCTTCGCGGCGGGTGGGCCCACgcccggttt                                                           >  2:900892/1‑143 (MQ=255)                                                                                           gatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTG‑CT‑CTT‑TCTTTTCCCC‑CCC‑CGGGAACCTGGT‑CGGGGCGATGC‑TTTTTTTTGCCGCTTATCCTCTTATCTTTTTCCTTCGCGTCGGCTGGCGCCACgcccgctctt                                                          >  2:487810/1‑137 (MQ=255)                                                                                                    gTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAg                                            <  2:1088429/151‑1 (MQ=255)                                                                                                      ttACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAgg                                           <  1:1099389/150‑1 (MQ=255)                                                                                                        aCTTTGTTGCCGATTATGGC‑GGGAATTTT‑CT‑CCT‑CCTTTTCCAG‑CCC‑GGGGAACCCGCT‑GGGGGCGGGTC‑TATTTTTTTCCGTTTTCCTTCTTTTCTTTTTTCTTCGCGGCGGCCGGGCCCACGTCGGGTTTTTTGTCGGCGGGAAgg                                           >  1:133276/1‑148 (MQ=255)                                                                                                            tgttgCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGcgatcg                                     <  2:242484/150‑1 (MQ=255)                                                                                                                  cGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGt                              <  2:898476/151‑1 (MQ=255)                                                                                                                   gATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTg                             <  2:606309/151‑1 (MQ=255)                                                                                                                    aTTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGc                            <  2:103765/151‑1 (MQ=255)                                                                                                                        tGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCg                           >  1:1151808/1‑148 (MQ=255)                                                                                                                         ggC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCgaga                        <  1:166399/150‑1 (MQ=255)                                                                                                                         ggC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCgaga                        <  2:257258/150‑1 (MQ=255)                                                                                                                           c‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCgagaga                      >  1:955821/1‑150 (MQ=255)                                                                                                                           c‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGaa                     <  2:264160/151‑1 (MQ=255)                                                                                                                               gAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTc                  >  1:610732/1‑151 (MQ=255)                                                                                                                                 aTTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc                >  2:620723/1‑151 (MQ=255)                                                                                                                                                  tGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGt  >  1:16908/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                          AATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTG‑CT‑CAT‑CATTTGCCAG‑CCA‑GGGAAACCTGCT‑GGGGCCGATTC‑TATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGT  >  NC_000913/3285711‑3286002

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