Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,431	0	A	C	41.0%	112.8 / 22.1	78	G34G (GGT→GGG)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (38/8);  new base C (1/31);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 6.18e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGAC  >  NC_000913/2248283‑2248578                                                                                                                                                     |                                                                                                                                                     tGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAcc                                                                                                                                                    >  2:789145/1‑151 (MQ=255) tGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAcc                                                                                                                                                    >  2:460663/1‑151 (MQ=255)   gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTg                                                                                                                                                  <  1:1146547/151‑1 (MQ=255)    cAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGt                                                                                                                                                 >  1:136390/1‑151 (MQ=255)    cAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGt                                                                                                                                                 <  1:465963/151‑1 (MQ=255)        cGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCg                                                                                                                                             >  1:1371535/1‑151 (MQ=255)         gAACTGGAAAACGTGATATATGCAGCCAGTTCCTCGAGACGGGTCATCAGGAAGCAAACCCTCAGACCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTGGGGCCAGAGGCAACAAAAAGACCTGTAcc                                                                                                                                              >  2:213629/1‑149 (MQ=255)         gAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCg                                                                                                                                             >  1:157930/1‑150 (MQ=255)              ggAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑aa                                                                                                                                      <  1:530317/151‑1 (MQ=255)               gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCAGCAGGCCAATCATCATATACGAGAGCAATGCCCCGCCCGGTCCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AAc                                                                                                                                     <  2:984510/151‑1 (MQ=255)                 aaCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCg                                                                                                                                   <  2:30056/151‑1 (MQ=255)                     ggCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCaa                                                                                                                               <  1:1029954/151‑1 (MQ=255)                          atatGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAAC‑AAAAGACCTGTACCGATGG‑AACCGCCAATGGCa                                                                                                                          >  2:1112277/1‑150 (MQ=255)                               cAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAatcat                                                                                                                     >  1:418503/1‑151 (MQ=255)                                aGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAatcatc                                                                                                                    >  2:436041/1‑151 (MQ=255)                                  ccAGTTCCCCGAGACGGGTCATCAGGAAGAAACCCATCAGGCAAATCAGAATATACGAGAGCAATCCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGCAAAAAAAAAGCCCTGTCCCGAGGA‑AACCGCCAATGGCAatcatc                                                                                                                    <  1:743833/149‑1 (MQ=255)                                  ccAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGt                                                                                                                  >  1:1041317/1‑151 (MQ=255)                                    aGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGt                                                                                                                  >  2:91028/1‑149 (MQ=255)                                      tcccccGGGACGGGTCATCAGAAAAAAACCCACCAGGCAAATCAGCATAAACGGGAACAACCCCCCCCCGGGCCCTCCCTAAGAACTGGTTGCGCAAGAGGCAAAAAAAAGCCTTGAACCGAGGG‑AACCGCCAATGGCAATCATCGTCAg                                                                                                               <  2:428/146‑1 (MQ=255)                                            gagaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                                                        >  2:286728/1‑151 (MQ=255)                                              gaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCCACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACgcgc                                                                                                      >  2:331764/1‑151 (MQ=255)                                              gaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACgcgc                                                                                                      >  1:720557/1‑151 (MQ=255)                                                  ggacatgaGGAAGAAAACCACCAGCCAATTCAACATAAAGGAGAAGAAGCCCCCGCCCGCACGTGCCGAAGAACTCGTGGCGCCAGAGGCAAAAAAAAGCCCGGTCCCGAGGA‑AGCCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                                  <  2:1036620/148‑1 (MQ=255)                                                    tcatcaGGAATTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTa                                                                                                 >  2:409294/1‑150 (MQ=255)                                                          ggAAGAAAACCACCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGCCCTGCCTGAGAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTaa                                                                                                <  1:1347142/145‑1 (MQ=255)                                                               tAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGt                                                                                      >  2:1267461/1‑150 (MQ=255)                                                                aacccaccaGCCAAATCAGCATAAACGAGAGAAACCCCCCGCCCGCCCCTGCCGGAAAAATCGTGGCGCCAGAGAAAAAAAAAAGCCCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTa                                                                                     <  2:1297070/147‑1 (MQ=255)                                                                aaaCCATCAGCCCAACCAGCATAAACGAGAGCAATCCCCCGCCGGGCCCGGCCGGAGAAATCTTCGCGCCAGAGGCAAAAAAAAGCCCGGTACGGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTaa                                                                                    <  1:592729/151‑1 (MQ=255)                                                                aaaCCACCAGGCCAAGCAGAATATACGAAACAAATGCCCCGCCCGGCCCTGCCTGAGAAATCGTTGCGCCAGGGGCAAAAAAAAGCCCTTTACCGAGGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTaa                                                                                    <  2:1185674/151‑1 (MQ=255)                                                                        aGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGCAATCGTTGCGCCAGAGGCAACAAACAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGCCAGGTGACGCGCCTTTACTTCACGGCGTAAGCCCgg                                                                              >  2:729046/1‑149 (MQ=255)                                                                        aGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgcg                                                                            >  1:444708/1‑151 (MQ=255)                                                                            cAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTc                                                                        >  2:104122/1‑151 (MQ=255)                                                                            cAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTc                                                                        >  1:397818/1‑151 (MQ=255)                                                                               tCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                                     >  2:1098635/1‑151 (MQ=255)                                                                                 aGCATATACGAGAGCAAGGCCCCGCCCGGACCGGCCTGAGAAATCGTTGCGCAAGAGGAAAAAAAAAGCCCTGTCCCGATGA‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                                     <  1:345072/149‑1 (MQ=255)                                                                                   cATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                                 >  2:109584/1‑151 (MQ=255)                                                                                   cATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                                 >  2:135662/1‑151 (MQ=255)                                                                                   cATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGACGTAAGCCCGGCGCTTCTGTGGtt                                                                 >  1:392749/1‑151 (MQ=255)                                                                                    atatACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGttt                                                                >  1:970541/1‑151 (MQ=255)                                                                                     tataCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                                               >  2:418470/1‑151 (MQ=255)                                                                                      ataCGAGAGCACTCCCCCGCCCGGCCCGGCCTGAGAAATCGTTGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGATGG‑AGCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                                 <  2:606692/148‑1 (MQ=255)                                                                                       taCGAGAGCAATGCCCTGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                                               >  1:1219499/1‑149 (MQ=255)                                                                                        aCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGt                                                            >  2:395865/1‑151 (MQ=255)                                                                                         cGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGtt                                                           >  2:1331431/1‑151 (MQ=255)                                                                                            gggCAATGCCCCCCCCAGCCCGGCCGGAGAAATGGTTGCGACAGGGGAAAAAAAAAGCCCTGTGCGGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTa                                                              <  2:971366/143‑1 (MQ=255)                                                                                             accAAGGCCCCGCCGGGACCTGCGGGAGAAATCTTGGCGCCAGAGGCAAAAAAAAGCCCTGTACGGAGGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTgg                                                                   <  2:223999/137‑1 (MQ=255)                                                                                              gCAATGCCCCGCCCGACCGTGCCTGAAAAACCTTTGCGCCAGAGGAAAAAAAAAGCCTGGCCCCGAGGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTc                                                         <  1:427814/148‑1 (MQ=255)                                                                                                aaTGCCCCGCCCGGCCCTCCCTGAAAAATCGTGGCGCAAGAGGAAAAAAAAAGCCCTGTCCCGATGA‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCAGCGCTTCTGTGGTTTTAGTTTCGGaa                                                     <  1:789710/150‑1 (MQ=255)                                                                                                aaTCCCCCGCCGGGACCGCCCTGAAAAATCGTTGCGCAAGAGACAAAAAAAAGCCCTGTCCGGAGGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaaa                                                    <  1:692839/151‑1 (MQ=255)                                                                                                aaTCCCCCCCCGGGACCTGCATGAAAACTCGTTGCGCCAGAGCAAAAAAAAAGCCCTGCACGGATGA‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaaa                                                    <  1:1123178/151‑1 (MQ=255)                                                                                                aaTCCCCCCCCCGCCCCTGCCAGAGAACTTTTTGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGATGA‑AGCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGttt                                                          <  2:1327565/145‑1 (MQ=255)                                                                                                 aTGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTa                                                              >  1:1045633/1‑140 (MQ=255)                                                                                                    ccccccccGCCCCTGCCTGAAAAATCGTTGCGCCAGAGCAAAAAAAAAGCCCTGTCCCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAc                                                   <  1:1089737/148‑1 (MQ=255)                                                                                                    ccccGCCCGGCCCTGCCTGAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCa                                                 <  1:774036/150‑1 (MQ=255)                                                                                                    ccccGCCCGGACCTGCCGGAAAAATCTTTGCGCCAGAGGCAACAAAAAGCCGTGCACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGttt                                                          <  2:421908/141‑1 (MQ=255)                                                                                                    ccccGCCCGCCCCTGCCGGAGAAATCTTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGAGGA‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTa                                                              <  2:1044976/137‑1 (MQ=255)                                                                                                        gcccgGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACa                                             >  2:762728/1‑150 (MQ=255)                                                                                                         cccgGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGCCCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGAGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaa                                            >  2:496703/1‑150 (MQ=255)                                                                                                           cgGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaaaa                                          >  2:976073/1‑150 (MQ=255)                                                                                                           cgGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaaaa                                          >  2:49899/1‑150 (MQ=255)                                                                                                             accctgccggAAAATCCGTGGCCCCAGAGGCAAAAAAAAGCCCTGTCCCGAGGA‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAAt                                        <  2:921731/148‑1 (MQ=255)                                                                                                               cctgccagAGAAATGGTGCCCCCAGAGCAAAAAAAAAGCCCTGCCCCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGc                                      <  2:1000336/150‑1 (MQ=255)                                                                                                               ccggcctgAGAATTCGTGCCGCAAGAGCAAAAAAAAAGCCCTGGCCAGAGGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                                       <  1:24293/146‑1 (MQ=255)                                                                                                                   cctgAGAAAGCGTTGCGCCGGAGGCAAAAAAAGGCCCTGTCCCGATGA‑ACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTATGTGGTTTTAGTTTCGGAAACCATACAAAAATg                                       <  2:170073/145‑1 (MQ=255)                                                                                                                     tgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGAGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                               >  2:228710/1‑151 (MQ=255)                                                                                                                       agaAATTGTGGCGCAAGAGGCAAAAAAAAGCCCTGTCCCAATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCa                                <  2:649496/148‑1 (MQ=255)                                                                                                                       aaaaaTCGTTGCCCCAGAGGAAACAAAAAGCCCTGTACCGTGGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                               <  1:585619/147‑1 (MQ=255)                                                                                                                         aaggcgtgccgcCAGAGGCAAAAAAAAGCCCTTTCCCGAGGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                            <  2:1370777/141‑1 (MQ=255)                                                                                                                               ttGCGCCAGAGGCAAAAAAAAGCCCTTTCCCGAGGA‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaa                         <  1:995517/147‑1 (MQ=255)                                                                                                                               ttGCCCAAGAGGCAAAAAAACGACCTGTACCGAGGA‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATcc                                 <  1:601195/139‑1 (MQ=255)                                                                                                                                     cAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcgc               >  2:405083/1‑151 (MQ=255)                                                                                                                                        agcaaaaaaaaaGCCATGTACGGACGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATg                     <  2:1016146/135‑1 (MQ=255)                                                                                                                                            aaaaaaaaGCCCTGT‑CCGAAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTa        <  2:531794/148‑1 (MQ=255)                                                                                                                                             acaaAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTAAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCCTCTGTGGTTTTAGTTTCGGAAAACATACAAAAATGCTATCCATATTAAAAATGAGGAGCGAtt           >  2:1001164/1‑147 (MQ=255)                                                                                                                                                aaaaGCCTTGTCCCGAGGG‑ACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAg    <  2:56045/151‑1 (MQ=255)                                                                                                                                                aaaaGACCTGCACCGAGGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAg    <  2:432433/151‑1 (MQ=255)                                                                                                                                                 aaaGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGc      >  2:631966/1‑148 (MQ=255)                                                                                                                                                  aaGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGAc  >  2:1296889/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGG‑AACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGAC  >  NC_000913/2248283‑2248578

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