Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	296,370	C→A	5.5%	A20S (GCG→TCG)	PP_0242 ←	TetR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	296,370	0	C	A	5.5%	294.4 / 4.8	90	A20S (GCG→TCG)	PP_0242	TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (37/48);  new base A (2/3);  total (39/51)
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.16e-01

GCAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGC  >  NC_002947/296222‑296514                                                                                                                                                       |                                                                                                                                                 gCAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGACGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGcc                                                                                                                                                 >  2:2323320/1‑150 (MQ=255)  cAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCt                                                                                                                                                <  2:584917/150‑1 (MQ=255)  cAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCt                                                                                                                                                <  2:168999/150‑1 (MQ=255)           ccTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCAACCCTTGTTCGGCCAACAGCTCACGCGCAGAC‑‑GCCAGGATCGCCTCGATcctgg                                                                                                                                       >  1:2299054/1‑150 (MQ=255)            cTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggc                                                                                                                                      >  1:2282087/1‑150 (MQ=255)            cTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGAGGCGATcctggc                                                                                                                                      <  1:1342979/150‑1 (MQ=255)             tCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAAGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCGAAGAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggcc                                                                                                                                     >  1:633384/1‑150 (MQ=255)                aGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggcctgg                                                                                                                                  >  2:1472030/1‑150 (MQ=255)                aGCAAGGCCGGCACGCTGACGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggcctgg                                                                                                                                  <  1:2323320/150‑1 (MQ=255)                     ggCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCCAGGCCAGGATCGAGGCGATcctggcggctgcg                                                                                                                               <  2:99505/150‑8 (MQ=38)                        cGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGgc                                                                                                                          <  2:184876/150‑1 (MQ=255)                         ggCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGgcg                                                                                                                         <  2:51679/150‑1 (MQ=255)                                     aaaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTgcgc                                                                                                             >  1:2395343/1‑150 (MQ=255)                                     aaaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTgcgc                                                                                                             <  1:100852/150‑1 (MQ=255)                                      aaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGcgcg                                                                                                            <  2:1222833/150‑1 (MQ=255)                                      aaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGcgcg                                                                                                            <  1:952632/150‑1 (MQ=255)                                       aaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                           >  2:766293/1‑150 (MQ=255)                                        aaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCggg                                                                                                          <  1:766293/150‑1 (MQ=255)                                        aaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCggg                                                                                                          <  1:2096317/150‑1 (MQ=255)                                            gATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtgg                                                                                                      <  2:2282087/150‑1 (MQ=255)                                             aTACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggt                                                                                                     <  2:2395343/150‑1 (MQ=255)                                             aTACACCGACGAAGGCGGAATCTGCGCCCGCTCAGCCACGCTGTAGCTATAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggt                                                                                                     >  1:792276/1‑150 (MQ=255)                                              tACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggtg                                                                                                    <  1:591501/150‑1 (MQ=255)                                                 acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATGCTGGGCTGACTGTTGTCTCGTGGTTTGCGCGGGgtggtggtg                                                                                                 >  2:998468/1‑150 (MQ=255)                                                 acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                           <  1:1563278/140‑1 (MQ=255)                                                 acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                           >  2:1563278/1‑140 (MQ=255)                                                 acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGgtggtggtg                                                                                                 >  2:1411737/1‑150 (MQ=255)                                                   cgacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCg                                                                                               <  2:1083112/150‑1 (MQ=255)                                                    gacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCg                                                                                               >  2:516477/1‑149 (MQ=255)                                                    gacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCg                                                                                               <  1:516477/149‑1 (MQ=255)                                                    gacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCGGTAGATAGAACGTGCTGCCACCCCTTATTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGACCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCgg                                                                                              >  1:425939/1‑150 (MQ=255)                                                          ggCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                        >  1:741419/1‑150 (MQ=255)                                                          ggCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                        <  1:868170/150‑1 (MQ=255)                                                           gCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatca                                                                                       >  2:1941509/1‑150 (MQ=255)                                                             ggAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAc                                                                                     <  2:393655/150‑1 (MQ=255)                                                               aagcTGCGCCCGCTCGGCCACGCTGTAGATAGAGCCTCCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGACCCTGGCCTGGCTGCTGGCGCGTGGTGTGCGCGGGGCGGTGGTGCGGTTCATCAActg                                                                                   <  2:633384/147‑1 (MQ=255)                                                                aTCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAActgc                                                                                  <  1:1533419/150‑1 (MQ=255)                                                                     cgcCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGCCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGACCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGCGGTGGTGGTGCGGTTCATCAACTGCTGAtt                                                                             <  2:425939/150‑1 (MQ=255)                                                                         cGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCa                                                                         <  1:956996/150‑1 (MQ=255)                                                                          gCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGGTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAg                                                                        <  1:924456/150‑1 (MQ=255)                                                                            tCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                        >  1:821716/1‑132 (MQ=255)                                                                            tCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                        <  2:821716/132‑1 (MQ=255)                                                                                ccACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCGATTTCTGATTGCCAGGATGCt                                                                  >  2:2361293/1‑150 (MQ=255)                                                                                ccACGCTGGAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCATCTGCTGATTGCCAGGATGCt                                                                  <  2:754605/150‑1 (MQ=255)                                                                                 cACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTg                                                                 <  2:608153/150‑1 (MQ=255)                                                                                   cGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGc                                                               <  1:1472030/150‑1 (MQ=255)                                                                                      tGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGtt                                                                                            >  2:1211865/1‑118 (MQ=255)                                                                                      tGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGtt                                                                                            <  1:1211865/118‑1 (MQ=255)                                                                                       gTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGAGGCGATcctggc                                                                                                                                      <  2:1125176/75‑1 (MQ=38)                                                                                       gTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGAGGCGATcctggc                                                                                                                                      >  1:1125176/1‑75 (MQ=38)                                                                                         agatagaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTa                                                         <  2:1504465/150‑1 (MQ=255)                                                                                           atagaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcat                                                                                         >  2:736724/1‑116 (MQ=255)                                                                                           atagaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcat                                                                                         <  1:736724/116‑1 (MQ=255)                                                                                            tagaAAGTGCTGCCACGCTTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGgc                                                      <  2:792276/150‑1 (MQ=255)                                                                                              gaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAct                                                                                    >  1:904468/1‑118 (MQ=255)                                                                                              gaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAct                                                                                    <  2:904468/118‑1 (MQ=255)                                                                                               aaaGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGcc                                                   <  2:1778904/150‑1 (MQ=255)                                                                                                  gTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCcgac                                                <  1:2059781/150‑1 (MQ=255)                                                                                                      tgcCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAAc                                            >  2:1498485/1‑150 (MQ=255)                                                                                                          aCCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGt                                        >  2:1562816/1‑150 (MQ=255)                                                                                                            cccTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCg                                      >  2:1887518/1‑150 (MQ=255)                                                                                                              cTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggtg                                                                                                    >  2:1197414/1‑86 (MQ=255)                                                                                                              cTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggtg                                                                                                    <  1:1197414/86‑1 (MQ=255)                                                                                                              cTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGaa                                    <  2:957091/150‑1 (MQ=255)                                                                                                                 gTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAc                                                                                     >  2:943738/1‑98 (MQ=255)                                                                                                                 gTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAc                                                                                     <  1:943738/98‑1 (MQ=255)                                                                                                                  ttCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCa                                >  1:2106263/1‑150 (MQ=255)                                                                                                                     ggCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGa                                     <  2:1580699/142‑1 (MQ=255)                                                                                                                     ggCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGa                                     >  1:1580699/1‑142 (MQ=255)                                                                                                                     ggCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCg                             >  2:278930/1‑150 (MQ=255)                                                                                                                           cAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCg                       >  2:2085567/1‑150 (MQ=255)                                                                                                                           cAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCg                       <  2:741419/150‑1 (MQ=255)                                                                                                                             gCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCGGTGACTTCgac                     <  2:471893/150‑3 (MQ=255)                                                                                                                             gCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGtt                     <  2:2106263/150‑1 (MQ=255)                                                                                                                              cTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGttt                    <  1:1411737/150‑1 (MQ=255)                                                                                                                               tCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGtttt                   >  1:1182031/1‑150 (MQ=255)                                                                                                                               tCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGGTCATCAACTGCTGATTGCCATGAGGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGtttt                   >  1:1741818/1‑150 (MQ=255)                                                                                                                                 aCGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGttt                    <  2:2105572/147‑1 (MQ=255)                                                                                                                                 aCGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCCGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGttt                    >  1:2105572/1‑147 (MQ=255)                                                                                                                                   gcgacgcc‑‑gctAGGATCGCTTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTc               >  2:600258/6‑150 (MQ=255)                                                                                                                                       aGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAg                                          <  2:1235782/119‑1 (MQ=255)                                                                                                                                       aGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAg                                          >  1:1235782/1‑119 (MQ=255)                                                                                                                                       aGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGc           >  2:2015883/1‑150 (MQ=255)                                                                                                                                         cc‑‑gccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTc         <  1:998468/150‑1 (MQ=255)                                                                                                                                             gccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGaa                                             >  2:1570321/1‑112 (MQ=255)                                                                                                                                             gccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGaa                                             <  1:1570321/112‑1 (MQ=255)                                                                                                                                             gccAGGATCGAGGCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCgg       >  1:183914/1‑150 (MQ=255)                                                                                                                                              ccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGa      <  1:1941509/150‑1 (MQ=255)                                                                                                                                                  gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                           <  2:1053330/45‑1 (MQ=255)                                                                                                                                                  gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                           >  1:1053330/1‑45 (MQ=255)                                                                                                                                                  gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGc  >  1:1042964/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                 GCAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGC  >  NC_002947/296222‑296514

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