Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	349,946	A→G	7.2%	intergenic (+374/‑66)	prpB → / → prpC	2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	349,946	0	A	G	7.2%	203.5 / 5.6	69	intergenic (+374/‑66)	prpB/prpC	2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base G (3/2);  total (35/34)
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.95e-01

GCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGC  >  NC_000913/349798‑350091                                                                                                                                                     |                                                                                                                                                  gcgcCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATTCCTGATGCGACGCTGACGCGAATTATCAGGCATACGGGGTGAACAGAAATGTAGGTCGGATaa                                                                                                                                                  >  2:1334829/1‑150 (MQ=11)  cgccgcTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAg                                                                                                                                                 >  1:51478/1‑150 (MQ=34)    ccgcTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGAc                                                                                                                                               <  2:51478/150‑1 (MQ=34)         gACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATa                                                                                                                                          <  1:1161027/150‑1 (MQ=34)          aCGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATGCGACGctatg                                                                                                                                         <  1:112375/150‑5 (MQ=255)           cgcgACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATGCGACGctatgc                                                                                                                                        >  2:1244208/1‑145 (MQ=255)                       ggCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATGCGACGCTATGCGTCTTATCCGac                                                                                                                            <  1:1244208/150‑1 (MQ=11)                           tACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGacaaca                                                                                                                        >  1:110486/1‑150 (MQ=32)                              ggggTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACACCAAtc                                                                                                                     >  2:1547373/1‑150 (MQ=21)                                   gCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGAc                                                                                                                >  1:326705/1‑150 (MQ=39)                                   gCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGAc                                                                                                                <  2:1005141/150‑1 (MQ=39)                                     aCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGAccc                                                                                                              >  1:493154/1‑150 (MQ=40)                                     aCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGAccc                                                                                                              >  2:267568/1‑150 (MQ=40)                                        gAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTAc                                                                                                           >  2:1558831/1‑150 (MQ=255)                                         aaCTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACa                                                                                                          >  2:184670/1‑150 (MQ=255)                                             gTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAg                                                                                                                                                 >  1:547712/1‑107 (MQ=32)                                             gTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAg                                                                                                                                                 <  2:547712/107‑1 (MQ=32)                                             gTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATg                                                                                                      >  2:1198534/1‑150 (MQ=255)                                               tggTGGGATAAGACGAATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGAt                                                                                                    <  1:47176/149‑1 (MQ=255)                                               aGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGAt                                                                                                    <  2:1604071/150‑1 (MQ=255)                                                  tCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACa                                                                                                          >  2:846714/1‑141 (MQ=255)                                                  tCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACa                                                                                                          <  1:846714/141‑1 (MQ=255)                                                      aTAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATg                                                                                             >  1:1155952/1‑150 (MQ=255)                                                            cGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAgg                                                                                       >  1:905752/1‑150 (MQ=255)                                                               atcgCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGaca                                                                                    <  1:1115561/147‑1 (MQ=255)                                                                 aGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                                                                                                                    >  2:590083/1‑100 (MQ=33)                                                                 aGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                                                                                                                    <  1:590083/100‑1 (MQ=33)                                                                 aGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGacaac                                                                                  <  1:1368847/150‑1 (MQ=255)                                                                  gCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGacaaca                                                                                 >  2:1234013/1‑150 (MQ=255)                                                                      cGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAg                                                                             <  2:1252234/150‑1 (MQ=255)                                                                         aTCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGa                                                                          >  2:369125/1‑150 (MQ=255)                                                                           ccGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGaca                                                                        >  1:337676/1‑150 (MQ=255)                                                                               cAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAAc                                                                    <  1:1558831/150‑1 (MQ=255)                                                                               cAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAAc                                                                    >  2:1129397/1‑150 (MQ=255)                                                                                 aCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGa                                                                  >  2:310620/1‑150 (MQ=255)                                                                                  cagATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGAt                                                                 >  1:413488/3‑150 (MQ=255)                                                                                    gATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATcc                                                               >  1:348015/1‑150 (MQ=255)                                                                                      tGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTg                                                             <  1:1291671/150‑1 (MQ=255)                                                                                         cTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCaa                                                          <  1:1234013/150‑1 (MQ=255)                                                                                           gATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTg                                                             >  1:1383667/1‑145 (MQ=255)                                                                                           gATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTg                                                             <  2:1383667/145‑1 (MQ=255)                                                                                              gcgcCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAg                                                     >  2:682114/1‑150 (MQ=255)                                                                                               cgccgcTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGt                                                    <  2:348015/150‑1 (MQ=255)                                                                                                gccgcTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGac                                                                                     >  1:366301/1‑116 (MQ=255)                                                                                                gccgcTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGac                                                                                     <  2:366301/116‑1 (MQ=255)                                                                                                 ccgcTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTAc                                                  >  2:109539/1‑150 (MQ=255)                                                                                                   gcTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGTTAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTAccc                                                >  2:638420/1‑150 (MQ=255)                                                                                                        cgcgACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATACCAATGCCGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTc                                           >  2:157741/1‑150 (MQ=255)                                                                                                        cgcgACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTTCAAAACAGTACCCATGTc                                           <  2:413488/150‑1 (MQ=255)                                                                                                        cgcgACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTc                                           >  1:1631084/1‑150 (MQ=255)                                                                                                         gcgACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCa                                          <  2:1631084/150‑1 (MQ=255)                                                                                                            aCTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTa                                       <  2:905752/150‑1 (MQ=255)                                                                                                            aCTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTa                                       <  1:1198534/150‑1 (MQ=255)                                                                                                            aCTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTa                                       >  1:1203571/1‑150 (MQ=255)                                                                                                            aCTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTa                                       <  1:753606/150‑1 (MQ=255)                                                                                                                aTCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAAcc                                   <  1:682114/150‑1 (MQ=255)                                                                                                                aTCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAAcc                                   <  1:1473647/150‑1 (MQ=255)                                                                                                                aTCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAAcc                                   <  1:458165/150‑1 (MQ=255)                                                                                                                         aCGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATc                          <  2:337676/150‑1 (MQ=255)                                                                                                                            gggTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCtgt                       >  1:311744/1‑150 (MQ=255)                                                                                                                                        ctgTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTgg           <  1:310620/150‑1 (MQ=255)                                                                                                                                          gTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCg         <  1:267568/150‑1 (MQ=255)                                                                                                                                            ttgTCGGATGCGACGCTATGCGTCTTATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGGCACAACGATCCTGCAAAACAGTGCCCATGTCATTAAGCCGAAAAAATCTGTGGCGCTTTCTGGCGtg       >  2:50639/3‑149 (MQ=255)                                                                                                                                            ttgTCGGATGCGACGCTATGCGTCTTATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGtt       >  1:872601/3‑150 (MQ=255)                                                                                                                                             tgTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGaaaaa                             <  2:1228030/126‑1 (MQ=255)                                                                                                                                             tgTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGaaaaa                             >  1:1228030/2‑127 (MQ=255)                                                                                                                                              gTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTcc     <  1:184670/150‑1 (MQ=255)                                                                                                                                              gTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTcc     <  2:1155952/150‑1 (MQ=255)                                                                                                                                                 ggATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCcggc  >  1:58333/1‑150 (MQ=255)                                                                                                                                                 ggATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCcggc  <  2:873855/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAATGACGAGGACAACATGAGCGACACAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGC  >  NC_000913/349798‑350091

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