Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,511,772	0	A	C	50.8%	74.4 / 99.5	59	V312G (GTG→GGG)	mntH	manganese/divalent cation transporter
Reads supporting (aligned to +/- strand):  ref base A (26/3);  new base C (0/30);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 9.23e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.60e-09

GAGCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGAAA  >  NC_000913/2511623‑2511919                                                                                                                                                      |                                                                                                                                                    gagCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATcaccac                                                                                                                                                    >  1:1061874/1‑151 (MQ=255)  agCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCaccacc                                                                                                                                                   <  2:1211392/151‑1 (MQ=255)  agCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCaccacc                                                                                                                                                   <  2:287281/151‑1 (MQ=255)   gCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCt                                                                                                                                                  >  1:310686/1‑151 (MQ=255)    cGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTg                                                                                                                                                 >  2:314189/1‑151 (MQ=255)         ccAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCg                                                                                                                                            >  1:422439/1‑151 (MQ=255)            aaaCTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGcc                                                                                                                                          >  1:396160/1‑150 (MQ=255)              aCTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAgt                                                                                                                                       >  2:144074/1‑151 (MQ=255)                           tGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCg                                                                                                                          >  1:842691/1‑151 (MQ=255)                                 aTAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGa                                                                                                                      >  1:1052100/1‑149 (MQ=255)                                   acccAGAAGCCGGGTGGGCTCTAATCCCTTTAGACGATCAAAAAAGGACGGCAACTTGGTGATTGACCGACCCACCCCCAGCGGCATTGGAAAGGGAATGAACCCCTGCATCACCACCGGCCCCGCCAGTGTCCCCACCACCTTTGAAGa                                                                                                                   <  1:875840/148‑1 (MQ=255)                                   aaCCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGAc                                                                                                                  >  2:823452/1‑151 (MQ=255)                                    aCCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGAc                                                                                                                  >  1:1053551/1‑150 (MQ=255)                                       aGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGcc                                                                                                              >  1:1165550/1‑151 (MQ=255)                                            ccGTGTCGGATCTAATCCCAGAAAAATAACAATAAATGAGGGCAACAGGGGGACTGTGCGACGCCCCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCg                                                                                                           <  2:697276/149‑1 (MQ=255)                                              gtgtCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCgcagc                                                                                                       >  2:609478/1‑151 (MQ=255)                                                gtCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCgcagca                                                                                                      >  2:956476/1‑150 (MQ=255)                                                   ggATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAAcc                                                                                                   >  1:1095964/1‑150 (MQ=255)                                                   ggATCTAACCCCATCAGAATGACAAAAAATACCGGCAACAGGGTGACTGTACGCCCCACCCACAGCGGGATAAAAAAGCAAATGAAGCCCTGCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCa                                                                                                  <  1:618927/151‑1 (MQ=255)                                                   ggATCTAACCCCATCAGAATGAAAAAAAAGGACGGCAACATGGTGACTGTACGACGCCCCCACAGCGGGATATAAAAGCGAATGAAGCCCTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCa                                                                                                  <  1:99030/151‑1 (MQ=255)                                                   ggATCTAACCCCACCAGAATGAAAAAAAATGCCGCCAACAGGGTGCCTGCACACCCCCCCCACAGCGGGATATAAAAGCAAATGAACCCCTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAAcc                                                                                                   <  1:804180/150‑1 (MQ=255)                                                         aaccccacaaGAATGACAATAAAGGCCGGCAACAGGGTGCCTGTCCGCCCCCCCCACAGGGGAATATGAAACCGAATGAAGCCCTGCATCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAg                                                                                                 <  1:1215018/137‑1 (MQ=255)                                                         aaccccATAAAAATAAAAAAAAATGACGGAAACATGGTGACTGTACGCCGCCCCCCCAGCGGGATATAAAAGCGAATGAACCCCTGCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAg                                                                                                 <  2:1035548/143‑1 (MQ=255)                                                           tCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAAt                                                                                          >  2:188606/1‑151 (MQ=255)                                                                   gAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGAc                                                                                  >  2:1150031/1‑151 (MQ=255)                                                                    aaTAACAATAAATAACGGCACCATGGTGGATGTACGACGCCCCCACAGCGGAATATGAAAGCGAAAGAGCCCCTGCCTCCCCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGAcc                                                                                 <  2:683031/151‑1 (MQ=255)                                                                     acgACAAAAAATTACGGAAACAAGGGGACGGTACGACGCCCCCACAGCGGTATAGGAAAGAAAAAGAAGCCCTACAAACCCCCCTACCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGAcc                                                                                 <  2:722834/148‑1 (MQ=255)                                                                     aaaaaaaTAAATGACGCCAACATGGTGCCTGTCCGACGCCCCCACGGCGGAATATGAAAGAGAATAAAGCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCg                                                                                <  2:303407/146‑1 (MQ=255)                                                                        aaaaaaaaggccGGCAACAGGGTGACTGTCCGCCCCACCCACAGCGGGATATAAAAGGGAATGAACCCCTGCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGt                                                                               <  2:702486/138‑1 (MQ=255)                                                                        aCAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTg                                                                             >  1:942360/1‑151 (MQ=255)                                                                          aaTAAATGACGGAAACATGGTGACTGTACGACGCCCCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGtt                                                                              <  2:36048/148‑1 (MQ=255)                                                                           aaaaaTGACGGCAACATGGGGACTGTCCGACGCCCCCACAGCGGGAAATAAAAGCGAAAGAATCCCTGCACCACCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCg                                                                          <  2:101689/149‑1 (MQ=255)                                                                                tGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCg                                                                     >  1:639912/1‑151 (MQ=255)                                                                                  aCGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTg                                                                   >  2:879064/1‑151 (MQ=255)                                                                                  aCGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTg                                                                   >  2:393869/1‑151 (MQ=255)                                                                                       aaCATGGTGCCTGTACGCCCCCCCCACAGCGGAATATGAAAGCGAATGAAGCCCGCCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCt                                                                <  1:618957/149‑1 (MQ=255)                                                                                       aaCATGGTGCCTGTACGACGCCCCCACAGCGGGATAGGAAAGCGAATGAACCCCTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTgg                                                                  <  1:1163406/147‑1 (MQ=255)                                                                                       aaCATGGTGATTGTACGGCCCCCCCACAGGGGGATAGGAAAGCAAATGAAGCCCTGCATCCCCCCCCGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGc                                                                      <  2:765293/143‑1 (MQ=255)                                                                                       aaCATGGTGACTGTACGCCCCCCCCACAGCGGGATATGAAAGCGAATGACCCCCTCCATCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCtt                                                               <  2:360716/150‑1 (MQ=255)                                                                                        aCAGGGTGCCTGTACGCCCCCCCCACAGGGGGAGAGGAAAGAGAATGAAGCCCCGCCCCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTaa                                                             <  1:190195/151‑1 (MQ=255)                                                                                            ggTGACTGTACGACGCCCCCACAGGGGGATATGAAAGCGAATGAACCCCTGAATCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAg                                                          <  2:1165679/150‑1 (MQ=255)                                                                                                   gtccgaccCCCCCACGGGGGGATAGAAAAGCGAATGAACCCCTGCACCCCCCCCCCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGc                                                   <  1:1055686/147‑1 (MQ=255)                                                                                                      caacgCCCCCACAGCGGGATATGAAAGCGAAGAAACCCCTGCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTTCAGc                                                <  1:641674/148‑1 (MQ=255)                                                                                                       gacccccccACGGCGGGTTTAGAAAGCGAATAAAGCCGTGCACCCCCCCCGCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGc                                                <  1:290498/143‑1 (MQ=255)                                                                                                           ccccccCAGCGGGTTATAAAAGCAAATGAACCCCTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAg                                           <  1:792316/148‑1 (MQ=255)                                                                                                           cACCAACAGGGGGATAAGAAAGCGAATGAAGCCCGGCATCACCACCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGt                                              <  2:471712/147‑1 (MQ=255)                                                                                                             cccACAGGGGAATATAAAAGCAAAGGAAGCCCTGCATCACCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAg                                           <  2:175794/148‑1 (MQ=255)                                                                                                             cccACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTc                                             <  2:1121524/146‑1 (MQ=255)                                                                                                               cacaGCGGGATATGAAAGCGAATGAACCCCTGCACCACCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAg                                      <  1:355686/151‑1 (MQ=255)                                                                                                                  aGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCt                                   >  1:1218115/1‑151 (MQ=255)                                                                                                                   gCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCtc                                  >  2:457548/1‑151 (MQ=255)                                                                                                                    cGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCtc                                  >  1:937499/1‑150 (MQ=255)                                                                                                                             aaaGCGAATGAACCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGa                           <  2:437290/148‑1 (MQ=255)                                                                                                                              aggCGAATGAACCCCTCCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCgg                       <  1:908481/149‑1 (MQ=255)                                                                                                                                gcaaaTGAAGCCCTGCATCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGc                      <  2:780338/147‑1 (MQ=255)                                                                                                                                   atagaggCCCTCCTTCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCa                     <  2:1118699/142‑1 (MQ=255)                                                                                                                                      gAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACAcc                >  1:1112962/1‑150 (MQ=255)                                                                                                                                          cccTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGt             <  2:199897/149‑1 (MQ=255)                                                                                                                                          cccTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTAt           >  2:982665/1‑151 (MQ=255)                                                                                                                                          cccTGCACCCCCCCCGCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGt             <  2:879120/149‑1 (MQ=255)                                                                                                                                                aTCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCgg     >  2:222544/1‑151 (MQ=255)                                                                                                                                                aTCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCgg     >  1:1076607/1‑151 (MQ=255)                                                                                                                                                     caccTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGaaa  >  1:830965/1‑149 (MQ=255)                                                                                                                                                     caccTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGaaa  >  2:857010/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GAGCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGAAA  >  NC_000913/2511623‑2511919

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