Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,960,966	A→T	6.5%	intergenic (‑148/+116)	fliE ← / ← atoC	flagellar hook‑basal body complex protein FliE/two‑component system DNA‑binding transcriptional activator AtoC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,960,966	0	A	T	6.5%	178.0 / 3.2	62	intergenic (‑148/+116)	fliE/atoC	flagellar hook‑basal body complex protein FliE/two‑component system DNA‑binding transcriptional activator AtoC
Reads supporting (aligned to +/- strand):  ref base A (25/33);  new base T (2/2);  total (27/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 = 7.26e-01

GACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGAC  >  NC_002947/4960819‑4961105                                                                                                                                                     |                                                                                                                                            gACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCa                                                                                                                                             >  1:45918/1‑148 (MQ=255)  aCCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCa                                                                                                                                             >  4:250254/1‑147 (MQ=255)   cccTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGAATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaga                                                                                                                                           <  1:265845/149‑1 (MQ=255)   cccTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGTCGCagag                                                                                                                                          >  8:82281/1‑149 (MQ=255)   cccTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCagag                                                                                                                                          >  5:262335/1‑149 (MQ=255)    ccTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGc                                                                                                                                         >  3:13207/1‑149 (MQ=255)           tctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCg                                                                                                                                  >  8:229031/1‑149 (MQ=255)           tctcCGGCGGCCCCATCTTTAGGACAGTCATGCGGCGAATGACCTACCTAAACCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGccc                                                                                                                                    <  6:252714/147‑1 (MQ=39)             tcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGACCGGCCCCGGc                                                                                                                                >  5:77377/1‑149 (MQ=255)               cTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCNTCGCGATAGGGCCGCAGAGCGGCCCCGGCg                                                                                                                               <  1:33813/148‑1 (MQ=255)               cTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGccctatc                                                                                                                                <  3:150722/147‑5 (MQ=255)                   ggCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCtt                                                                                                                          <  6:175994/149‑1 (MQ=255)                     ccGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTTGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTg                                                                                                                        <  1:53463/149‑1 (MQ=255)                     ccGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTg                                                                                                                        <  7:229031/149‑1 (MQ=255)                         aTTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGcccta                                                                                                                                  <  7:167586/135‑3 (MQ=255)                         aTTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGcccta                                                                                                                                  >  8:167586/1‑133 (MQ=255)                          tttttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                   >  7:80075/1‑149 (MQ=255)                          tttttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                   >  1:243439/1‑149 (MQ=255)                          tttttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                   <  1:79529/149‑1 (MQ=255)                          tttttttGACAGTCATGCGGCGAATGACCTACCAAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                   <  2:45918/149‑1 (MQ=255)                             ttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGCTCTTTGCGGCGGAg                                                                                                                >  6:73183/1‑149 (MQ=255)                                 gACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTg                                                                                                             <  5:86077/148‑1 (MQ=255)                                     gTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGAc                                                                                                        >  4:313583/1‑149 (MQ=255)                                         tGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTgg                                                                                                    <  3:250254/149‑1 (MQ=255)                                           cggcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGGCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTggg                                                                                                   <  7:82281/148‑1 (MQ=255)                                               gAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGc                                                                                              <  5:122857/149‑1 (MQ=255)                                                    acctacctAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTcgc                                                                                         <  3:262032/149‑1 (MQ=255)                                                    acctacctAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTcg                                                                                          <  3:313583/148‑1 (MQ=255)                                                    acctacctAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGAACTTTGCGGCGGTGGTGAAGACATGGGGCAGCTAcgc                                                                                         >  4:262032/1‑149 (MQ=38)                                                     cctacctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTcgcg                                                                                        >  7:91255/1‑149 (MQ=255)                                                     cctacctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTcgcg                                                                                        >  6:303058/1‑149 (MQ=255)                                                      ctacctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCgcgc                                                                                       >  5:369974/1‑149 (MQ=255)                                                           taAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCAt                                                                                  >  1:189816/1‑149 (MQ=255)                                                            aaaGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGCCCGGGGGCCGCGTAGCGCCCCAtt                                                                                 >  2:9599/1‑149 (MQ=255)                                                            aaaGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCGTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCAt                                                                                  >  7:302636/1‑148 (MQ=255)                                                                caacaaAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGAGGTTTGCGGGGGCGCTGGAGACCTGGGGCCGCTTGGCGgccccgtggcg                                                                             >  8:303970/1‑142 (MQ=38)                                                                  acaaAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGa                                                                            >  4:303758/1‑148 (MQ=255)                                                                  acaaAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCATGCGTCGCGATAGGGCCGCACAGCGGCCCCGGCGAGCATTGCGGCGGAGCTGAAGACCTGGGGCCCCTTAGCGACCCATTcgacac                                                                           >  4:220312/1‑145 (MQ=38)                                                                    aaaGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGGGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGacac                                                                         <  8:91255/149‑1 (MQ=255)                                                                          gCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGcc                                                                   <  5:258016/149‑1 (MQ=255)                                                                               ccTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTcc                                                              <  3:220312/149‑1 (MQ=255)                                                                                 ttCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCca                                                            >  8:246415/1‑149 (MQ=255)                                                                                    gCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACag                                                         <  2:189816/149‑1 (MQ=255)                                                                                     cAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                                        <  7:267961/149‑1 (MQ=255)                                                                                            tttCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCcgc                                                 <  5:73183/149‑1 (MQ=255)                                                                                                 cAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAg                                            <  7:246415/149‑1 (MQ=255)                                                                                                  agag‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCag                                                                                  <  6:37388/110‑2 (MQ=255)                                                                                                  agag‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCag                                                                                  >  5:37388/1‑109 (MQ=255)                                                                                                       aTCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCa                                       <  8:260383/149‑1 (MQ=255)                                                                                                         cacacaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                                        <  1:157793/130‑1 (MQ=255)                                                                                                           cacaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                                        >  7:230258/1‑128 (MQ=255)                                                                                                           cacaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                                        <  8:230258/128‑1 (MQ=255)                                                                                                               aTTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGAttt                               <  2:243439/149‑1 (MQ=255)                                                                                                               aTTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGAttt                               <  3:303758/149‑1 (MQ=255)                                                                                                                  ccTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGc                             <  5:303058/148‑1 (MQ=255)                                                                                                                      cGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCCCGCCAGTTGCAGATGATTTGCCAGTc                        <  1:9599/149‑1 (MQ=255)                                                                                                                         gAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGc                     >  7:35330/1‑149 (MQ=255)                                                                                                                               ccTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGaa                 >  5:342672/1‑147 (MQ=255)                                                                                                                                 ttGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGAttt                               <  7:328749/131‑1 (MQ=255)                                                                                                                                 ttGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGAttt                               >  8:328749/1‑131 (MQ=255)                                                                                                                                 ttGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGt                         >  1:164445/1‑137 (MQ=255)                                                                                                                                    cGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCt           <  8:35330/148‑1 (MQ=255)                                                                                                                                    cGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTg          <  6:369974/149‑1 (MQ=255)                                                                                                                                      tCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCgcgc                                               >  7:137548/1‑110 (MQ=255)                                                                                                                                      tCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCgcgc                                               <  8:137548/110‑1 (MQ=255)                                                                                                                                         cgATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTc     <  6:342672/149‑1 (MQ=255)                                                                                                                                            cggggCCGCTCTGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCGGGCCTCGAc  >  8:326612/3‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                            GACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAG‑ATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGAC  >  NC_002947/4960819‑4961105

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