Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,285,853	0	G	T	42.9%	71.9 / 29.1	56	A126S (GCC→TCC)	agaD	N‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base G (5/27);  new base T (24/0);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 2.79e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGT  >  NC_000913/3285703‑3286002                                                                                                                                                         |                                                                                                                                                          gggTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTg                                                                                                                                                           <  1:607400/151‑1 (MQ=255)     tattaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTgcca                                                                                                                                                        <  1:475656/150‑1 (MQ=255)        taATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAAATTTGTCCTCATTTttcag‑c                                                                                                                                                     >  1:155877/1‑149 (MQ=255)         aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTgccag‑cc‑a                                                                                                                                                  <  1:111560/150‑1 (MQ=255)         aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTgccag‑cc‑a                                                                                                                                                  <  1:381907/150‑1 (MQ=255)                 tGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GGAATTTCTCCTTCTTTTCCAC‑CC‑C‑gggg                                                                                                                                             >  2:263097/1‑145 (MQ=255)                      gaaAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GGATTTTCTCATCATTTTCCAG‑CC‑A‑GGGGAAAct‑gc                                                                                                                                     >  2:275643/1‑148 (MQ=255)                           aGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCg                                                                                                                             <  1:328253/151‑1 (MQ=255)                              agaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGTGAATTTG‑TCATCATTTTTCAG‑CC‑AGGGGAAACCT‑CTTGGGGCCgtt                                                                                                                           >  2:432266/1‑149 (MQ=255)                              agaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCCTCATTTGCCCG‑CC‑C‑GGGGAAACC‑GCTGGGGGCggt                                                                                                                           >  2:834799/1‑148 (MQ=255)                                aaaaCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTTTTCATCTTTTTTCAG‑CC‑A‑GGGGAACCC‑GCTGGGGCCGGTTCt                                                                                                                        >  1:718982/1‑151 (MQ=255)                                 aaaCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTa                                                                                                                       <  2:200394/151‑1 (MQ=255)                                 aaaCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTa                                                                                                                       <  2:813347/151‑1 (MQ=255)                                     cgcgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTTCTCCTCATTTTCCAC‑CC‑C‑GGGGAACCT‑GCTGGGGGCGGTTCTAtttt                                                                                                                   >  1:660894/1‑151 (MQ=255)                                     cgcgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTAtttt                                                                                                                   <  2:70640/151‑1 (MQ=255)                                       cgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTc                                                                                                                 <  2:883322/151‑1 (MQ=255)                                        gACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTCACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GTGAAACCT‑GCTGGGGTCGATTCTATTTTTCg                                                                                                                <  2:259841/151‑1 (MQ=255)                                              aTTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GGATTTTCTCCTCATTTTCCAA‑CC‑A‑GGGGAACCT‑GCTGGGGGCGATTCTAtttttttcccttt                                                                                                          >  1:188993/1‑143 (MQ=255)                                              aTTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGttt                                                                                                          <  2:245247/151‑1 (MQ=255)                                               ttAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTTTTCATTATTTTTCAG‑CC‑A‑GGGGAACATGGCTGGGG‑CGGTTCTATTTTTTGCCGttt                                                                                                          >  1:614833/1‑150 (MQ=255)                                                tAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTa                                                                                                         <  2:496933/150‑1 (MQ=255)                                                       cTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTTCTCATCATTTTCCCG‑CC‑C‑GGGGAAACT‑GCTGGGGGCGGTTCTATTTTTTGCCGTTTTCCTGCtt                                                                                                  >  1:31586/1‑150 (MQ=255)                                                        tCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GGATTTGCTCTTCATTTTCCTG‑CC‑A‑GGGAAACCT‑ACTGGGGGCGGGTCTTTTTTTTTCCCCTTAACTGCtttt                                                                                                >  1:691170/1‑149 (MQ=255)                                                           aaGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTTTTTATCATTTTCCAG‑CC‑C‑GGGGAACCT‑GCTGGGGGCGATTCTTTTTTTTGCCCGTTTCCTGCTTTTCtt                                                                                             >  1:819153/1‑151 (MQ=255)                                                            aGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGG‑GAATTTG‑TCCTCCTTTTCCCG‑CC‑A‑GGGGAACCC‑GCTGGGGCCGCTTTTTTTTTTTGCCCTTTTCCTTCtttttttt                                                                                            >  1:413637/1‑151 (MQ=255)                                                             gTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCtttt                                                                                           <  2:77315/151‑1 (MQ=255)                                                                  aCTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTTCTCATCATTTTCCAG‑CC‑C‑GGGGAAACT‑GCTGGGGGCGGTTCTATTTTTTGCCGTTTTCCTGCTTTTCTTTTTTCt                                                                                       >  1:444253/1‑150 (MQ=255)                                                                     gTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTgcgc                                                                                   <  2:138683/151‑1 (MQ=255)                                                                         ggCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCg                                                                               <  2:534107/151‑1 (MQ=255)                                                                            ccAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTTCTCATCATTTTTCCG‑CC‑C‑GGGGAACCC‑TCTGGGGGCGGTTTTTTTTTTTTCCGTTTTCCTGCTTTTTTTTTTTCTGCGcggcgg                                                                              >  2:664473/1‑149 (MQ=255)                                                                                     gggATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GGATTTTTTCATCCTTTTTCAG‑CC‑A‑GGGGAAACT‑GCTGGGGGCGATTTTTTTTTTTTCCGTTTTCCTTCTTTTTTTTTTTCTTCGCGTCGGCTGGGccccc                                                                    >  1:772052/1‑148 (MQ=255)                                                                                      ggATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GGATTTTCTCCTCCTTTTTCAG‑CC‑A‑GGGGAACCT‑TCTGGGGCCGGTTCTTTTTTTTGCCGCTTTTCTGCTTTTTTTTTTTTTGCGCGGCGGCCGGGccc                                                                      >  2:428057/1‑147 (MQ=255)                                                                                       gATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTc                                                                 <  2:327865/151‑1 (MQ=255)                                                                                        aTTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAGTTTTCTCCTCCTTTTTCCG‑CCAA‑GGGAACCCG‑GCGGGGGCCG‑TTTTTTTTTTCGCCGGTTTTCTTCTTATCTTTTTTCTTCGCGTCCGGTGGGGCCAc                                                                    >  1:349654/1‑147 (MQ=255)                                                                                        aTTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCg                                                                <  2:960908/151‑1 (MQ=255)                                                                                          tGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GGATTTTCTCATCATTTTCCAC‑CC‑C‑GGGGAACCC‑GCTGGGGGCGGTTTTTTTTTTTGCCGTTTACCTGCTTTTTTTTTTTCTTCGCCGCGGCTGGGCCCCCGTCgg                                                               >  1:931629/1‑150 (MQ=255)                                                                                          tGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GGATTTGCTCATTATTTTTCAG‑CC‑A‑GGGGAACCT‑GCTGGGGGCGGTTCTATTTTTTGCCGTTTTTCTGCTTTTCTTTTTTCTGCGCGTCGGCCGGACCCCCGTCgg                                                               >  2:545747/1‑150 (MQ=255)                                                                                          tGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCgg                                                               <  1:924372/150‑1 (MQ=255)                                                                                                     ttttCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GTATTGTCTTCTTCTTTTTCAC‑CC‑A‑GGGGAAACC‑GCTGGGGGCGGTTCTTTTTTTTTCCGTTTGTCTGCCTGTTTTTTTCCTGCGCGTCGGGTGGACCCCCGCGGGTTAttttgttgg                                                   >  2:127392/1‑144 (MQ=255)                                                                                                      tttCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GGATTTTCTCTTCATTTTCCAGCCC‑G‑GGGAAACCT‑‑CTGGGGGCGCTTTTTTTTTTTTCCGTTTTCCTGCTTTTCTTTTTTCTGCGCGGCCGCTGGGCCCCCGCCGGGTAttttg                                                       >  2:172218/1‑143 (MQ=255)                                                                                                      tttCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCgg                                                   <  1:777001/150‑1 (MQ=255)                                                                                                       ttCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGc                                                  <  1:255822/150‑1 (MQ=255)                                                                                                        tCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCg                                                 <  2:801988/150‑1 (MQ=255)                                                                                                           ggTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTTCTCATCCTTTTCCAG‑CC‑A‑GGGGAAACT‑GCTGGGGGCGGTTCTATTTTTTTCCGTTTTACTTCTTATCTTTTTTCTGCGCGTCGGCTGGGCCCCCCTCGGGTTTTCAGt                                                      >  1:860980/1‑142 (MQ=255)                                                                                                           ggTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGa                                              <  2:66740/150‑1 (MQ=255)                                                                                                                aCTTTGTTGCCGATTATGGC‑GG‑GAATTTTTTCATTATTTTTCAG‑CC‑A‑GGGGAAACT‑GCTGGGGGCGGTTTTATTTTTTGCCGTTTTCCTGCTTTTCTTTTTTCTTCGCGTCGGCTGGGCCCCCCTCGGGTTTTTAGTCGGCGTGAAAGcg                                         >  1:571580/1‑150 (MQ=255)                                                                                                                 cTTTGTTGCCGATTATGGC‑GG‑GAAATTTCTCATCCTTTTCCCG‑CC‑A‑GGGGAAACT‑GCTGGGGGCGGTTTTTTTTTTTTCCGTTTTCCTTCTTTTCTTTTTTCTTCGCGGCGGCTGGGCCCCCCGCGGGTTTTTAGTCGGCGTGCAGGcg                                         >  2:131208/1‑149 (MQ=255)                                                                                                                    tgttgCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGcgatcg                                     <  1:576028/150‑1 (MQ=255)                                                                                                                    tgttgCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGcgatcg                                     <  2:380074/150‑1 (MQ=255)                                                                                                                        gCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATaa                                 <  1:649503/150‑1 (MQ=255)                                                                                                                           gATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTg                             <  1:810910/151‑1 (MQ=255)                                                                                                                              tATGGC‑GG‑GGAATTGTTCATCATTTTTCAG‑CC‑C‑GGGGAAACT‑GCTGGGGGCGGTTTTATTTTTTCCCGTTTTCCTGCTTTTTTTTTTTCTGCGCGGCGGCTGGACCCCCGTCGGGTTTTTCGTCGGCGTGAAGGCGATCGATaaa                                >  2:456843/1‑145 (MQ=255)                                                                                                                                     gg‑gAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAAc                    >  2:41524/1‑151 (MQ=255)                                                                                                                                      g‑gAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACt                   <  1:275667/151‑1 (MQ=255)                                                                                                                                         aaTTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCg                 >  1:414383/1‑151 (MQ=255)                                                                                                                                          aTTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc                <  1:849960/151‑1 (MQ=255)                                                                                                                                            ttGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCa               <  1:696454/150‑1 (MQ=255)                                                                                                                                             tGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc                >  1:614335/1‑148 (MQ=255)                                                                                                                                              gCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc                >  1:268085/1‑147 (MQ=255)                                                                                                                                                 catcatTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAgatg           <  1:549883/149‑1 (MQ=255)                                                                                                                                                    catTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGc      <  2:444214/151‑1 (MQ=255)                                                                                                                                                         gccag‑cc‑a‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGt  <  1:491813/150‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                          GGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GG‑GAATTTGCTCATCATTTGCCAG‑CC‑A‑GGGAAACCT‑GCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGT  >  NC_000913/3285703‑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: ‑