Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	42.2%	100.9 / 32.2	64	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (32/5);  new base C (0/27);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 4.76e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCC  >  NC_000913/2984294‑2984586                                                                                                                                                     |                                                                                                                                                  atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                                                                >  1:57294/1‑151 (MQ=255)   atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt                                                                                                                                               <  2:896764/150‑1 (MQ=255)        aTGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcccc                                                                                                                                         <  2:481660/151‑1 (MQ=255)          gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccc                                                                                                                                        <  2:403521/150‑1 (MQ=255)          gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccc                                                                                                                                        <  1:439004/150‑1 (MQ=255)              aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgccg                                                                                                                                    >  2:538962/1‑150 (MQ=255)                ggAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGAt                                                                                                                                  >  1:400779/1‑150 (MQ=255)                     aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACaa                                                                                                                              >  2:551151/1‑149 (MQ=255)                                 cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                >  1:633138/1‑151 (MQ=255)                                 cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                >  2:35889/1‑151 (MQ=255)                                  aaaaaaTGTATGCCGCTTTGCCACCAATCTTTATCAAAAAAATAAATGACGGCATATTTTTATGGGAAAGGAGAAAATGCA‑GATTTTTTCTTCCCTTTAAATGGTCAACCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                 <  2:568762/147‑1 (MQ=255)                                     aattgtatgcCGATTTGCCACCCATTATTTTTCAACAGATTAAGGCCGGCAGATTTTTACGGGTAGGGAGAAAATGCA‑GAGCTTTTCTCCACTTTAAGTGTAAAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                            <  1:283158/148‑1 (MQ=255)                                     aaaTGTATGTCGACTGGCCAGCCATCATTATTAAAAAAATTAAGGAGGGCAGATATTTATCGGTAGGGAGGAAATGAA‑GAGTTTTTCGTCACTTTAAATGGTAAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                            <  1:932541/151‑1 (MQ=255)                                          tgtgtCGACTTCCCAGCCATCATTATACAACAGATAAAGGACGGCAGATTTTTATCGGTAAGGAGAAAATGCA‑GAGTTTTTCGTCACTTTAAATGGTAAACGCCACCCCGGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                                                                        <  2:624660/148‑1 (MQ=255)                                           atgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                       >  2:313479/1‑150 (MQ=255)                                              tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                    >  2:78537/1‑150 (MQ=255)                                               cGACTTCCCAGCCATTATTATTCAACAGATTAATGACGGAAAATATTTAGGGGAAAGGAGAAAAGCCA‑GAGTTTTTCTTCAGGTACATGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                  <  2:881139/151‑1 (MQ=255)                                                 aCTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt                                                                                                >  1:754594/1‑151 (MQ=255)                                                     gccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                             >  1:156223/1‑150 (MQ=255)                                                      ccccccaTCATTATAAAACAGATTAAGGAGGGCAGATTTTTCTGGGAAAGGAGAAAATGCA‑GAGTTTTTCGCCACTTTAAAGGTTAAACCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                              <  1:753224/144‑1 (MQ=255)                                                        acccaTCATTATAAAAAAAATCAATGCCGATAGAAATTAATGGGTAAGG‑GGAAATTAATGAGCTATTCCACCCTTTCAATGGTAATGGCTACCCCGGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                            <  1:502911/146‑1 (MQ=255)                                                          ccaTAATTATTCAACATATTAATGCCGGAAGATATTTATGGGTAAGGAGAAAATGCA‑GACCTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                          <  2:480729/148‑1 (MQ=255)                                                               attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGGTTCACCATTTCATg                                                                                  >  1:398448/1‑151 (MQ=255)                                                               attatcCAAAAATTTAAGGAGGGCAGATTTTTACCGGTAGGGAGAAAATGGG‑GAGCTTTTCGCCGGGTTAAGGGTTCAACGCCCCCCCCGGGCCCCCGCCGATACAAGGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                                    <  1:851705/149‑1 (MQ=255)                                                                ttattCAACAGATTAATGACGGCAGATTTTTACCGGTAAGGAGAAAATGCA‑GAGTTATTCGTCACGTTAAATGGTCAATCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                                                                         <  1:190748/143‑1 (MQ=255)                                                                ttatcaaaCAAATTAATACCGGAAGATTTTTATCGGTAAGGAGAAAATGCA‑GATCTATTAGTCACTTACAATGGCAAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                                    <  2:263442/142‑1 (MQ=255)                                                                  attCAACAGATTAATGACGGCAGATATTTATCGGAAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGTTCAATGCCACACCCTGACCCCCGCCGATACAAAGCGGTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGacc                                                                               >  2:528566/1‑149 (MQ=255)                                                                   ttaaaCAGATAAATGGCGGCAGATATTTAGCGGTAAGGAGGAAATGCA‑GAGCTTTTCGTCACTTTCAATGGTAAAGGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                                <  1:362077/146‑1 (MQ=255)                                                                     aaaaaaataaaTGACGGCAAATTTTTAGGGGAAAGGAGAAAATGCA‑GAGTTTTTCGTCACGTTCAGTGGTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                             <  2:951734/141‑1 (MQ=255)                                                                           aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                      >  2:911364/1‑151 (MQ=255)                                                                           aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                      >  2:820791/1‑151 (MQ=255)                                                                             tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                     >  2:612573/1‑150 (MQ=255)                                                                                    ggCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                             >  2:323800/1‑151 (MQ=255)                                                                                      cAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                           >  1:488366/1‑151 (MQ=255)                                                                                      cAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTAGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                           >  1:305395/1‑151 (MQ=255)                                                                                       aaatttttGTCGGTAGGGAGGAAAGGAG‑GACTTATTCGTCCCGTACAATGTCAAACGCCCCCCCATGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                              <  1:857511/142‑1 (MQ=255)                                                                                        aatttttagggtaaaGGAGAAAAGGCA‑GACCTTTTCGTCCCGTAAAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt                                                               <  2:737962/132‑1 (MQ=255)                                                                                          tagttaggggaaaGGAGAAAAGGCA‑GACCTATTCGTACCGTTCAAGGGTCAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc                                                        <  2:258773/139‑1 (MQ=255)                                                                                           atTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                      >  2:871188/1‑151 (MQ=255)                                                                                             ttATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                     >  1:709034/1‑150 (MQ=255)                                                                                              tATGGGTAAGGAGGATATGCA‑GAGTTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                    <  1:554922/150‑1 (MQ=255)                                                                                                tCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                                 >  1:871106/1‑151 (MQ=255)                                                                                                  ggTGAGGAGGAAATGCA‑GTGTTATTCTTCTTTTTCAATGTTAAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                                 <  1:489202/149‑1 (MQ=255)                                                                                                     aaGGAGGAAATGCA‑GAGCTTTTCGTCACGTTAAATGGTCAACCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTa                                            <  2:257788/151‑1 (MQ=255)                                                                                                     aaGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                             >  2:141028/1‑150 (MQ=255)                                                                                                      aggaggAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                           >  2:766658/1‑151 (MQ=255)                                                                                                       ggaggaAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                          >  2:832654/1‑151 (MQ=255)                                                                                                         aagaAATGCA‑GATTTTTTAGTTACGTTAAAGGGTAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                          <  1:138489/147‑1 (MQ=255)                                                                                                           gaAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                       >  2:222694/1‑150 (MQ=255)                                                                                                            aaaTGCA‑GGGCTATTCGTCACGTTAAATGGTCAGTGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt                                      <  2:458046/150‑1 (MQ=255)                                                                                                            aaaTGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                     >  2:727099/1‑151 (MQ=255)                                                                                                               tGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                  >  1:388806/1‑151 (MQ=255)                                                                                                                gCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                 >  2:781256/1‑151 (MQ=255)                                                                                                                     agCTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                 <  1:222885/147‑1 (MQ=255)                                                                                                                            cttcgcGTTAAATGCTCAACCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                        <  2:708256/144‑1 (MQ=255)                                                                                                                            cGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                      >  2:874995/1‑151 (MQ=255)                                                                                                                             gTCCCGTTCAAGGGTCAACCCCCCACCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg                     <  1:556138/151‑1 (MQ=255)                                                                                                                              tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                    >  2:638605/1‑151 (MQ=255)                                                                                                                               cccGTTAAATGGTCAACGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                        <  1:557327/144‑1 (MQ=255)                                                                                                                                aCGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                   >  1:44168/1‑150 (MQ=255)                                                                                                                                         ggTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattga           <  1:871475/149‑1 (MQ=255)                                                                                                                                         ggTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattga           <  1:470886/149‑1 (MQ=255)                                                                                                                                                gCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAccc  >  1:292913/1‑151 (MQ=255)                                                                                                                                                 ccACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAcc   <  1:425220/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  ATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCA‑GAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCC  >  NC_000913/2984294‑2984586

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