Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,431	0	A	C	46.4%	90.5 / 27.4	70	G34G (GGT→GGG)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (29/8);  new base C (0/32);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 1.68e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGA  >  NC_000913/2248282‑2248577                                                                                                                                                        |                                                                                                                                                    gTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGAc                                                                                                                                                    >  1:1019542/1‑151 (MQ=255)   ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCt                                                                                                                                                  >  1:967556/1‑151 (MQ=255)    gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑g                                                                                                                                                <  2:238417/151‑1 (MQ=255)    gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑g                                                                                                                                                >  1:952085/1‑151 (MQ=255)        aCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTAcc                                                                                                                                            <  2:91314/151‑1 (MQ=255)          gAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCg                                                                                                                                           >  1:793695/1‑150 (MQ=255)           aaCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGAt                                                                                                                                         >  1:914120/1‑151 (MQ=255)            aCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATg                                                                                                                                        >  1:388978/1‑151 (MQ=255)               ggAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGaa                                                                                                                                     <  2:1035599/151‑1 (MQ=255)                gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAAc                                                                                                                                    <  2:384192/151‑1 (MQ=255)                   aCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGc                                                                                                                                 <  2:125015/151‑1 (MQ=255)                     cGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCa                                                                                                                               <  2:742700/151‑1 (MQ=255)                      ggCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCa                                                                                                                               >  1:890113/1‑150 (MQ=255)                      ggCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCa                                                                                                                               <  2:689726/150‑1 (MQ=255)                         atatatGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCaa                                                                                                                              >  2:278807/1‑148 (MQ=255)                                cAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAatcat                                                                                                                    >  2:367310/1‑151 (MQ=255)                                         ccccgggccTGATCATAAGGAAGAAAACCAACAGCCAAATCA‑GCATATACGAGAGAAATGCCCCCCCGGGCCCTGCCGGAAAAATCGTTGCGCAA‑GGGGCAAAAAAAAACCCT‑GTCCGGATGAAACCGCCAATGGCAATCATCGTCAgg                                                                                                             <  1:1055526/141‑1 (MQ=255)                                            cGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAc                                                                                                         >  2:578776/1‑150 (MQ=255)                                             gagaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                                                       >  1:1012164/1‑151 (MQ=255)                                                   ggacatcaGGAAAAAAACCATCA‑GCCAAACAGGCAAAAACGAGAGCAACGCCCCGCCCGGACGGGCCAGAGAAGGCGTGCCGCAA‑GAGACAAAAAAAAGCCCT‑GTCCGGAGGGAACCGCCAATGGCAATCATCGTCAg                                                                                                              <  2:831384/135‑1 (MQ=255)                                                   ggTCATCAGGAAGAAAACCATCGGCCCAATCA‑GCATAAACGAGAGCAATCCCCCCCCCGGCCCTGCCTGAAAAACCGTTGCGCCA‑GAGGCAAAAAAAAGCCCT‑GTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                                 <  1:719531/151‑1 (MQ=255)                                                          agaaaGTAACCCACCAGGCAAACCA‑GCATATACGAGAGCAACCCCCCCCCCGACCCTCCCTGAAAAATCG‑TGCCCCAGGAGGAAAAAAAAACCCCG‑TTACCAATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCa                                                                                           <  2:908930/147‑1 (MQ=255)                                                           ggAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACg                                                                                         >  2:893712/1‑151 (MQ=255)                                                             aaGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGACCCGCACGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCCATCATCGTCAGGTGACGCGCCTTTAATTCACgg                                                                                        >  2:402383/1‑150 (MQ=255)                                                             aaGAAACCCATCAGCCAAATCA‑GCAAAAACGAGACCAAGCCCCCCCCCGGACCTCCCTGAAAAATCGTGGCCCCA‑GAGGAAAAAAAAAGCCCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGGGACGCGCCTTTAATTCACgg                                                                                        <  2:261538/150‑1 (MQ=255)                                                                aaaaCCACCGGGCCAACCA‑GCATATCCGAGAGAAATCCCCCGCCCGGCCCTGCCTGAGAAATCGTTGCGCCA‑GAGGAAAAAAAAAGCCCT‑GTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTa                                                                                    <  1:333682/150‑1 (MQ=255)                                                                 aaaCCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGAGACGCGCCTTTAATTCACGGCGTa                                                                                    >  2:1060857/1‑150 (MQ=255)                                                                 aaaCCATCAGGCCAATCA‑GCATATACGAGAAAAAGGCCCCGCCCGGCCCTGCCTGAAAAATCGTTGCGCCA‑GGGGCAAAAAAAAGCCCT‑GTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGt                                                                                     <  1:283309/149‑1 (MQ=255)                                                                      agcAGCCCAACCA‑GCAAATACGAGAGAAATGCCCCGCCCGGGCCTGCCGGAGAAAGCGTGGCGCCA‑GAGGAAAAAAAAAGCCCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGccc                                                                               <  2:222087/148‑1 (MQ=255)                                                                         agcccAAGCA‑GCAAATAGGAGACCAAGCCCCCCCCCGGCCCTCCCAGAAAACTCGTGGCGCAA‑GAGCAAAAAAAAAGCCCT‑GTCCTGATGGATCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgc                                                                            <  2:117538/147‑1 (MQ=255)                                                                             cAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTc                                                                       >  2:516770/1‑151 (MQ=255)                                                                              aaccA‑GCATATACGAGAGAAACCCCCCCCCCGGCCCTGCTTGAGAAATCGTTGCGCCA‑GAGGAAAAAAAAAGACCT‑GTCCCAAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgc                                                                            <  2:294445/142‑1 (MQ=255)                                                                                       tataGGAGAGCAATCCCCCCCCCGGCCCTCCCTAAAAAATCGTTGCGCAA‑GAGGAAAAAAAAAGCCCT‑GAACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGttt                                                               <  1:362600/150‑1 (MQ=255)                                                                                        ataCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                                              >  1:416115/1‑150 (MQ=255)                                                                                            gagaGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGttt                                                         >  2:167609/1‑151 (MQ=255)                                                                                            gagaGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGttt                                                         >  1:592779/1‑151 (MQ=255)                                                                                             agagaaaccccccccCGGGCCCTGCCTGAAAAATCGTTGCGCCA‑GAGGCAAAAAAAAGCCCT‑GTCCCAATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAg                                                            <  1:862544/138‑1 (MQ=255)                                                                                             agagCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGACGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTc                                                        >  2:864825/1‑151 (MQ=255)                                                                                              gagaAATGCCCCGCCCGGCCCTGCCGGAAAATTCTTTGCGCAG‑GAGGCAAAAAAAAGCCTT‑GTCCGGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTc                                                        <  1:383739/150‑1 (MQ=255)                                                                                                  attcccccccCGGGCCCTGCCTGAGAAATCGTTCCGCCA‑GAGGCAAAAAAAAGCCCT‑GTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCgg                                                      <  1:944973/144‑1 (MQ=255)                                                                                                  aaTGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                                    >  2:1056485/1‑150 (MQ=255)                                                                                                  aaTGCCCCCCCCGGCCCTGCCTGAAAAATCGTGGCGCGA‑GAGGCAAAAAAAGGCCCT‑GTACCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTc                                                        <  1:734991/146‑1 (MQ=255)                                                                                                  aaTCCCCCGCCCGGCCCTGCCTGAAAAATCGTTGCGCCA‑GAGGCAAAAAAAAGCCCT‑GTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGttt                                                               <  1:193567/139‑1 (MQ=255)                                                                                                     cccccGCCCCCACCTCCAGGAAAACCCGTGGCCCAA‑GGGGCAAAAAAAACCCCG‑GTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAcc                                                 <  2:668314/149‑1 (MQ=255)                                                                                                      ccccGCCGGGCCCGCCCTGAGAACTCGTTGCCCCA‑GAGGCAAAAAAAAGCCCT‑GTCCGGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                                    <  2:1042523/146‑1 (MQ=255)                                                                                                       cccgcccgGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCAt                                               >  1:450888/1‑150 (MQ=255)                                                                                                         cgcccgGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAACGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATAc                                             >  2:674096/1‑150 (MQ=255)                                                                                                           cccgGCCCTGCCTGAAAAATCGTTGCGCCA‑GAGGAAAAAAAAAGCCCT‑GTACCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACa                                            <  1:849803/149‑1 (MQ=255)                                                                                                           cccgGACCTGCCTGAGAAATCGTTGCGCCA‑GAGCAAAAAAAAAACCCT‑GTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaaa                                          <  1:572702/151‑1 (MQ=255)                                                                                                           cccgGACCTGCCTGAGAAATAGTTGCGCCA‑GAGGCAACAACAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaa                                           >  2:302785/1‑150 (MQ=255)                                                                                                               gACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                                      >  1:57123/1‑151 (MQ=255)                                                                                                                aCCTGCTGGAGAATTCGTTGCGCGA‑GAGGCAA‑AAAAAGGCCTGGAACCAAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                                      <  1:991051/150‑1 (MQ=255)                                                                                                                aCCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGCTGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGc                                     >  2:1061627/1‑151 (MQ=255)                                                                                                                     ccggaaaaaTGGTGGCGCAA‑GAGCAAAAAAAAAGCCCT‑GTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATcc                                <  1:665523/145‑1 (MQ=255)                                                                                                                         aaaaaTCGTTGCGCAA‑GAGGCAAAAAAAAGCCCT‑GTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATc                             <  1:927209/148‑1 (MQ=255)                                                                                                                          gaAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                           >  2:943811/1‑151 (MQ=255)                                                                                                                           aaaTCGTTCCCCAA‑GAGCAAAAAAAAAGCCTT‑GTCCCGATGAAGCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTa                          <  1:936346/151‑1 (MQ=255)                                                                                                                            aaTCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTa                          >  2:763227/1‑150 (MQ=255)                                                                                                                                 ttGCGCCA‑GAGGCAAAAAAAAGCCCT‑TTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAAt                     <  2:970651/150‑1 (MQ=255)                                                                                                                                    cgcCA‑GAGGCAAAAAAAAGCCCT‑GTACCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATg                    <  1:525984/148‑1 (MQ=255)                                                                                                                                      ccA‑GAGCCAACAAAAACCCCT‑CTGCCAAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATg                    <  1:278986/146‑1 (MQ=255)                                                                                                                                          gagGCAAAAAAAAGCCCT‑GTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgc                <  1:68461/147‑1 (MQ=255)                                                                                                                                               aaaaaaaaGCCCT‑GTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt        <  1:790647/147‑1 (MQ=255)                                                                                                                                                acaaAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt        >  2:725931/1‑149 (MQ=255)                                                                                                                                                  aaaaaGCCTT‑GTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTg         <  1:48741/146‑1 (MQ=255)                                                                                                                                                  aaaaaGCCGT‑TTCCCGATAAACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCa    <  2:793399/151‑1 (MQ=255)                                                                                                                                                  aaaaaGCCCT‑GTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAg                  <  2:143735/137‑1 (MQ=255)                                                                                                                                                  aaaaaGCCCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCa    <  1:727996/151‑1 (MQ=255)                                                                                                                                                  aaaaaGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGc     >  1:898625/1‑150 (MQ=255)                                                                                                                                                    aaaGCCCT‑GTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGa  <  1:450209/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                    GTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCA‑GCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCA‑GAGGCAACAAAAAGACCT‑GTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGA  >  NC_000913/2248282‑2248577

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