Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,431	0	A	C	43.9%	83.5 / 12.0	59	G34G (GGT→GGG)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (22/10);  new base C (0/25);  total (22/37)
Fisher's exact test for biased strand distribution p-value = 2.62e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGAT  >  NC_000913/2248284‑2248581                                                                                                                                                    |                                                                                                                                                       ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAcc                                                                                                                                                      <  2:411136/150‑1 (MQ=255) ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCt                                                                                                                                                     <  2:924269/151‑1 (MQ=255)  gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTg                                                                                                                                                    <  1:1029781/151‑1 (MQ=255)          aCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATg                                                                                                                                            <  1:36699/151‑1 (MQ=255)           ccGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATgg                                                                                                                                           >  1:452026/1‑151 (MQ=255)           ccGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATgg                                                                                                                                           >  1:407886/1‑151 (MQ=255)               aaaCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAAcc                                                                                                                                       >  1:664486/1‑151 (MQ=255)                aaCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCg                                                                                                                                      <  2:106294/151‑1 (MQ=255)                 aCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGc                                                                                                                                     >  1:825663/1‑151 (MQ=255)                  ccGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGcc                                                                                                                                    >  1:250462/1‑151 (MQ=255)                             gCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAatc                                                                                                                          >  2:465684/1‑150 (MQ=255)                                   aGTCCACCGAGACGGGCCATCAGGAAGTAAACCATCGGCCAAATCAACAAATACGAGAGAAATGCCCCCCCCGAACCTGCCGGAGAAATCGTGGCGCGAGAGGCAACAAAAAGCCCTGCCCAGAGGAAACCGCCAATGGCAATCATCg                                                                                                                      <  1:595261/148‑1 (MQ=255)                                    gTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTc                                                                                                                    >  2:153691/1‑149 (MQ=255)                                       cACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTg                                                                                                               >  2:52186/1‑151 (MQ=255)                                        aCCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTg                                                                                                               >  1:78130/1‑150 (MQ=255)                                         ccgggcttGGTCATCAGGAAGAAAACCAGCAGGCAAACAAGCTAATAGGAGACAAACGCCCGGCCGGCACCTGCTTGAAAATTCGTTGCGCCGGAGGCAACAAAAAAACCTTTACCGATGGAACCGCCAATGGCAATCATCGTc                                                                                                                    <  1:790818/137‑1 (MQ=255)                                          cGAGACGGGTAATCAGGAAGAAAACCATCAGCCCAATCGGCATACACGAGAGCAATGCCCCGCCCGGCCCTGCCGGAGAAATCGTTGCGCCGGAGGCAAAAAAAGGACCTGTCCGGATGAAACCGCCAATGGCAATCATCGTCAGGt                                                                                                                <  1:808022/147‑1 (MQ=255)                                            agaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgcg                                                                                                          >  1:495333/1‑151 (MQ=255)                                                  gTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                                     >  2:44140/1‑150 (MQ=255)                                                   tcatcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTa                                                                                                    >  1:115533/1‑150 (MQ=255)                                                           aaGAAAACCATCAGGCCAATCAGAATATACGAGAGCAATCCCCCGCCCGGCCCTGCCTAAGAAATCGTGGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGc                                                                                           <  2:10268/151‑1 (MQ=255)                                                               aaaCCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTa                                                                                        >  1:567730/1‑150 (MQ=255)                                                                acccaccaGGCCAATCAGCATATAGGAGAACAATCCCCCGCCCGGCCCTGCCTGAGAATTTGTTGCGCGGGAGTAAACAAAAAGGCTTTGACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTaa                                                                                       <  2:299398/148‑1 (MQ=255)                                                                            aaTCAGAATATACGAGAGCAATCCCCCGCCCGACCCTGCCTGAAAAATCTTTGCGCCGGAGGCAAAAAAAAGCCCTGCCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCt                                                                          <  1:169799/151‑1 (MQ=255)                                                                              tCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                                        >  2:608214/1‑151 (MQ=255)                                                                               cAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAAAAATCGTTGCGCCAGAGGCAACAAAAAGCCCTGTACCGATGAAACCGCCAATGCCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgtg                                                                       <  2:511892/151‑1 (MQ=255)                                                                                aGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATAGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTgg                                                                      >  1:500546/1‑151 (MQ=255)                                                                                  cATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                                    >  2:581146/1‑151 (MQ=255)                                                                                   atatACGAGAGCAATCCCCCCCCCGGCCCTGCTTAAAAAATCGTTGCGCCAGAGGCAACAAAAAGCCCTGCCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGttt                                                                   <  2:365515/151‑1 (MQ=255)                                                                                    tataCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                                                  >  2:911052/1‑151 (MQ=255)                                                                                    tataCGAGAGCAAGCCCCCGCCCGGCCCTGCCTGAAAAATCGTTGCGCAAGAGGCAAAAAAAAGCCCTGCCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                                    <  2:1027210/149‑1 (MQ=255)                                                                                     ataCGAGAGCAATCCCCCGCCCGGCCCTGCCGAAGAAATCGTTGCGCCAGAGGCAAAAAAAAGACCTGTACCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTgg                                                                      <  1:563674/146‑1 (MQ=255)                                                                                     ataCGAGAGCAATCCCCCGCCCGCACCTGCCGGAAAAATCGTGGCGCAAGAGGAAAAAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                                    <  1:11812/148‑1 (MQ=255)                                                                                     ataCGAGAGAAATCCCCCGCCCGGCCCTGCCGAAAAAATCGTTGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                                    <  2:992314/148‑1 (MQ=255)                                                                                          agagCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTc                                                            >  1:993132/1‑151 (MQ=255)                                                                                          agaccaaccccccGCCCGCCCCCGCCTGAAAAATCGTTGCCCCAGAGGCAAAAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                                                  <  2:122204/136‑1 (MQ=255)                                                                                               aaTGCCCCGCCCGGACCTGCCGGAAAAATCTTTGCCCCAGAGGCAACAAAAAGACCTGTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCg                                                           <  1:897765/147‑1 (MQ=255)                                                                                                   ccccGCCGGGCCCTGCCGGAAAAATCGTGGCGCCAGAGACAAAAAAAAGCCCGGCCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCAt                                                   <  1:502771/151‑1 (MQ=255)                                                                                                     ccgcccgGACCTGCCTGAGAAATCGTTGCGCCAGAGCCAAAAAAAAGACCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATAc                                                 <  2:1011859/151‑1 (MQ=255)                                                                                                           gaccctccctgAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAAt                                           <  2:859933/148‑1 (MQ=255)                                                                                                              cctgcctgAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCt                                        <  2:522030/151‑1 (MQ=255)                                                                                                              cctcccggAAAAATCGTGGCGCAAGAGGAAACAAAAAGCCCTGTCCGGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                                          <  2:629100/149‑1 (MQ=255)                                                                                                                   ctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATcc                                    >  1:32503/1‑150 (MQ=255)                                                                                                                      aaaaaTCGTTGCGCCAGGGGCAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                                  <  2:291503/147‑1 (MQ=255)                                                                                                                      aaaaaTCGTGGCCCAAGAGGCAAAAAAAAGCCCTGTACCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                                  <  1:21458/147‑1 (MQ=255)                                                                                                                        aaagcttggcgcAAGAGCCAAAAAAAAGCCCGTCACCTATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTAt                                      <  1:268279/135‑1 (MQ=255)                                                                                                                        aaaTCGTTGCGCCAGAGGAAAAAAAAAGCCCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                                  <  1:22316/147‑1 (MQ=255)                                                                                                                         actctctGCGCGAGAGGCAAAAAAAAGCCCAGTCCATTGGGAGCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                               <  1:1015874/143‑1 (MQ=255)                                                                                                                            cGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaaa                           >  1:257094/1‑150 (MQ=255)                                                                                                                               cccgcCAGGGGAAAAAAAAAGCCCGGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaaa                           <  1:591284/145‑1 (MQ=255)                                                                                                                                   ccAGAGGCAAAAAAAAGCCGTGTCCCGGTGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGa                       <  1:90390/147‑1 (MQ=255)                                                                                                                                     agagGCAAAAAAAAGCCCTCTCCAGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcg                   <  1:669896/149‑1 (MQ=255)                                                                                                                                      gagaaaaaaaaaaGGCCTGTACGGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaaa                           <  1:387537/132‑1 (MQ=255)                                                                                                                                       agacaaaaaAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcgc                  <  2:924493/145‑1 (MQ=255)                                                                                                                                           aacaaAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt            <  1:236774/150‑1 (MQ=255)                                                                                                                                           aaaaaaaaGCCCTGTACGGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGcgcg                 <  1:1013654/142‑1 (MQ=255)                                                                                                                                             caaAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGAGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGc         >  1:232602/1‑151 (MQ=255)                                                                                                                                               aaaaGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAg       <  2:894526/151‑1 (MQ=255)                                                                                                                                                   gACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACg    >  2:420591/1‑150 (MQ=255)                                                                                                                                                    aCCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGAt  >  2:302716/1‑151 (MQ=255)                                                                                                                                                    aCCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGAt  >  2:642335/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                       GGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGAT  >  NC_000913/2248284‑2248581

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