Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,399,148	Δ1 bp	6.6%	coding (981/1302 nt)	tolB →	Tol biopolymer transport system protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,399,148	0	A	.	6.6%	188.0 / 7.2	61	coding (981/1302 nt)	tolB	Tol biopolymer transport system protein
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base . (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCCAGATCAGCCGTGTTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGACTTGCAACGCGGAAGTGTAAAGATTCTGTC  >  NC_002947/1399000‑1399288                                                                                                                                                      |                                                                                                                                              gcCAGATCAGCCGTGTTACTGCAGGCCCGGGGATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑Ta                                                                                                                                               <  4:220936/149‑1 (MQ=255)   cAGATCAGCCGTGTTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAAc                                                                                                                                             <  3:163870/149‑1 (MQ=255)   cAGATCAGCCGTGTTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACAGTGGCGGCAAACCGCAGAACTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAAc                                                                                                                                             >  8:306620/1‑149 (MQ=255)         aGCCGTGTTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGt                                                                                                                                       <  4:260168/149‑1 (MQ=255)         aGCCGTGTTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGt                                                                                                                                       <  6:220134/149‑1 (MQ=255)          gCCGTGTTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCg                                                                                                                                        >  3:360931/1‑147 (MQ=255)          gCCGTGTTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTa                                                                                                                                      <  5:149481/149‑1 (MQ=255)           ccGTGTTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAg                                                                                                                                     <  3:185923/149‑1 (MQ=255)           ccGTGTTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAg                                                                                                                                     >  3:122402/1‑149 (MQ=255)             gtgtTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAgg                                                                                                                                    >  1:87106/1‑148 (MQ=255)             gtgtTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGt                                                                                                                                   <  2:126633/149‑1 (MQ=255)                ttACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAAc                                                                                                                                >  5:172155/1‑149 (MQ=255)                      cAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACATCACTTACGACAGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTTCCTAGCGCG‑TAAATTTCGTAGGTAACTACa                                                                                                                            >  5:95400/1‑147 (MQ=255)                               gTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑T‑ACGTTCGTAGGTAACTACAACGCCAACCCg                                                                                                                 >  6:143359/1‑148 (MQ=255)                               gTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑T‑ACGTTCGTAGGTAACTACAACGCCAACCCg                                                                                                                 >  6:149787/1‑148 (MQ=255)                                  tCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGaaa                                                                                                              <  3:95577/149‑1 (MQ=255)                                   cAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAAc                                                                                                             >  3:253219/1‑149 (MQ=255)                                             ccGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGa                                                                                                   <  3:155151/149‑1 (MQ=255)                                                   tGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGaa                                                                                               >  8:368876/1‑147 (MQ=255)                                                     gggTAAAGATGGCAATACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGaaaa                                                                                             <  7:297244/147‑1 (MQ=255)                                                     gggTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGTCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGaaaa                                                                                             <  4:63070/147‑1 (MQ=255)                                                          aaGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTg                                                                                      >  6:214140/1‑149 (MQ=255)                                                           aGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGGGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTgg                                                                                     <  1:404063/149‑1 (MQ=255)                                                             aTGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGt                                                                                    >  2:300370/1‑148 (MQ=255)                                                                gCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAAt                                                                                 <  5:214140/148‑1 (MQ=255)                                                                 cAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑a                                                                              <  6:172155/149‑1 (MQ=255)                                                                    cacCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑Att                                                                            >  7:285111/1‑148 (MQ=255)                                                                      cccTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑T‑ACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATc                                                                        <  5:143359/149‑1 (MQ=255)                                                                      cccTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑T‑ACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATc                                                                        <  5:149787/149‑1 (MQ=255)                                                                       ccTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATc                                                                        >  7:85943/1‑149 (MQ=255)                                                                        ctctACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCT‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCg                                                                       <  6:95400/149‑1 (MQ=255)                                                                         tctACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCCACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGc                                                                      >  5:35797/1‑149 (MQ=255)                                                                         tctACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGc                                                                      >  7:267069/1‑149 (MQ=255)                                                                                cacttcCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACa                                                                 <  4:360931/147‑1 (MQ=255)                                                                                 acttcCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCACCCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCACCCggg                                                              >  6:111137/1‑149 (MQ=255)                                                                                 acttcCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAggg                                                              <  8:85943/149‑1 (MQ=255)                                                                                  cttcCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGc                                                             >  7:246268/1‑149 (MQ=255)                                                                                      cGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCCACAGGGCTTTa                                                         >  4:110045/1‑149 (MQ=255)                                                                                       gACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑a                                                                              <  1:192872/127‑1 (MQ=255)                                                                                       gACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑a                                                                              >  2:192872/1‑127 (MQ=255)                                                                                       gACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTAc                                                        <  8:285111/149‑1 (MQ=255)                                                                                        aCCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTAcc                                                       <  3:91384/149‑1 (MQ=255)                                                                                             ggcggcAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACtt                                                  >  1:366372/1‑149 (MQ=255)                                                                                                ggcAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACGACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCa                                                >  3:44310/1‑148 (MQ=255)                                                                                                   ataCCGCAGATCTATAAACAGTCGGTCATTGGTGGTTGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGt                                            <  7:306620/147‑1 (MQ=255)                                                                                                   aaaCCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCGTTAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAg                                             >  1:299089/1‑149 (MQ=255)                                                                                                       cGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTggcg                                        <  4:122402/149‑1 (MQ=255)                                                                                                          aGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGca                                     <  2:335401/149‑1 (MQ=255)                                                                                                              cTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAg                                  >  6:374360/1‑148 (MQ=255)                                                                                                               tataAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGa                                >  8:211438/1‑149 (MQ=255)                                                                                                                ataAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGAc                               >  2:245099/1‑149 (MQ=255)                                                                                                                 taAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAg                                  <  2:299089/145‑1 (MQ=255)                                                                                                                   aaCAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGACtt                             <  8:267069/148‑1 (MQ=255)                                                                                                                      aGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATGTATTATGCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGACTTGc                           >  1:335384/1‑148 (MQ=255)                                                                                                                          ggTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGACTTGCAAcgcg                     >  6:145153/1‑149 (MQ=255)                                                                                                                          ggTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGACTTGCAAcgc                      <  7:211438/148‑1 (MQ=255)                                                                                                                           gTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGTGGACGAAAAGACCCTg                                                                                      <  7:384367/84‑1 (MQ=255)                                                                                                                           gTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGTGGACGAAAAGACCCTg                                                                                      >  8:384367/1‑84 (MQ=255)                                                                                                                              tttgggtgtggtgCCGTTCGTG‑TTTCGTTTTTATGTACCTACAACGCCAACCCTAATCTTTCGGCGGACGAAAAGACCCTGGTTATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGACTTGCAACGCGGaa                  <  2:366372/141‑1 (MQ=255)                                                                                                                                tggtggtggtgCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGACTTGCAACGCGGAAgtg               >  4:309299/1‑149 (MQ=255)                                                                                                                                            cGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGACTTGCAACGCGGAAGTGTAAAGATtctg    >  7:191040/1‑148 (MQ=255)                                                                                                                                             gAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGACTTGCAACGCGGAAGTGTAAAGATtctgtc  <  8:246268/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                              GCCAGATCAGCCGTGTTACTGCAGGCCCGGGTATCAATACCGAGCCGTTCTGGGGTAAAGATGGCAACACCCTCTACTTCACTTCCGACCGTGGCGGCAAACCGCAGATCTATAAACAGTCGGTCAGTGGTGGTGGTGCCGAGCGCG‑TAACGTTCGTAGGTAACTACAACGCCAACCCGAAACTTTCGGCGGACGAAAAGACCCTGGTAATG‑ATTCATCGCCAACAGGGCTTTACCAACTTCAAAGTGGCGGCACAGGACTTGCAACGCGGAAGTGTAAAGATTCTGTC  >  NC_002947/1399000‑1399288

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