Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,071,673	C→T	6.2%	P31P (CCG→CCA)	rutD ←	putative reactive intermediate detoxifying aminoacrylate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,071,673	0	C	T	6.2%	203.0 / 2.9	65	P31P (CCG→CCA)	rutD	putative reactive intermediate detoxifying aminoacrylate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base T (2/2);  total (35/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 = 8.75e-01

CTCAATCCCTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTC  >  NC_000913/1071526‑1071817                                                                                                                                                    |                                                                                                                                                 cTCACTCCCTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGc                                                                                                                                                  <  5:118630/148‑1 (MQ=255) cTCAATCCCTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCg                                                                                                                                                 <  1:330488/149‑1 (MQ=255)      tCCCTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAAc                                                                                                                                            <  7:391384/149‑1 (MQ=255)        ccTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCa                                                                                                                                          <  2:211197/149‑1 (MQ=255)        ccTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCGGGGCGATACTGTAATCTGCTGCCAGCGTGTCGGGATTATTGCCGGCGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAAc                                                                                                                                            <  1:105552/147‑1 (MQ=255)        ccTGCGGCTACCAGCGCCTGATGCAGTTACGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGGCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGAGACTCCTGCTCCAGCACCGACAGGTGCGGTAAc                                                                                                                                            >  5:74226/1‑147 (MQ=255)                 aCCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGTGGTAACCAGTAACTGcc                                                                                                                                 >  5:365391/1‑149 (MQ=255)                   cAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAAGCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGATGGTCGTAACAGACTACCTGATAGTCATGCTCCAGTACCGCCACTTGCGGTAACCAGTAACTGCCGc                                                                                                                               >  1:106781/1‑149 (MQ=255)                         cTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTAcccc                                                                                                                         >  8:413193/1‑149 (MQ=255)                            aTGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTAccccc                                                                                                                        <  2:20699/147‑1 (MQ=255)                              gCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTAcc                                                                                                                           >  8:199709/1‑142 (MQ=255)                              gCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTAcc                                                                                                                           <  7:199709/142‑1 (MQ=255)                                    ccGCTGCCATCTGTGCGATACTGTAATCTTCTTCCAGCGTGTCGGGATTATTGCGGGTGCCGCGCTGGTCTTAACAGACTACCTGATACTCCTTCTCCAGCACCGCCAGTTGCGGTAAACAGTAACTGCAGCTACCCCCAAGACCCGa                                                                                                               <  6:193722/148‑1 (MQ=255)                                        tGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATc                                                                                                           >  2:327167/1‑148 (MQ=255)                                          ccATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGTGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAac                                                                                                        <  6:365391/149‑1 (MQ=255)                                          ccATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAac                                                                                                        >  4:294339/1‑149 (MQ=255)                                                         tGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGTGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCa                                                                                         >  3:92344/1‑149 (MQ=255)                                                         tGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCa                                                                                         >  2:300753/1‑149 (MQ=255)                                                             aTCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCa                                                                                      <  7:413193/148‑1 (MQ=255)                                                                ttctGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCa                                                                                   >  5:386957/1‑148 (MQ=255)                                                                  ctGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGc                                                                                    <  1:50570/145‑1 (MQ=255)                                                                       aGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAggggg                                                                           >  6:257547/1‑149 (MQ=255)                                                                        gCGTTTCGGGATTATTGACGGTGCCGCTCTGGTCGTAACATACTACCTGATAATCCTGCTCCAGCACCGCCAGTTGCGGTAAATAGTAATTGACTCTAACCCCAAGAAACGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGa                                                                          <  6:386957/149‑1 (MQ=255)                                                                        gCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTACCCAGTAACTGCCGCTACCCCCCAGCCCCGAAATCAACCCCACGACGGGCGCAGCAGCATAAgggg                                                                            >  2:381060/1‑147 (MQ=255)                                                                          gtgtCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAgg                                                                        <  1:327167/149‑1 (MQ=255)                                                                             tCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGt                                                                       >  6:86303/1‑147 (MQ=255)                                                                             tCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGt                                                                       >  7:188250/1‑147 (MQ=255)                                                                             tCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTg                                                                      <  1:300753/148‑1 (MQ=255)                                                                               gggATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGTGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTgaga                                                                   <  4:92344/149‑1 (MQ=255)                                                                                   ttattGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTg                                                                >  4:123034/1‑148 (MQ=255)                                                                                      ttGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGaaa                                                             <  3:294339/148‑1 (MQ=255)                                                                                      ttGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGaaa                                                             <  5:86303/148‑1 (MQ=255)                                                                                        gCCTGTGCCGCGCTGGTCGTAACAGACTACCTGATAGTCCTGCTCCAGTTCAGCGAGTTGAGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATGATCACCAATACGGGCTCATCAGCATAAGGGGGAGGTGAGAGTGaaatt                                                           <  2:106781/148‑3 (MQ=255)                                                                                              gCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATcg                                                    >  1:264486/1‑149 (MQ=255)                                                                                                cgcgCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATcgcg                                                  >  1:342834/1‑149 (MQ=255)                                                                                                     tGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGc                                                <  4:308341/146‑1 (MQ=255)                                                                                                     tGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGc                                                >  3:308341/1‑146 (MQ=255)                                                                                                             aCAGCCTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAACTCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGc                                      <  1:381060/148‑1 (MQ=255)                                                                                                               aGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGa                                    >  4:205435/1‑148 (MQ=255)                                                                                                                   tACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCaa                                                                                                         <  7:153135/75‑1 (MQ=255)                                                                                                                   tACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCaa                                                                                                         >  8:153135/1‑75 (MQ=255)                                                                                                                      cTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGgc                                                                                            <  1:184767/85‑1 (MQ=255)                                                                                                                      cTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGgc                                                                                            >  2:184767/1‑85 (MQ=255)                                                                                                                         aTACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCCCTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAAtt                         >  1:222559/1‑149 (MQ=255)                                                                                                                         aTACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCaa                           <  2:264486/147‑1 (MQ=255)                                                                                                                            cTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATtgtg                      <  8:188250/149‑1 (MQ=255)                                                                                                                             tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGACATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTgg                     >  7:381598/1‑149 (MQ=255)                                                                                                                             tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAggg                                                                             <  5:314875/93‑1 (MQ=255)                                                                                                                             tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAggg                                                                             >  6:314875/1‑93 (MQ=255)                                                                                                                             tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACt                                          <  2:173943/128‑1 (MQ=255)                                                                                                                             tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACt                                          >  1:173943/1‑128 (MQ=255)                                                                                                                             tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTgg                     >  6:294152/1‑149 (MQ=255)                                                                                                                               cTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGc                    >  6:218876/1‑148 (MQ=255)                                                                                                                               cTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACATGGCGATATCCGCAATTGTTGc                    >  6:336385/1‑148 (MQ=255)                                                                                                                                tGCTCCCGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGt                                                                 <  4:10690/102‑1 (MQ=255)                                                                                                                                tGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGa                  >  4:173917/1‑149 (MQ=255)                                                                                                                                 gCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCTTCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGAt                 >  6:203644/1‑149 (MQ=255)                                                                                                                                  cTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGAtt                <  5:218876/149‑1 (MQ=255)                                                                                                                                     cAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGAc                                                                                                                   >  6:234993/1‑47 (MQ=255)                                                                                                                                     cAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGAc                                                                                                                   <  5:234993/47‑1 (MQ=255)                                                                                                                                          ccGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGTTTTCCACCAg        >  1:153242/1‑149 (MQ=255)                                                                                                                                               aGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGAttt               <  4:103251/137‑1 (MQ=255)                                                                                                                                               aGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGAttt               >  3:103251/1‑137 (MQ=255)                                                                                                                                                gTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGt   >  6:150869/1‑148 (MQ=255)                                                                                                                                                gTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTc  <  3:205435/149‑1 (MQ=255)                                                                                                                                                 ttGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGt   >  2:128109/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CTCAATCCCTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTC  >  NC_000913/1071526‑1071817

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