Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,750,363	G→A	8.0%	A419V (GCC→GTC)	thiI ←	tRNA sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,750,363	0	G	A	8.0%	203.7 / 7.5	75	A419V (GCC→GTC)	thiI	tRNA sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (38/31);  new base A (3/3);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

GTGCAAACGGCTCATCACACCCTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGC  >  NC_002947/5750221‑5750500                                                                                                                                               |                                                                                                                                          gtgcAAACGGCTCATCACACCCTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCggg                                                                                                                                    <  2:938146/150‑1 (MQ=255)     aaaCGGCTCATCACACCCTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTgg                                                                                                                                >  1:1162153/1‑150 (MQ=255)     aaaCGGCTCATCACACCCTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTgg                                                                                                                                <  1:907673/150‑1 (MQ=255)         ggCTCATCACACCCTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGaa                                                                                                                            >  1:11852/1‑150 (MQ=255)         ggCTCATCACACCCTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGaa                                                                                                                            >  1:186711/1‑150 (MQ=255)         ggCTCATCACACCCTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGaa                                                                                                                            >  2:317013/1‑150 (MQ=255)              atcaCACCCTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCg                                                                                                                       <  1:121547/150‑1 (MQ=255)                cacaCCCTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGAt                                                                                                                     <  1:317013/150‑1 (MQ=255)                  cacCCTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGa                                                                                                                   <  2:1162153/150‑1 (MQ=255)                     ccTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGa                                                                                                                >  2:861604/1‑150 (MQ=255)                       ttGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATC‑GGGTGGCGAATGTCGATGACGATct                                                                                                             >  1:595147/1‑150 (MQ=255)                             cAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATctggcct                                                                                                        >  1:971868/1‑150 (MQ=255)                              aaTACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATC‑GGGTGGCGAATGTCGATGACGATctggcctgg                                                                                                      <  2:595147/150‑1 (MQ=255)                               aTACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATctggcctgg                                                                                                      <  1:861604/150‑1 (MQ=255)                                     cAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATctggcc                                                                                                         <  2:793517/141‑1 (MQ=255)                                     cAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATctggcc                                                                                                         >  1:793517/1‑141 (MQ=255)                                     cAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAgcgc                                                                                                >  2:278895/1‑150 (MQ=255)                                        gTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTc                                                                                             >  1:1402869/1‑150 (MQ=255)                                           cTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGc                                                                                          <  1:930752/150‑1 (MQ=255)                                                gtgTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACtt                                                                                     >  2:420810/3‑150 (MQ=255)                                                  gTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCgg                                                                                           <  1:1332101/142‑1 (MQ=255)                                                  gTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCgg                                                                                           >  2:1332101/1‑142 (MQ=255)                                                   tGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCt                                                                                  <  2:186711/150‑1 (MQ=255)                                                    ggggTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTc                                                                                 <  2:971868/150‑1 (MQ=255)                                                    ggggTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTc                                                                                 <  1:278895/150‑1 (MQ=255)                                                       gTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATct                                                                                                             >  1:605715/1‑119 (MQ=255)                                                       gTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATct                                                                                                             <  2:605715/119‑1 (MQ=255)                                                       gTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGAt                                                                              >  1:367848/1‑150 (MQ=255)                                                            ggTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGa                                                                         >  1:875345/1‑150 (MQ=255)                                                            ggTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGa                                                                         <  1:861132/150‑1 (MQ=255)                                                                 ttGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCAc                                                                                       >  1:470696/1‑131 (MQ=255)                                                                 ttGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCAc                                                                                       <  2:470696/131‑1 (MQ=255)                                                                 ttGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTcc                                                                    >  2:276547/1‑150 (MQ=255)                                                                   gAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCtt                                                                  >  2:912107/1‑150 (MQ=255)                                                                     tcTTGCTGTTGATGGCGTAGACCGGACTGGCGTAAATTTAGGAATTTACTATCCCGAGAGGTTCTTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGCTCTGGCCTGCCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGc                                                                <  2:667204/149‑1 (MQ=255)                                                                       ttGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGACAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGcca                                                              >  1:667204/1‑149 (MQ=255)                                                                          cTGTTGATGGCGTAGAACGGCATTGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGc                                                           <  2:110196/150‑1 (MQ=255)                                                                             ttGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCa                                                        <  1:276547/150‑1 (MQ=255)                                                                                     gTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATcac                                                >  1:166890/1‑150 (MQ=255)                                                                                       aGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACAt                                              <  2:1402869/150‑1 (MQ=255)                                                                                       aGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACAt                                              >  2:116371/1‑150 (MQ=255)                                                                                                      ggACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGa                               >  2:1130726/1‑150 (MQ=255)                                                                                                        aCTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGAtt                                                                             >  2:703756/1‑102 (MQ=255)                                                                                                        aCTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGAtt                                                                             <  1:703756/102‑1 (MQ=255)                                                                                                          ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTc                                                                                                                        >  1:1013393/1‑57 (MQ=255)                                                                                                          ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTc                                                                                                                        <  2:1013393/57‑1 (MQ=255)                                                                                                          ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACtt                           <  1:1130726/150‑1 (MQ=255)                                                                                                          ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACtt                           >  1:822740/1‑150 (MQ=255)                                                                                                          ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACtt                           >  1:1285789/1‑150 (MQ=255)                                                                                                          ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGACATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGcc                                                               >  1:712640/1‑114 (MQ=255)                                                                                                          ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGACATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGcc                                                               <  2:712640/114‑1 (MQ=255)                                                                                                          ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGACATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGcc                                                               <  2:712411/114‑1 (MQ=255)                                                                                                          ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGACATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGcc                                                               >  1:1426216/1‑114 (MQ=255)                                                                                                          ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGACATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGcc                                                               <  2:1426216/114‑1 (MQ=255)                                                                                                          ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGACATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGcc                                                               >  1:712411/1‑114 (MQ=255)                                                                                                                  cTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCTTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCt                   >  2:85394/1‑150 (MQ=255)                                                                                                                     tcttcACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCaa                <  1:849914/146‑1 (MQ=255)                                                                                                                      cAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAag               >  1:1398758/1‑150 (MQ=255)                                                                                                                        gCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAagag             >  1:726059/1‑150 (MQ=255)                                                                                                                           ccAGAGGTTCGTCTTCCTGGGCATGGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAgcgc          >  2:85896/1‑150 (MQ=255)                                                                                                                           ccAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAgcgc          >  1:1193041/1‑150 (MQ=255)                                                                                                                           ccAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAgcgc          >  1:271217/1‑150 (MQ=255)                                                                                                                             agagGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCTTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgc        <  1:85394/150‑1 (MQ=255)                                                                                                                               agGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTc      >  2:879315/1‑150 (MQ=255)                                                                                                                                 gTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCaa                >  2:284907/1‑138 (MQ=255)                                                                                                                                 gTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCaa                <  1:284907/138‑1 (MQ=255)                                                                                                                                  ttCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAaga              >  2:1195385/1‑139 (MQ=255)                                                                                                                                  ttCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAaga              <  1:1195385/139‑1 (MQ=255)                                                                                                                                   tCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGc  <  1:912107/150‑1 (MQ=255)                                                                                                                                      tCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCa                                                                                        >  2:370502/1‑61 (MQ=255)                                                                                                                                      tCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCa                                                                                        <  1:370502/61‑1 (MQ=255)                                                                                                                                      tCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGAt                                       <  1:1282471/110‑1 (MQ=255)                                                                                                                                      tCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGAt                                       >  2:1282471/1‑110 (MQ=255)                                                                                                                                        ttCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGt                                                                      <  2:961241/77‑1 (MQ=255)                                                                                                                                        ttCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGt                                                                      >  1:961241/1‑77 (MQ=255)                                                                                                                                               |                                                                                                                                          GTGCAAACGGCTCATCACACCCTTGTCGCAATACAGCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGC  >  NC_002947/5750221‑5750500

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