Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,852,105	G→T	K383N (AAG→AAT)	relA →	ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,852,105	0	G	T	96.0%	301.7 / ‑0.7	99	K383N (AAG→AAT)	relA	ATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (50/45);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

ACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGC  >  NC_002947/1851959‑1852241                                                                                                                                                   |                                                                                                                                          acTGCGGTAATCGGACCCGAGGGCAAGGTGCTGGAGGTGCAGGTCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATAtt                                                                                                                                        >  2:1372816/1‑150 (MQ=255)   tGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCATATGCGCCCACTGGCGCTACAAGGGCACCGACGTCACGCCCAGCTCCAACCACTACGAAGAAAATATTTc                                                                                                                                      >  1:573473/1‑150 (MQ=255)      ggTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTg                                                                                                                                   >  2:1969990/1‑150 (MQ=255)      ggTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTg                                                                                                                                   <  1:85981/150‑1 (MQ=255)       gTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGTATGACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTgg                                                                                                                                  >  1:1987219/1‑150 (MQ=255)           tCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGc                                                                                                                              >  2:258314/1‑150 (MQ=255)            cGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCg                                                                                                                             >  1:259858/1‑150 (MQ=255)              gCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTc                                                                                                                           >  1:544335/1‑150 (MQ=255)                 ccGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAgg                                                                                                                        >  1:1102166/1‑150 (MQ=255)                 ccGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAgg                                                                                                                        >  1:1686461/1‑150 (MQ=255)                 ccGAGGGCAAGGTGCTGGAGGTGCAGATCCGCACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAGGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAgg                                                                                                                        >  2:596679/1‑150 (MQ=255)                      ggCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTg                                                                                                                   <  1:1969990/150‑1 (MQ=255)                        cAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTGCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGa                                                                                                                 >  1:1016486/1‑150 (MQ=255)                        cAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGa                                                                                                                 >  2:719577/1‑150 (MQ=255)                        cAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGa                                                                                                                 <  1:1937975/150‑1 (MQ=255)                          aGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGt                                                                                                               <  1:198131/150‑1 (MQ=255)                           ggTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTg                                                                                                              >  2:869808/1‑150 (MQ=255)                                tGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACg                                                                                                         >  1:803107/1‑150 (MQ=255)                                 ggAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGa                                                                                                        >  1:494119/1‑150 (MQ=255)                                   aGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑ag                                                                                                     <  2:2031890/150‑1 (MQ=255)                                     gTGCAGATTCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGGACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTg                                                                                                                   >  2:1184864/1‑135 (MQ=255)                                     gTGCAGATTCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTg                                                                                                                   <  1:1184864/135‑1 (MQ=255)                                     gTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑agag                                                                                                   <  2:494119/150‑1 (MQ=255)                                      tGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCATGGTTGGGTCAGGTGCTGGAGTGGCACGA‑AGAGc                                                                                                  >  1:457791/1‑150 (MQ=255)                                      tGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGAACTAAGGTCAAGCCCAGCTGCAACTACGACGTAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGc                                                                                                  <  1:258314/150‑1 (MQ=255)                                          gATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTggg                                                                                              <  2:1780638/150‑1 (MQ=255)                                              cGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGAc                                                                                          >  1:2032748/1‑150 (MQ=255)                                              cGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCAAACCGCAAAATTCTTCGAAGAAAGGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGAc                                                                                          >  1:1882006/1‑150 (MQ=255)                                                   ccACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTTCAATTCATGAGTGGCACGA‑AGAGCTGGGCGACATCgg                                                                                     >  2:700940/1‑150 (MQ=255)                                                   ccACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCgg                                                                                     <  2:1968923/150‑1 (MQ=255)                                                     aCGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑ag                                                                                                     >  1:1359665/1‑132 (MQ=255)                                                     aCGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑ag                                                                                                     <  2:1359665/132‑1 (MQ=255)                                                         cATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGtggcct                                                                               <  2:573473/150‑1 (MQ=255)                                                           tGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGtggcctgg                                                                             <  1:635363/150‑1 (MQ=255)                                                            gCACGAAGAGGCCGAACTTGGCGTATGCGTCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGGGCTGGGCGACATCGAtggcctggc                                                                            >  2:764989/1‑150 (MQ=255)                                                            gCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATTGGtggcctggc                                                                            <  2:259858/150‑1 (MQ=255)                                                              aCGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCg                                                                          >  1:845135/1‑150 (MQ=255)                                                                gAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAg                                                                        >  1:1177406/1‑150 (MQ=255)                                                                    agGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCGCCAACCACTACGAAGTAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGGGCTGGGCGACATCGGTGGCCTGGCCGAGCAGt                                                                    <  2:213660/150‑1 (MQ=255)                                                                      gCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTg                                                                  >  2:1783794/1‑150 (MQ=255)                                                                       ccGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGc                                                                 >  2:467245/1‑150 (MQ=255)                                                                         gAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCgg                                                               <  1:467245/150‑1 (MQ=255)                                                                         gAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCgg                                                               >  2:474715/1‑150 (MQ=255)                                                                         gAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCgg                                                               <  2:561846/150‑1 (MQ=255)                                                                            cTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTc                                                            <  2:2032748/150‑1 (MQ=255)                                                                            cTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTc                                                            <  2:1987219/150‑1 (MQ=255)                                                                            cTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTc                                                            >  2:2101780/1‑150 (MQ=255)                                                                            cTTGGCGTATGCGCCCACTGGCGATCCAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTc                                                            >  2:1305608/1‑150 (MQ=255)                                                                             ttGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCgg                                                               <  1:2015744/146‑1 (MQ=255)                                                                             ttGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCgg                                                               >  2:2015744/1‑146 (MQ=255)                                                                             ttGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCg                                                           >  1:61489/1‑150 (MQ=255)                                                                                gCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACa                                                        >  2:1549545/1‑150 (MQ=255)                                                                                   caTGCGCCCACTGGAGCTAAAAGCAACCACCCGTCCAGCCCAGCCACCACCCCCACCCAGCAAATATTTCCTGGTTGCGTCAGGTTCGGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGAcctag                                                     <  1:1372816/149‑5 (MQ=255)                                                                                      gcgcCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGc                                                  <  2:61489/150‑1 (MQ=255)                                                                                        gcCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGAGCGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCCGGGCGAGCAGCGGTGACTCTACATCGAGCCt                                                >  1:2109556/1‑150 (MQ=255)                                                                                          ccACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGa                                              >  2:1201634/1‑150 (MQ=255)                                                                                               ggCGCTACAAGGGCACCGACGTCAAGCCCAGCTGCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCggg                                         <  2:1016486/150‑1 (MQ=255)                                                                                               ggCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCggg                                         <  2:1177406/150‑1 (MQ=255)                                                                                               ggCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCggg                                         <  1:719577/150‑1 (MQ=255)                                                                                                    tACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGTCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTAt                                    <  2:1102166/150‑1 (MQ=255)                                                                                                    tACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTAt                                    >  1:2111331/1‑150 (MQ=255)                                                                                                    tACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTAt                                    <  1:764989/150‑1 (MQ=255)                                                                                                    tACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTAt                                    <  2:1882006/150‑1 (MQ=255)                                                                                                      cAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGa                                                                         <  2:1630648/111‑1 (MQ=255)                                                                                                      cAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGa                                                                         >  1:1630648/1‑111 (MQ=255)                                                                                                       aaGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTAtgtg                                 <  1:869808/150‑1 (MQ=255)                                                                                                         gggCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGa                                                                                           >  1:624681/1‑90 (MQ=255)                                                                                                         gggCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGa                                                                                           <  2:624681/90‑1 (MQ=255)                                                                                                         gggCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGtt                               >  2:183726/1‑150 (MQ=255)                                                                                                         gggCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTACGGGTCGACATCGAGCCTGACCGGGTTTATGTGtt                               >  1:1009548/1‑150 (MQ=255)                                                                                                              ccGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCAccc                          <  1:1958205/150‑1 (MQ=255)                                                                                                               cGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCAcccc                         >  1:922289/1‑150 (MQ=255)                                                                                                                  cGTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGa                      <  2:922289/150‑1 (MQ=255)                                                                                                                   gTCAAGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGAc                     <  2:544335/150‑1 (MQ=255)                                                                                                                      aaGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCa                             <  1:37274/139‑1 (MQ=255)                                                                                                                      aaGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCa                             >  2:37274/1‑139 (MQ=255)                                                                                                                      aaGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGc                  >  1:1824621/1‑150 (MQ=255)                                                                                                                      aaGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGc                  >  2:1651585/1‑150 (MQ=255)                                                                                                                      aaGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGc                  >  2:54004/1‑150 (MQ=255)                                                                                                                       aGCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGcc                 <  2:1189832/150‑1 (MQ=255)                                                                                                                        gCCCAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCa                >  2:2144758/1‑150 (MQ=255)                                                                                                                         cccAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGCCCGGGTTTCTGTGTTCACCTCCGACGGCCAc               >  2:148872/1‑150 (MQ=255)                                                                                                                         cccAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGAAATCTGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCAc               <  1:786912/150‑1 (MQ=255)                                                                                                                         cccAGCTCCAACCACTACGAAGAAAATATTTCCTGGTTGAGTCAGGTGCTGGAGTGGCACGA‑AGAGCCGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGACGACAACGAGCCGGACCCGGTTTATGTGTTCACCCCCAACGGCCAc               >  1:244346/1‑150 (MQ=255)                                                                                                                          ccAGCTTCAAACAGTAGGAAGAAAATATTGCCTGGTTGCGTCAGGTGCTGGAGTACCACCATGGAAATGGG‑GGTAGCAGCGGCCTGGCTGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACg              <  1:148872/150‑1 (MQ=255)                                                                                                                                 cAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGc                                                  >  2:1697428/1‑107 (MQ=255)                                                                                                                                 cAACCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGc                                                  <  1:1697428/107‑1 (MQ=255)                                                                                                                                  aaCCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGTGTTTATGTGTTCACCCCCGACGGCCACGCCATCGAc      <  1:596679/150‑1 (MQ=255)                                                                                                                                  aaCCACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGAc      <  1:2101780/150‑1 (MQ=255)                                                                                                                                    ccACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCt    <  1:1201634/150‑1 (MQ=255)                                                                                                                                    ccACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCt    >  1:1135429/1‑150 (MQ=255)                                                                                                                                    ccACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCt    >  1:1059510/1‑150 (MQ=255)                                                                                                                                     cACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGc                                                  >  1:406458/1‑103 (MQ=255)                                                                                                                                     cACTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGc                                                  <  2:406458/103‑1 (MQ=255)                                                                                                                                      aCTACGAAGAAAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTACGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGc  <  2:1009548/150‑1 (MQ=255)                                                                                                                                           gaagaaAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACAt                                                       <  1:811764/92‑1 (MQ=255)                                                                                                                                           gaagaaAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACAt                                                       >  2:811764/1‑92 (MQ=255)                                                                                                                                           gaagaaAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGc  >  2:1770267/1‑145 (MQ=255)                                                                                                                                           gaagaaAATATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGc  <  1:1770267/145‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                          ACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGA‑AGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGC  >  NC_002947/1851959‑1852241

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