Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,594,068	A→C	13.7%	intergenic (‑266/+35)	PP_1399 ← / ← kgtP	hypothetical protein/alpha‑ketoglutarate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,594,068	0	A	C	13.7%	188.6 / 24.5	73	intergenic (‑266/+35)	PP_1399/kgtP	hypothetical protein/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base C (7/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

TGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCC  >  NC_002947/1593929‑1594213                                                                                                                                            |                                                                                                                                                    tGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGgcgc                                                                                                                                            <  8:321013/149‑1 (MQ=255)   tCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCACCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGcc‑‑c                                                                                                                                        >  2:42464/1‑149 (MQ=255)     gTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGcc‑‑cc                                                                                                                                       <  7:31007/148‑1 (MQ=255)         ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCttt                                                                                                                                    >  5:188293/1‑147 (MQ=255)         ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTa                                                                                                                                  >  7:159985/1‑149 (MQ=255)         ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTa                                                                                                                                  >  7:53658/1‑149 (MQ=255)            aGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTa                                                                                                                               >  8:226930/1‑149 (MQ=255)            aGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGCATGGAgcgc                                                                                                                                            >  5:16229/1‑138 (MQ=255)            aGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGCATGGAgcgc                                                                                                                                            <  6:16229/138‑1 (MQ=255)             ggCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTAc                                                                                                                              >  6:251775/1‑149 (MQ=255)              gCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACt                                                                                                                             <  5:206344/149‑1 (MQ=255)                 gCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTgcc                                                                                                                          <  8:159985/149‑1 (MQ=255)                 gCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTgc                                                                                                                           >  2:73971/1‑148 (MQ=255)                   aaGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGATTTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTgccag                                                                                                                        >  3:335385/1‑149 (MQ=255)                       ctggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTgccaggcca                                                                                                                    <  5:251775/149‑1 (MQ=255)                        tggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTgccaggccag                                                                                                                   >  2:291956/1‑149 (MQ=255)                         ggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGc                                                                                                                  <  3:332785/149‑1 (MQ=255)                           ctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGt                                                                                                                >  5:206049/1‑149 (MQ=255)                               ggAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAATAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTAc                                                                                                             <  6:247535/148‑1 (MQ=255)                               agaacGTTGGCGCCCGCGCGAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACt                                                                                                            >  7:210316/5‑149 (MQ=255)                                gAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACt                                                                                                            >  7:243556/1‑148 (MQ=255)                                 aGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCa                                                                                                          <  1:42464/149‑1 (MQ=255)                                 aGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCa                                                                                                          <  5:200381/149‑1 (MQ=255)                                    gTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                       <  1:73971/149‑1 (MQ=255)                                     ttCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTgg                                                                                                      <  8:53658/149‑1 (MQ=255)                                     ttCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTgg                                                                                                      >  8:157274/1‑149 (MQ=255)                                     ttCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTgg                                                                                                      >  8:152285/1‑149 (MQ=255)                                                 ccAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGt                                                                                          <  7:157274/149‑1 (MQ=255)                                                 ccAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGt                                                                                          <  7:152285/149‑1 (MQ=255)                                                  cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCa                                                                                                          <  5:43908/132‑1 (MQ=255)                                                  cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCa                                                                                                          >  6:43908/1‑132 (MQ=255)                                                  cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                       <  8:344609/135‑1 (MQ=255)                                                  cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                       >  7:344609/1‑135 (MQ=255)                                                  cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtg                                                                                         >  2:49019/1‑149 (MQ=255)                                                   aaCGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCt                                                                                                                                      <  1:60735/103‑1 (MQ=255)                                                   aaCGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCt                                                                                                                                      >  2:60735/1‑103 (MQ=255)                                                     cGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgctg                                                                                      <  7:226930/149‑1 (MQ=255)                                                           gccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                                >  6:294583/1‑149 (MQ=255)                                                            ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgtg                                                                               <  4:335385/149‑1 (MQ=255)                                                            ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                                <  5:294583/148‑1 (MQ=255)                                                            ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                                <  8:249542/148‑1 (MQ=255)                                                                ccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCTCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTc                                                                             <  1:291956/149‑1 (MQ=255)                                                                      cAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATg                                                                      <  6:206049/148‑1 (MQ=255)                                                                            gATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtg                                                                                         <  8:167101/123‑1 (MQ=255)                                                                            gATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtg                                                                                         >  7:167101/1‑123 (MQ=255)                                                                              tAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTa                                                              >  3:140078/1‑148 (MQ=255)                                                                                  gTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAAc                                                           >  3:264308/1‑147 (MQ=255)                                                                                     gAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagca                                                      <  8:243556/149‑1 (MQ=255)                                                                                           aTAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTg                                                                                  <  4:315383/115‑1 (MQ=255)                                                                                                 cGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGa                                          <  6:188293/149‑1 (MQ=255)                                                                                                 cGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGa                                          >  7:13996/1‑149 (MQ=255)                                                                                                   ttAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATa                                        <  3:310983/149‑1 (MQ=255)                                                                                                        aagGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAg                                                 >  1:114520/1‑135 (MQ=255)                                                                                                        aagGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAg                                                 <  2:114520/135‑1 (MQ=255)                                                                                                          ggAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCgg                                 <  4:140078/149‑1 (MQ=255)                                                                                                            aaaaGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                                <  4:109668/148‑1 (MQ=255)                                                                                                             cctctCTGAATCGAACAGGGCCGCTCCATGCAAGGGGGCGCC‑‑CCTTTTAGTACTGCCACGCCAGCGCTACTCCGACGCAATCCGCGAGTGCCGCTTGGTGTCTTTCCTGGTCGCGTAAACCAGCAGCGAGCAGGCGACACAGGCGGTGa                              <  8:318957/145‑1 (MQ=255)                                                                                                              aaGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGAc                             <  1:49019/149‑1 (MQ=255)                                                                                                                   tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGa                              >  7:318957/1‑143 (MQ=255)                                                                                                                   tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                         >  6:325012/1‑148 (MQ=255)                                                                                                                     aaTAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCGGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTGCCAg                      >  1:7218/1‑149 (MQ=255)                                                                                                                     aaTAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                         <  5:325012/146‑1 (MQ=255)                                                                                                                        ctaaaaGGGGCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                         <  6:180545/141‑1 (MQ=39)                                                                                                                        ctaaaaGGGGCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                         >  5:180545/3‑143 (MQ=39)                                                                                                                          aaaaGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAg                   >  2:201115/1‑147 (MQ=255)                                                                                                                           aaaGGGGCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAg                >  4:171885/1‑149 (MQ=255)                                                                                                                             aGGGGCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGcc              >  4:24378/1‑149 (MQ=255)                                                                                                                              ggggCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCg             <  3:24378/149‑1 (MQ=255)                                                                                                                                ggCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGct           >  3:42649/1‑149 (MQ=255)                                                                                                                                ggCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGct           >  7:325065/1‑149 (MQ=255)                                                                                                                                ggCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGct           >  5:18034/1‑149 (MQ=255)                                                                                                                                 gcgcTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCATCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGctc          >  7:142030/1‑149 (MQ=255)                                                                                                                                    cTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGTCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCtctc        >  4:361339/1‑148 (MQ=255)                                                                                                                                            aaTGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATccc  <  3:326117/146‑1 (MQ=255)                                                                                                                                            aaTGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATccc  >  4:326117/1‑146 (MQ=255)                                                                                                                                            |                                                                                                                                                    TGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCC  >  NC_002947/1593929‑1594213

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