Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,654,883	G→T	9.5%	intergenic (+35/‑223)	aceA → / → PP_4118	isocitrate lyase/membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,654,883	0	G	T	9.5%	197.4 / 9.4	63	intergenic (+35/‑223)	aceA/PP_4118	isocitrate lyase/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (25/32);  new base T (3/3);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

ACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGG  >  NC_002947/4654736‑4655030                                                                                                                                                    |                                                                                                                                                    aCCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGTCGCTTCGTCGGTGACTGCGCTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTtgt                                                                                                                                                    <  7:51116/149‑1 (MQ=255) aCCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTtgt                                                                                                                                                    <  4:130576/149‑1 (MQ=255) aCCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTtgt                                                                                                                                                    >  3:130576/1‑149 (MQ=255)      cagGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAg                                                                                                                                               >  7:363177/1‑149 (MQ=255)        ggAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTtgt                                                                                                                                                    >  7:118389/1‑142 (MQ=255)        ggAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTtgt                                                                                                                                                    <  8:118389/142‑1 (MQ=255)        ggAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCt                                                                                                                                             >  3:83032/1‑149 (MQ=255)          aaGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTc                                                                                                                                            <  8:132039/148‑1 (MQ=255)           aGTGGGCACTGGCTACTTCGACGACATGACGACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAATAAGAGCAGTTCCAGTGATAGTGTGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCg                                                                                                                                          <  3:241456/149‑1 (MQ=255)            gTGGCCACTGGCTACATCGACGACATGAACACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCgg                                                                                                                                         <  6:240160/149‑1 (MQ=255)                  aCTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCttt                                                                                                                                   <  4:174885/149‑1 (MQ=255)                         aCTTCCACGACATCACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTa                                                                                                                            <  5:316800/149‑1 (MQ=255)                          cTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACAGACCGCTTCGACCGAAGAAGAGCAGTTCCACTGATAGAGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTTAAAGCTCCGGGCCTTTGTTTTTAc                                                                                                                           >  7:354292/1‑149 (MQ=255)                              gacgacATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCga                                                                                                                       <  4:254380/149‑1 (MQ=255)                                   cATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAg                                                                                                                  >  2:233554/1‑149 (MQ=255)                                    aTGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgg                                                                                                                 <  3:186722/149‑1 (MQ=255)                                    aTGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgg                                                                                                                 >  8:61221/1‑149 (MQ=255)                                      gACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGtt                                                                                                                                <  3:242281/132‑1 (MQ=255)                                      gACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGtt                                                                                                                                >  4:242281/1‑132 (MQ=255)                                           ccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccgggc                                                                                                          <  4:83032/149‑1 (MQ=255)                                                  ccAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGc                                                                                                                                              >  3:357900/1‑106 (MQ=255)                                                  ccAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGc                                                                                                                                              <  4:357900/106‑1 (MQ=255)                                                    aGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGc                                                                                                  >  5:13097/1‑148 (MQ=255)                                                       gcggcgCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAAt                                                                                              <  2:267906/149‑1 (MQ=255)                                                        cggcgCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAAtt                                                                                             <  8:363177/149‑1 (MQ=255)                                                         ggcgCTTCGTCGGTGCCTGCACTGTCCGTTTCGACCTAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACTAATGCCCTGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGACAGGGCCTGGCCTGCTGGCCAAttt                                                                                            <  8:354292/149‑1 (MQ=255)                                                           cgcTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAAtttt                                                                                           >  6:258383/1‑148 (MQ=255)                                                             cTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTAc                                                                                        <  1:233554/149‑1 (MQ=255)                                                                  tcgGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGcc                                                                                   >  4:24263/1‑149 (MQ=255)                                                                     gTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGtt                                                                                <  6:13097/149‑1 (MQ=255)                                                                               tGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTgg                                                                                                   <  1:322700/120‑1 (MQ=255)                                                                               tGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTgg                                                                                                   >  2:322700/1‑120 (MQ=255)                                                                                 aCCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTcc                                                                      <  5:258383/147‑1 (MQ=255)                                                                                         gACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTAcc                                                             <  7:61221/148‑1 (MQ=255)                                                                                          aCCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGt                                                                                                                                 <  3:3746/79‑1 (MQ=255)                                                                                          aCCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGt                                                                                                                                 >  4:3746/1‑79 (MQ=255)                                                                                           ccGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCga                                                                                                                       <  6:389598/88‑1 (MQ=255)                                                                                           ccGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCga                                                                                                                       >  5:389598/1‑88 (MQ=255)                                                                                            cGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTAt                                                         >  2:115710/1‑149 (MQ=255)                                                                                              aagaagAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTc                                                       >  2:168224/1‑149 (MQ=255)                                                                                                gaagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcgc                                                     <  4:216128/149‑1 (MQ=255)                                                                                                 aagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcgcg                                                    <  4:333800/149‑1 (MQ=255)                                                                                                     gCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTa                                                <  3:24263/149‑1 (MQ=255)                                                                                                       aGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCttt                                                                                                                                   >  7:292946/1‑64 (MQ=14)                                                                                                       aGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCttt                                                                                                                                   <  8:292946/64‑1 (MQ=14)                                                                                                        gTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAAc                                             >  7:146138/1‑149 (MQ=255)                                                                                                            cACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCg                                          >  3:267715/1‑148 (MQ=255)                                                                                                             aCTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCa                                        <  1:19824/149‑1 (MQ=255)                                                                                                             aCTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCa                                        >  6:118733/1‑149 (MQ=255)                                                                                                                         ttGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATc                            <  1:168224/149‑1 (MQ=255)                                                                                                                              gTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCTGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                       <  5:234596/149‑1 (MQ=255)                                                                                                                               ttGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                       >  6:358564/1‑148 (MQ=255)                                                                                                                                 gAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                       >  4:177728/1‑146 (MQ=255)                                                                                                                                  aaCAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                   <  3:177728/149‑1 (MQ=255)                                                                                                                                  aaCAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCa                                                                              <  7:328712/90‑1 (MQ=21)                                                                                                                                  aaCAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCa                                                                              >  8:328712/1‑90 (MQ=21)                                                                                                                                   aCAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGac                  <  3:298128/149‑1 (MQ=255)                                                                                                                                   aCAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCgg                                                   >  7:319315/1‑116 (MQ=21)                                                                                                                                   aCAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCgg                                                   <  8:319315/116‑1 (MQ=21)                                                                                                                                    cAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTc                                                                       >  2:335931/1‑95 (MQ=255)                                                                                                                                    cAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTc                                                                       <  1:335931/95‑1 (MQ=255)                                                                                                                                     aaaGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGaca                 >  7:205309/1‑148 (MQ=255)                                                                                                                                      aaGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGttt                                                                        >  2:355205/1‑92 (MQ=255)                                                                                                                                      aaGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGttt                                                                        <  1:355205/92‑1 (MQ=255)                                                                                                                                           ccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCg          <  8:205309/149‑1 (MQ=255)                                                                                                                                              gggCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCAGCTACCGTATTCGCGGGTAAACCCGCAACTACACTGATCGGAAGCAGACACACATCGGc        >  6:98698/1‑148 (MQ=255)                                                                                                                                                 cTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGAtt    <  8:146138/149‑1 (MQ=255)                                                                                                                                                  ttGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTg   >  8:6764/1‑149 (MQ=255)                                                                                                                                                   tgtgAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTgg  >  6:341227/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGG  >  NC_002947/4654736‑4655030

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