Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	157,877	G→A	5.9%	intergenic (‑152/‑113)	citN ← / → PP_0148	citrate transporter/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	157,877	0	G	A	5.9%	201.2 / 2.3	67	intergenic (‑152/‑113)	citN/PP_0148	citrate transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base A (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TCCAGGCGTGGCGCGGCGAAAGGACAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGC  >  NC_002947/157734‑158016                                                                                                                                                 |                                                                                                                                            tCCAGGCGTGGCGCGGCGAAAGGACAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGcgcg                                                                                                                                         >  4:98842/1‑148 (MQ=255) tCCAGGCGTGGCGCGGCGAAAGGACAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCgg                                                                                                                                        >  3:39360/1‑149 (MQ=255)   cAGGCGTGGCGCGGCGAAAGGACAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGcgcg                                                                                                                                         <  3:98842/146‑1 (MQ=255)   cAGGCGTGGCGCGGCGAAAGGACAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTa                                                                                                                                      <  3:108100/149‑1 (MQ=255)   cAGGCGTGGCGCGGCGAAAGGACAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGACGCGGTa                                                                                                                                      >  6:219849/1‑149 (MQ=255)    aGGCGTGGCGCGGCGAAAGGACAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTAt                                                                                                                                     >  2:232905/1‑149 (MQ=255)                  gAAAGGACAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCt                                                                                                                       >  3:111198/1‑149 (MQ=255)                    aaGGACAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGACGCGGTATTGACCGGGTGCCCTgg                                                                                                                     <  5:219849/149‑1 (MQ=255)                       gACAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCt                                                                                                                   >  8:262306/1‑148 (MQ=255)                         cAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGCGCGCGGTATTGACCGGGTGCCCTGGCTGAc                                                                                                                <  1:28599/149‑1 (MQ=255)                           gTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCg                                                                                                              <  7:262306/149‑1 (MQ=255)                              gATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTc                                                                                                           >  1:193600/1‑149 (MQ=255)                                       gTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCgg                                                                                                  <  8:5539/149‑1 (MQ=255)                                       gTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTCGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCgg                                                                                                  <  1:15633/149‑1 (MQ=255)                                            aGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGt                                                                                              <  4:111198/148‑1 (MQ=255)                                            aGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTg                                                                                             >  4:199092/1‑149 (MQ=255)                                             gcgcGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTgg                                                                                            >  1:235914/1‑149 (MQ=255)                                             gcgcGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTgg                                                                                            >  4:17285/1‑149 (MQ=255)                                               gcgTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACAGGGGATAAATCAGAATGACCATTGTTGTTGTTGATTTGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTTACCCGTCTTGTGCCGGGAGTggg                                                                                           >  1:267735/1‑149 (MQ=255)                                                cgTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGgcg                                                                                         <  4:39360/149‑1 (MQ=255)                                                 gTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGc                                                                                                     >  4:146532/1‑136 (MQ=255)                                                 gTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGc                                                                                                     <  3:146532/136‑1 (MQ=255)                                                     ggCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTtgt                                                                                                       >  2:197632/1‑130 (MQ=255)                                                     ggCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTtgt                                                                                                       <  1:197632/130‑1 (MQ=255)                                                          agtgCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACAATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTaa                                                                                >  2:245862/2‑148 (MQ=255)                                                                gggAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGc                                                                          >  5:252853/1‑148 (MQ=255)                                                                    gCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCaa                                                                        >  1:155203/1‑146 (MQ=255)                                                                    gCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCaa                                                                        <  2:155203/146‑1 (MQ=255)                                                                      aTGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGACGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCtt                                                                    >  8:226760/1‑148 (MQ=255)                                                                             gggggACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTaaaa                                                                              >  5:227236/1‑131 (MQ=255)                                                                             gggggACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTaaaa                                                                              <  6:227236/131‑1 (MQ=255)                                                                                   cTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCg                                                      <  2:193600/149‑1 (MQ=255)                                                                                   cTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCg                                                      >  1:202321/1‑149 (MQ=255)                                                                                   cTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCg                                                      <  3:17285/149‑1 (MQ=255)                                                                                   cTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCg                                                      <  6:252853/149‑1 (MQ=255)                                                                                    tACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTgg                                                                                                                     >  5:119759/1‑85 (MQ=255)                                                                                    tACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTgg                                                                                                                     <  6:119759/85‑1 (MQ=255)                                                                                    tACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCg                                                      >  5:279042/1‑148 (MQ=255)                                                                                      cGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCg                                                      <  6:279042/146‑1 (MQ=255)                                                                                      cGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCt                                                    <  2:202321/148‑1 (MQ=255)                                                                                       gAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTg                                                   >  2:33548/1‑148 (MQ=255)                                                                                       gAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTg                                                   >  1:47840/1‑148 (MQ=255)                                                                                        aaaaGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGaa                                                 <  1:232905/149‑1 (MQ=255)                                                                                             gCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCtt                                                                                                         >  6:245834/1‑88 (MQ=255)                                                                                             gCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCtt                                                                                                         <  5:245834/88‑1 (MQ=255)                                                                                                        gAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACga                                    >  4:320380/1‑147 (MQ=255)                                                                                                        gAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACga                                    >  4:315356/1‑147 (MQ=255)                                                                                                         aaTACATCAGAATCACCCTTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGagag                                 <  3:199092/149‑1 (MQ=255)                                                                                                              aTCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGaggtagg                            >  4:95405/1‑149 (MQ=255)                                                                                                                        ccATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGa                  <  6:57062/149‑1 (MQ=255)                                                                                                                         cATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCTCTGAAGCAAAGGTGAGGAGAGGTAGGTGATGCAGGAt                 >  8:178045/1‑149 (MQ=255)                                                                                                                         cATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATg                       <  1:277087/143‑1 (MQ=255)                                                                                                                         cATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATg                       >  2:277087/1‑143 (MQ=255)                                                                                                                         cATTGTTGTTGTTGATTGGGCCGGACGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGAt                 <  7:226760/149‑1 (MQ=255)                                                                                                                          aTTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGTCGAGAGGTAGGTGATGCAGGATg                >  7:210759/1‑149 (MQ=255)                                                                                                                              ttgttttttTTTGGGCCGGGCCCGGTATGGACCGGGGGCCTTGGCTGACGGGTCTTGTGCCGGTAGTGGGCGTTATCCTAAAAGGCCATGCTTTCAGCAAGCATTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTaa             <  1:245862/148‑1 (MQ=255)                                                                                                                              ttgttgttgATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAggg                                                                           <  7:222200/86‑1 (MQ=255)                                                                                                                              ttgttgttgATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAggg                                                                           >  8:222200/1‑86 (MQ=255)                                                                                                                                    ttgATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAg       <  3:315356/148‑1 (MQ=255)                                                                                                                                    ttgATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAg       <  3:320380/148‑1 (MQ=255)                                                                                                                                      gATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAg                                           >  4:40305/1‑110 (MQ=255)                                                                                                                                      gATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAg                                           <  3:40305/110‑1 (MQ=255)                                                                                                                                       aTTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTg    <  2:235914/148‑1 (MQ=255)                                                                                                                                       aTTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTg    <  2:267735/148‑1 (MQ=255)                                                                                                                                        ttGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGc  >  5:48517/1‑149 (MQ=255)                                                                                                                                           ggCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGaggt                               >  4:320626/1‑117 (MQ=255)                                                                                                                                           ggCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGaggt                               <  3:320626/117‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                            TCCAGGCGTGGCGCGGCGAAAGGACAGTCGATTCGGGCGTATCAGCGCGTAGGGCGCTGTGCGGGGAGCATGCGGAGGGGGACTACGAAAAGGCGGGCACA‑GGAATACATCAGAATCACCATTGTTGTTGTTGATTGGGCCGGGCGCGGTATTGACCGGGTGCCCTGGCTGACCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGC  >  NC_002947/157734‑158016

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: ‑