Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,671,863	G→A	6.3%	A152V (GCC→GTC)	acnB ←	aconitate hydratase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,671,863	0	G	A	6.3%	186.2 / 2.4	63	A152V (GCC→GTC)	acnB	aconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base A (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

AATGTCAGGGCGCGACCAGGCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGCCCAGTTCGGCGTCATCCAGCAGTG  >  NC_002947/2671736‑2672007                                                                                                                                 |                                                                                                                                                  aaTGTCAGGGCGCGACCAGGCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTCCTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGAGGCCA‑AGGATTGCTGGActg                                                                                                                             >  4:76215/1‑147 (MQ=255)    gTCAGGGCGCGACCAGGCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCt                                                                                                                          >  7:139/1‑149 (MQ=255)         ggCGCGACCAGGCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGgcg                                                                                                                     >  2:322347/1‑149 (MQ=255)          gcgcGACCAGGCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGgcgt                                                                                                                    >  5:94029/1‑149 (MQ=255)          gcgcGACCAGGCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGgcgt                                                                                                                    <  8:107211/149‑1 (MQ=255)            gcgACCAGGCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGgcgtgc                                                                                                                  <  3:181901/149‑1 (MQ=255)                 cAGGCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGgcgtgc                                                                                                                  <  1:94624/144‑1 (MQ=255)                 cAGGCGTCCGGGGCAGGCGACAGGTAGTCGGAGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGgcgtgc                                                                                                                  >  2:94624/1‑144 (MQ=255)                  aGGCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTg                                                                                                            >  8:23910/1‑149 (MQ=255)                    gCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGcc                                                                                                          >  5:88277/1‑149 (MQ=255)                    gCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGcc                                                                                                          >  3:152962/1‑149 (MQ=255)                         cGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGcc                                                                                                      <  2:203243/148‑1 (MQ=255)                           gggCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTg                                                                                                   >  2:253627/1‑149 (MQ=255)                            ggCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCGTTGGCGTGCGTGTTGCCGGCCTTgg                                                                                                  <  4:1653/149‑1 (MQ=255)                             gCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGc                                                                                                 >  3:284503/1‑149 (MQ=255)                             gCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AAGATTCCAGGACGGCCTTGGCGTGCTTGTTGCCGGCCTTGGc                                                                                                 <  4:96157/149‑1 (MQ=255)                                  cGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCttttt                                                                                            <  1:322347/149‑1 (MQ=255)                                  cGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCttttt                                                                                            <  6:88277/149‑1 (MQ=255)                                   gagaGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCCGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑ATGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCtt                                                                                                    <  1:64411/137‑1 (MQ=255)                                   gaaaGGTCGTCGGTTTTGGTTTCGCCAGTAACCTTGAACACGGTCAGGCTGTACTTGTCGGCGATGGCCTGGCGCG‑AGGTGAAAAACTCGCCGGCAGACC‑AGGAATCCAGGACGGCCTTGGCGTGAGTGTTGCCCGCCta                                                                                                    >  2:64411/4‑139 (MQ=255)                                    aCAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCg                                                                                          <  6:94029/149‑1 (MQ=255)                                      aGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGc                                                                                        <  5:246022/149‑1 (MQ=255)                                        gtcgtcGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTg                                                                                                   <  6:245134/136‑1 (MQ=255)                                        gtcgtcGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTg                                                                                                   >  5:245134/1‑136 (MQ=255)                                           gtcgGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCAcgt                                                                                   >  7:200071/1‑149 (MQ=255)                                                tgtTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTg                                                                               >  3:265265/1‑148 (MQ=255)                                                  ttGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCCTCGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTgg                                                                              <  1:253627/148‑1 (MQ=255)                                                    ggTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAgg                                                                          >  4:222797/1‑149 (MQ=255)                                                      tttCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAGCCACTCGCCGGCAGCCC‑CGGAGTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGTAggga                                                                        >  4:247443/1‑147 (MQ=255)                                                        tCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCAt                                                                       <  8:139/148‑1 (MQ=255)                                                           ccAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTgg                                                                                                                       >  2:67009/1‑97 (MQ=255)                                                           ccAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGaa                                                                   <  8:200071/149‑1 (MQ=255)                                                           ccAGGCACCTTGAACAAGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCGCAGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTgg                                                                                                                       <  1:67009/98‑1 (MQ=255)                                                            cAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGaa                                                                   >  1:209952/1‑148 (MQ=255)                                                              ggCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACAt                                                                <  2:209952/149‑1 (MQ=255)                                                                    ttGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAg                                                          >  7:286926/1‑149 (MQ=255)                                                                    ttGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGACTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAg                                                          >  5:299937/1‑149 (MQ=255)                                                                        acacGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCg                                                                     <  8:116912/134‑1 (MQ=255)                                                                        acacGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCg                                                                     >  7:116912/1‑134 (MQ=255)                                                                          acGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGtgtg                                                    >  7:243485/1‑149 (MQ=255)                                                                          acGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGtgtg                                                    >  2:73885/1‑149 (MQ=255)                                                                            ggTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCt                                                  >  4:224614/1‑149 (MQ=255)                                                                                aGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGa                                               <  8:286926/148‑1 (MQ=255)                                                                                aGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGa                                               <  1:129510/148‑1 (MQ=255)                                                                                        cTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCa                                                              >  3:217349/1‑125 (MQ=255)                                                                                        cTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCa                                                              <  4:217349/125‑1 (MQ=255)                                                                                        cTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCgg                                      <  3:224614/149‑1 (MQ=255)                                                                                           gTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGcc                                    <  1:73885/148‑1 (MQ=255)                                                                                                 gATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCtt                                                 >  7:209595/1‑129 (MQ=255)                                                                                                 gATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCtt                                                 <  8:209595/129‑1 (MQ=255)                                                                                                 gATGGCCGGGCGCG‑AGGGGAACCACTCGCCGGCAGCCC‑AGGATTCCCGGCCGGCCTTGGCGCGCGTGGTGCCGGCCTTGGCTTTTTCGGCCACCTCGGGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATgg                             <  8:51151/149‑1 (MQ=255)                                                                                                      ccGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGcc                         >  3:3422/1‑148 (MQ=255)                                                                                                      ccGGGCGCG‑AGGTGAACCACTCGCCGACAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGcc                         >  6:14975/1‑148 (MQ=255)                                                                                                      ccGGGCGCG‑AGGTGAACCACTCGCCGACAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGcc                         >  8:260051/1‑148 (MQ=255)                                                                                                          gcgcG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGCCCAGtt                    <  4:152962/149‑1 (MQ=255)                                                                                                           cgcg‑AGGTGAACCACTCGCCGACAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGCCCAGTTc                   <  7:260051/149‑1 (MQ=255)                                                                                                           cgcg‑AGGTGAACCACTCGCCGACAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGCCCAGTTc                   <  5:14975/149‑1 (MQ=255)                                                                                                            gcg‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGCCCAGTTCg                  <  6:299937/149‑1 (MQ=255)                                                                                                              g‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCGACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGCCCAGTTCGGc                >  1:103049/1‑149 (MQ=255)                                                                                                                       ccACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGCCCAGTTCGGCGTCATCca        >  4:46352/1‑149 (MQ=255)                                                                                                                             gCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGc                                                   >  8:216370/1‑100 (MQ=255)                                                                                                                             gCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGc                                                   <  7:216370/100‑1 (MQ=255)                                                                                                                             gCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGCCCAGTTCGGCGTCATCCAGCAGTg  >  7:120880/1‑149 (MQ=255)                                                                                                                                 |                                                                                                                                                  AATGTCAGGGCGCGACCAGGCGTCCGGGGCAGGCGACAGGTCGTCGGTGTTGGTTTCGCCAGGCACCTTGAACACGGTCAGGGTGTACTTGTCGGCGATGGCCGGGCGCG‑AGGTGAACCACTCGCCGGCAGCCC‑AGGATTCCAGGACGGCCTTGGCGTGCGTGTTGCCGGCCTTGGCTTTTTCGGCCACGTCGTGGAAGGCATCGAACATCAGTAGGGTGTGCTTGAGCTGTTCGGCCGCGATGGCGCCCAGTTCGGCGTCATCCAGCAGTG  >  NC_002947/2671736‑2672007

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