Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,592,024	C→A	21.7%	intergenic (‑89/‑114)	PP_1395 ← / → PP_1396	AraC family transcriptional regulator/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,592,024	0	C	A	21.7%	134.2 / 27.3	65	intergenic (‑89/‑114)	PP_1395/PP_1396	AraC family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base A (9/5);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCTCGGCCAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCTTTC  >  NC_002947/1591875‑1592173                                                                                                                                                      |                                                                                                                                                      gCTCGGCCAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACAGGCTCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGcc                                                                                                                                                       >  1:373578/1‑150 (MQ=255)         aGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCa                                                                                                                                               >  2:211017/1‑150 (MQ=255)            tCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCCGATCCTGCAGCGCGGCAGGa                                                                                                                                            <  2:636805/150‑1 (MQ=255)              cAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATc                                                                                                                                          >  2:338716/1‑150 (MQ=255)               gggtgcTGCAGTTCCAGGGTCGGTTCGGGGGATGTGGCGGTGGTCATGCTGGCTTGCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCa                                                                                                                                         <  2:421611/146‑1 (MQ=38)               aGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCa                                                                                                                                         <  2:1072952/150‑1 (MQ=255)                 gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCAgg                                                                                                                                       <  1:154333/150‑1 (MQ=255)                 gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGACACGGCCGCGCTGCAGGATCAgg                                                                                                                                       <  2:127099/150‑1 (MQ=255)                 cagcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGCTGGCTTTCCTCTTCTCACCGACGCAGCATAGGTTTGTGCAAAGACGCGCGCTAGTCGAAAGCTGCATGTCGATTGCCTGATCCTGCCGCTCTGCAGGTTCAgg                                                                                                                                       <  2:266683/148‑1 (MQ=255)                          ttCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACa                                                                                                                              >  2:521390/1‑150 (MQ=255)                             cAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGc                                                                                                                           >  2:59333/1‑150 (MQ=255)                             cAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGc                                                                                                                           <  2:147664/150‑1 (MQ=255)                             cAGGGTCGGTTCGGGGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAgg                                                                                                                                       >  2:1293190/1‑138 (MQ=255)                             cAGGGTCGGTTCGGGGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAgg                                                                                                                                       <  1:1293190/138‑1 (MQ=255)                                ggTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGt                                                                                                                        >  1:1256922/1‑150 (MQ=255)                                          gtGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTaa                                                                                                              >  1:1183204/1‑150 (MQ=255)                                              aCGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAAcgcg                                                                                                          <  1:521390/150‑1 (MQ=255)                                                      gtggtCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCgg                                                                                                  <  2:1256922/150‑1 (MQ=255)                                                            atgatgGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGAc                                                                                                                               <  2:867356/115‑1 (MQ=255)                                                            atgatgGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGAc                                                                                                                               >  1:867356/1‑115 (MQ=255)                                                               atgGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCa                                                                                                                                         <  2:507612/102‑1 (MQ=37)                                                               atgGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCa                                                                                                                                         >  1:507612/1‑102 (MQ=37)                                                                   cTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAAc                                                                                     >  1:1137532/1‑150 (MQ=255)                                                                      tCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCtt                                                                                  <  1:662853/150‑1 (MQ=255)                                                                        ctctGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTgg                                                                                <  2:289991/150‑1 (MQ=255)                                                                              tgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCTGCGCTGCAGGATCAgg                                                                                                                                       >  2:153632/1‑89 (MQ=255)                                                                              tgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCTGCGCTGCAGGATCAgg                                                                                                                                       <  1:153632/89‑1 (MQ=255)                                                                               gACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGGGGGTGAGCGCCTAACCTTGGCAGCCTg                                                                         <  2:270047/150‑1 (MQ=255)                                                                                      gcagcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCggg                                                                                                 >  2:462414/1‑119 (MQ=255)                                                                                      gcagcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCggg                                                                                                 <  1:462414/119‑1 (MQ=255)                                                                                      gcagcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtc                                                                  >  1:154136/1‑150 (MQ=255)                                                                                         gcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgt                                                               <  1:338716/150‑1 (MQ=255)                                                                                           aTAGGTTTGGGCCAAGACGGGCGTTAGTCCCAAGCCGCATTTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgtcc                                                             <  2:373578/150‑1 (MQ=255)                                                                                           aTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgtcc                                                             <  2:154136/150‑1 (MQ=255)                                                                                           aTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtcc                                                                 <  1:829732/146‑1 (MQ=255)                                                                                           aTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtcc                                                                 >  2:829732/1‑146 (MQ=255)                                                                                               gTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTg                                                                         >  1:480692/1‑134 (MQ=255)                                                                                               gTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTg                                                                         <  2:480692/134‑1 (MQ=255)                                                                                                              ggCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATtg                                          <  2:1137532/150‑1 (MQ=255)                                                                                                                              tGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGa                                                                                              >  2:722039/1‑82 (MQ=255)                                                                                                                              tGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGa                                                                                              <  1:722039/82‑1 (MQ=255)                                                                                                                                 aTGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaacaa                       >  2:956931/1‑150 (MQ=255)                                                                                                                                 aTGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaacaa                       >  1:1313417/1‑150 (MQ=255)                                                                                                                                 aTGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaacaa                       >  1:307168/1‑150 (MQ=255)                                                                                                                                  tgtgGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaacaac                      >  1:503705/1‑150 (MQ=255)                                                                                                                                   gtgGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAAcc                     >  1:1014005/1‑150 (MQ=255)                                                                                                                                           ccTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccc             >  2:76449/1‑150 (MQ=255)                                                                                                                                            cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgtc                                                              <  2:933728/100‑1 (MQ=255)                                                                                                                                            cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgtc                                                              >  1:933728/1‑100 (MQ=255)                                                                                                                                            cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca            >  2:326235/1‑150 (MQ=255)                                                                                                                                            cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca            >  2:892500/1‑150 (MQ=255)                                                                                                                                            cTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGTAAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca            >  2:711507/1‑150 (MQ=255)                                                                                                                                            cTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca            >  2:1394673/1‑150 (MQ=255)                                                                                                                                            cTGATCCTGCAGCGCGGCAGGATCAGGCAATCCACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca            >  2:788060/1‑150 (MQ=255)                                                                                                                                             tGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATTCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcacccac           >  1:824214/1‑150 (MQ=255)                                                                                                                                             tGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtc                                                                  <  2:1376525/95‑1 (MQ=255)                                                                                                                                             tGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtc                                                                  >  1:1376525/1‑95 (MQ=255)                                                                                                                                             tGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGAtt                                           <  2:1109932/118‑1 (MQ=255)                                                                                                                                             tGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGAtt                                           >  1:1109932/1‑118 (MQ=255)                                                                                                                                             tGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCt                                                                          >  2:1396410/1‑87 (MQ=255)                                                                                                                                             tGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCt                                                                          <  1:1396410/87‑1 (MQ=255)                                                                                                                                             tGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcacccac           >  1:966018/1‑150 (MQ=255)                                                                                                                                              gATCCTGCAGCCCGGCAGGATCAGGCAATCCACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcacccacc          >  2:172060/1‑150 (MQ=255)                                                                                                                                               aTCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCt         >  1:415294/1‑150 (MQ=255)                                                                                                                                                    gCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCtt    <  2:503705/150‑1 (MQ=255)                                                                                                                                                    gCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCtt    <  2:1313417/150‑1 (MQ=255)                                                                                                                                                      cgcgCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCAGTCTGCTTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCTTTc  <  2:194407/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GCTCGGCCAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCTTTC  >  NC_002947/1591875‑1592173

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