Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,915,514	C→T	6.1%	intergenic (‑314/+39)	glmU ← / ← atpC	fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/F1 sector of membrane‑bound ATP synthase, epsilon subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,915,514	0	C	T	6.1%	233.4 / 4.1	84	intergenic (‑314/+39)	glmU/atpC	fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/F1 sector of membrane‑bound ATP synthase, epsilon subunit
Reads supporting (aligned to +/- strand):  ref base C (42/35);  new base T (3/2);  total (47/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

ATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTAATGTGCT  >  NC_000913/3915366‑3915662                                                                                                                                                      |                                                                                                                                                     aTTAAGTCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAagccagc                                                                                                                                                      <  1:718189/149‑1 (MQ=255) aTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAagccagc                                                                                                                                                      >  2:731895/1‑149 (MQ=255)  ttAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAagccagcc                                                                                                                                                     >  2:1537245/1‑149 (MQ=255)     aGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGt                                                                                                                                                  <  2:587139/149‑1 (MQ=255)             ccTTAATTTAAGCGGCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGAc                                                                                                                                          <  2:682788/149‑1 (MQ=255)               ttAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTg                                                                                                                                        >  2:1508921/1‑149 (MQ=255)               ttAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCACCCTTTTTCCAGACTg                                                                                                                                        >  2:1423223/1‑149 (MQ=255)               ttAATTTAAGCGTCAGGTGGATGTTTTAGCTCTAATTTCGATCAATGAATAAACAGGAAAACAGGGGTGTTTTATATCAGCATGATCTATGTGAACGCAAGTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCGAGAcgg                                                                                                                                        >  2:177479/1‑147 (MQ=255)                 aaTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGc                                                                                                                                      <  1:306085/149‑1 (MQ=255)                 aaTTTAAGCGTCAGCTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTCAACGCTAATCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGCCAGGc                                                                                                                                      >  2:659176/1‑149 (MQ=255)                    ttAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCttt                                                                                                                                   >  2:1524224/1‑149 (MQ=255)                       aGCGTCAGGTGTATGTTTTTGCTCTTATTTCGCTCAATGAATAAACAGGATAACACGGGGGTTTTATATCAGCAGGTTCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCTGCCTGTTTCCAGACTGGCTTTtgt                                                                                                                                <  1:1173906/149‑1 (MQ=255)                        gCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTtgtg                                                                                                                               <  1:721983/149‑1 (MQ=255)                         cGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAATCCAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGTCTGGAAACAGGCTGGCtttttttt                                                                                                                              <  1:823349/149‑5 (MQ=255)                          gTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCt                                                                                                                             >  2:1093343/1‑149 (MQ=255)                          gTCAGGTGGATGTTTTTCCTCTTATTTCGATCAATTAATAAACTGTAAAACAGGGGGGTTTTATATCAGCTGTATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGTCTGGAAACAGGCTGGCttttttttg                                                                                                                             <  2:185025/149‑6 (MQ=37)                             aGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCtttt                                                                                                                          <  1:583501/149‑1 (MQ=255)                             aGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCtttt                                                                                                                          <  1:1439432/149‑1 (MQ=255)                                   aTGTTTTTGCTCTTATTTCGATCAATGACTAAACAGGACAACAGGGGGGTTTTATATCAGCACGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGcc                                                                                                                    <  2:881549/149‑1 (MQ=255)                                     gtgtttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATAGCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCgg                                                                                                                  <  1:1111446/146‑1 (MQ=255)                                     gTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCgg                                                                                                                  <  1:368286/149‑1 (MQ=255)                                     gTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGATTTTGTGCTTTTCAAGCCgg                                                                                                                  >  2:309118/1‑149 (MQ=255)                                      tttttGCTCCTATTTCGCTCAATGACTCAACAGGACAACAGGGGGGTTTTATATCAGCAGGATCTCTGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGgt                                                                                                                 <  1:1067097/149‑1 (MQ=255)                                        tttGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑GGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGtt                                                                                                              <  1:524323/149‑1 (MQ=255)                                         ttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGtt                                                                                                              <  1:748770/149‑1 (MQ=255)                                           gCTCTTATTTCGATCAATGAATAAACTGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATTTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCTGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCTGTGTTAc                                                                                                            <  1:908489/149‑1 (MQ=255)                                             tctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACAt                                                                                                          >  1:1362642/1‑149 (MQ=255)                                             tctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACAt                                                                                                          >  1:1263056/1‑149 (MQ=255)                                             tctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACAt                                                                                                          >  2:106478/1‑149 (MQ=255)                                                   ttCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTAGATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGGTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCttt                                                                                                    >  2:846568/1‑149 (MQ=255)                                                         cAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGt                                                                                              <  1:1317954/149‑1 (MQ=255)                                                          actGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTc                                                                                             <  1:1117357/147‑1 (MQ=255)                                                            tGAATAAACACGAACACAGGGGGGTTTTATATCACCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCaa                                                                                           <  2:24999/149‑1 (MQ=255)                                                             gAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAAc                                                                                          >  2:1521793/1‑149 (MQ=255)                                                                  aaCAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGAt                                                                                     >  2:233596/1‑149 (MQ=255)                                                                  aaCAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGACTGTTTCCAGACTGGCTTTTGTGCTATTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGAt                                                                                     >  1:961778/1‑149 (MQ=255)                                                                       gacaaCAGCGGTGTTTTATATCACCCTGATCTTTGTGAACGCGATTCAGGACGGGTCCCACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAAcgc                                                                                <  2:1060019/146‑1 (MQ=255)                                                                       gAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGACAGCTATTCTGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAAcgc                                                                                <  1:379031/149‑1 (MQ=255)                                                                       gAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAAcgc                                                                                >  1:1530449/1‑149 (MQ=255)                                                                        aaaaCAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTGTTTTGGTCAACTCGATAAcgcg                                                                               >  2:933360/1‑149 (MQ=255)                                                                          aaCAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGCCTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCa                                                                             >  1:1419582/1‑149 (MQ=255)                                                                                 gggTTTTATATCAGAAGGATCTATGTGAACGCTATTCAGGACGGGGCACTCGCGC‑AAAAAAAAGCCAGCCTGTTTCGAGACTGGCTATTGTGCTTTTGAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTgcgc                                                                      <  2:1391156/149‑1 (MQ=255)                                                                                    ttttATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGGCTGTTTCCAGACTGGCTTTTGTGCTTATCAAGCCTGTGTTACATCGCTTTTTTGGTCAACTCGTTTACGCGCAGCTGCGCGAt                                                                   >  1:318070/1‑149 (MQ=255)                                                                                    ttttATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGAt                                                                   >  2:468639/1‑149 (MQ=255)                                                                                     tttATATCAGCAGGATCTATGTGAACGCTATTGAGGACGGGGCACACGCGC‑AAAAAAATGCCAGCCTGTTTCCACACTGGCTTTTGTGCTTTTCAAGCCGGTTTTACATCGCTTTTTTGGTCAAGTCGATAACGCGCAGCTGCGCGTTc                                                                  <  1:493687/149‑1 (MQ=255)                                                                                     tttATAGCACCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTATCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATc                                                                  <  1:1054793/149‑1 (MQ=255)                                                                                     tttAAATCAGCAGGATATATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTAAACTCGATAACGCGCAGCTGCGCGATc                                                                  <  1:87522/149‑1 (MQ=255)                                                                                         tatCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCtt                                                              >  2:46937/1‑149 (MQ=255)                                                                                          atCAGCAGGATCTATGTGAACGCTATTCAGGACGGGGCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCttt                                                             <  2:309522/149‑1 (MQ=255)                                                                                           tCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCCGCCTGTTTCCAGACTGGCTTTTGTGCGTTTCAAGCCGGTGTTACATCGCTTTTTTGGTAAACTCGATAACGCGCAGCTGCGCGATCGCTTTg                                                            >  2:1538667/1‑149 (MQ=255)                                                                                            cagcctgATCTATGTGCACGCTAATCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCGGCGCGATCGCTTTgg                                                           <  2:321575/143‑1 (MQ=255)                                                                                             cgcagGATCTATGTGAACGCTATTCAGGACGGGTCCCACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGc                                                          <  2:673696/148‑1 (MQ=255)                                                                                             agcagGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGt                                                                                              <  1:1081994/113‑1 (MQ=255)                                                                                             agcagGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGt                                                                                              >  2:1081994/1‑113 (MQ=255)                                                                                             agcagGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTGTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGc                                                          >  1:562134/1‑149 (MQ=255)                                                                                              gcagGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTATTCACGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAAAGCGCAGCTGCGCGATCGCTTTGGcc                                                         >  1:341257/1‑149 (MQ=255)                                                                                               cagGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCa                                                        <  2:587010/149‑1 (MQ=255)                                                                                                 ggATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCATTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGt                                                      <  1:618655/149‑1 (MQ=255)                                                                                                  gATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGtt                                                     >  2:1196496/1‑149 (MQ=255)                                                                                                       aTGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCa                                                >  1:1544457/1‑149 (MQ=255)                                                                                                       aTGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGCTCAACTCGATAACGCTCAGCTGCGCGATCGCTTTGGCCAGTTCCGCa                                                >  2:1164226/1‑149 (MQ=255)                                                                                                        tgtgAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAg                                               >  1:708744/1‑149 (MQ=255)                                                                                                          tgAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGAc                                             >  1:1414141/1‑149 (MQ=255)                                                                                                          tgAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTATTCAAGCCGGTGTTACATCGCTTTTTTCGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGCCCAGTTCCGCAGAc                                             >  1:179992/1‑149 (MQ=255)                                                                                                           gAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACg                                            >  1:1305004/1‑149 (MQ=255)                                                                                                                cTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCATACGCCTGa                                       <  1:5927/149‑1 (MQ=255)                                                                                                                 tATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAg                                      <  1:1537245/149‑1 (MQ=255)                                                                                                                  aTTCTGGACGGGTCACACGGGC‑AAAAAAAAGCCAGCCTGTTACCAGACTGGCTTTTGTGCTTTTCAAGCCGGAGTTACATCGCTTTTTTGGGCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGc                                     <  1:1187558/149‑1 (MQ=255)                                                                                                                  aTTCAGGACGGGGCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGc                                     <  2:9003/149‑1 (MQ=255)                                                                                                                   ttCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCATCTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCg                                    >  2:1252462/1‑149 (MQ=255)                                                                                                                      aGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGCCTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTATTTGGTCAACTCGATAACGCGCAGCTGCGCGCTCGCTTGTGCCAGTTCCGCAGACGCCTGCGCGTaa                                 >  1:874047/1‑149 (MQ=255)                                                                                                                      aGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTATTGTGCTTATCAAGCCGGTGTTACATCGCTTATTTGGTCAACTAGAAAACGCGCCCCTGCTCGATACCTTTGGCCACTTCCGCACACGCATGAGCGTaa                                 >  1:91164/1‑149 (MQ=255)                                                                                                                       ggACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTTCTTTTCAAGCCGGTGGTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCAAGTTCCGCAGACGCCTGAGCGTact                                >  1:865477/1‑147 (MQ=255)                                                                                                                       ggACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACAGGCGTTTGTGCTTTTCAAGCCGGGGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCAGCGCGATAGCTTTGGCCAGTTCCGAACAAGCCTGAGCGTaaa                                >  2:1407372/1‑148 (MQ=255)                                                                                                                        gACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCACCTCGATAACGCGCCGCTGCGCGCTCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATc                               >  1:1513199/1‑149 (MQ=255)                                                                                                                              tCACACGCGCAAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTATTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGt                          >  1:311028/1‑149 (MQ=255)                                                                                                                              tCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTc                         >  1:1401104/1‑149 (MQ=255)                                                                                                                                acacGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGc                       >  1:959501/1‑149 (MQ=255)                                                                                                                                 cccgcgC‑TAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGcc                      <  1:134264/147‑1 (MQ=255)                                                                                                                                  acGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCg                     >  2:655123/1‑149 (MQ=255)                                                                                                                                    gcgc‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTg                   >  2:697044/1‑149 (MQ=255)                                                                                                                                       c‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTgaga                >  1:825259/1‑149 (MQ=255)                                                                                                                                         aaaaaaaaGCAAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTATATGCTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTCagag               >  2:637614/1‑149 (MQ=255)                                                                                                                                          cacaaaaGCCAGTCTGGAAACAGGCTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAgc              >  1:769570/4‑149 (MQ=255)                                                                                                                                           acaaaaGCCAGTCTGGAAACAGGCTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAgct             >  1:114743/3‑149 (MQ=255)                                                                                                                                             aaaaGCCAGTCTGGAAACAGGCTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCTATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACTCCTGACCGTAATCTACGTCGCCGTGAGAgctgc           >  2:263355/1‑149 (MQ=255)                                                                                                                                             aaaaGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAgctgc           <  2:1414141/149‑1 (MQ=255)                                                                                                                                                agccagccTGTTTCCTGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTGGATAACGCTCAGCTGCGCGATCGCTTAGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTaa        <  2:341257/149‑1 (MQ=255)                                                                                                                                                agccagccTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTaa        <  1:731895/149‑1 (MQ=255)                                                                                                                                                 gccagccTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTAAt       <  1:106478/149‑1 (MQ=255)                                                                                                                                                 gccagccTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTAAt       <  1:1508921/149‑1 (MQ=255)                                                                                                                                                  ccagccTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTAAtg      >  2:1479921/1‑149 (MQ=255)                                                                                                                                                   cagccTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTAAtgt     >  1:998211/1‑149 (MQ=255)                                                                                                                                                   cagccTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTATTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTAAtgt     <  1:939249/149‑1 (MQ=255)                                                                                                                                                     gccTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTAATGTGc   <  2:1401104/149‑1 (MQ=255)                                                                                                                                                      ccTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGCCGCCTGAGCGGAATCTACGTCGCCGTGAGAGCTGCTAATGTGca  >  1:278830/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                     ATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGC‑AAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTAATGTGCT  >  NC_000913/3915366‑3915662

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