Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	379,603	C→T	100%	intergenic (‑32/+3)	frmA ← / ← frmR	alcohol dehydrogenase class III; glutathione‑dependent formaldehyde dehydrogenase/regulator protein that represses frmRAB operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	379,603	0	C	T	90.2%	241.2 / 11.8	89	intergenic (‑32/+3)	frmA/frmR	alcohol dehydrogenase class III; glutathione‑dependent formaldehyde dehydrogenase/regulator protein that represses frmRAB operon
Reads supporting (aligned to +/- strand):  ref base C (5/4);  new base T (42/37);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGG  >  NC_000913/379506‑379699                                                                                                  |                                                                                                 aaCGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAg                                                                                                <  1:378626/100‑1 (MQ=255)    gTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTa                                                                                             >  1:1218250/1‑100 (MQ=255)    gTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTa                                                                                             >  2:1879805/1‑100 (MQ=255)     tCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTAt                                                                                            >  1:721227/1‑100 (MQ=255)      cAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTAtt                                                                                           >  2:1084903/1‑100 (MQ=255)        aTTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTa                                                                                         <  1:1446587/100‑1 (MQ=255)        aTTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTa                                                                                         <  1:2083878/100‑1 (MQ=255)         tttCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTaa                                                                                        >  1:853889/1‑100 (MQ=255)          ttCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAg                                                                                       >  1:392381/1‑100 (MQ=255)            cAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGAt                                                                                     <  1:961707/100‑1 (MQ=255)            cAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGAt                                                                                     <  2:1841139/100‑1 (MQ=255)             aaCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATa                                                                                    <  2:1088696/100‑1 (MQ=255)              aCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGAttg                                                                                   <  1:1765412/100‑3 (MQ=255)              aCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAg                                                                                   <  1:557257/100‑1 (MQ=255)              aCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAg                                                                                   >  1:2119760/1‑100 (MQ=255)               cGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGCTTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAgg                                                                                  >  2:1928130/1‑100 (MQ=255)                 aTTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCa                                                                                <  2:1928219/100‑1 (MQ=255)                 aTTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCa                                                                                >  1:1774943/1‑100 (MQ=255)                    tCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGa                                                                             >  1:1983931/1‑100 (MQ=255)                    tCCAGCGGTTTACCGGGACCAAATGCAACGTCAGCACGTGATTTCATCTCTCGCTCTTCTTCAATATGGTAATAGATTTAGCTATTTAAGATAgg                                                                                  <  2:298131/95‑1 (MQ=255)                    tCCAGCGGATTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGATATTTAAGATAGGCACGa                                                                             >  1:2001803/1‑100 (MQ=255)                      cAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAAc                                                                           >  1:1997863/1‑100 (MQ=255)                      cAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAAc                                                                           >  2:2324869/1‑100 (MQ=255)                      cAGCGCTTTACCGGGAGCAAATGCCACAGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTAGTTAAGATAGGCACGAAc                                                                           <  2:640615/100‑1 (MQ=255)                       aGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAAcc                                                                          <  1:474903/100‑1 (MQ=255)                       aGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGAGTTCATCTCTAGCTATTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAAcc                                                                          <  2:1004208/100‑1 (MQ=255)                         cGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAg                                                                        <  1:1066026/100‑1 (MQ=255)                         cGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAg                                                                        >  1:2432376/1‑100 (MQ=255)                            tttACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTc                                                                     >  1:1301166/1‑100 (MQ=255)                              tagcGGGAGCAAATGCAACGGCAGTACGTGATTTCATCGCTCGCTCTGCCTCAAGATGGTAATAGAGTTAGCTATTTAAGGTTGGCACGAACCAGTTCa                                                                    <  2:746545/96‑1 (MQ=255)                              tACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCaa                                                                   <  2:2029654/100‑1 (MQ=255)                                ccGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTACTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATa                                                                 <  2:2277921/100‑1 (MQ=255)                                  gggAGCAAATGCAACGGCAGCCCGTGATTTCATCTCACGCTCTTCCACAATCTGGCCCTAGATCTAGACATTTACGTTAGCCCCGAACCAGTTCaaca                                                                 >  2:1618248/1‑96 (MQ=255)                                  gggAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATact                                                               >  2:1574611/1‑98 (MQ=255)                                   ggAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAgtg                                                              <  1:2145922/100‑1 (MQ=255)                                     aGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtc                                                            >  1:1192232/1‑100 (MQ=255)                                      gCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcg                                                           >  2:1020545/1‑100 (MQ=255)                                      gCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcg                                                           <  2:1837026/100‑1 (MQ=255)                                      gCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcg                                                           <  2:472867/100‑1 (MQ=255)                                        aaaTGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcgtc                                                         >  2:467833/1‑100 (MQ=255)                                           tGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAAc                                                      <  2:721222/100‑1 (MQ=255)                                             cAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACgg                                                    <  2:1702462/100‑1 (MQ=255)                                                 ggCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTg                                                >  2:122071/1‑100 (MQ=255)                                                  gcagcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTgg                                               <  2:780456/100‑1 (MQ=255)                                                     gcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTg                                            <  1:2320421/100‑1 (MQ=255)                                                        cgtactctCATTTCTCGCTCTTCCTAAAAAGGGAAATCGAGTACGGTGTCTAAGATAGGCACGAAACGATTCAATCGTGTGGTAAACGGATTGGCTGAc                                          >  2:1067722/8‑99 (MQ=255)                                                        cGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCGCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGAc                                          <  2:1060937/99‑1 (MQ=255)                                                        cGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACt                                         >  2:265198/1‑100 (MQ=255)                                                         gTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACtt                                        >  2:100022/1‑100 (MQ=255)                                                          tGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTc                                       >  2:2070021/1‑100 (MQ=255)                                                            aTTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTcgc                                     >  1:1872614/1‑100 (MQ=255)                                                            aTTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGAATTGGCTAACTTcg                                      >  2:964577/1‑99 (MQ=255)                                                             tttCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTcgcg                                    >  2:859863/1‑100 (MQ=255)                                                             tttCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTcgcg                                    >  2:402858/1‑100 (MQ=255)                                                              ttCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAATTCAATAGTGTCGTCAACGTATCGGCTGACTTcgcg                                    <  2:430170/99‑1 (MQ=255)                                                               tCATCTCTCGCTCTTCCTCAATATGGTATTAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGc                                  <  2:438239/100‑1 (MQ=255)                                                               tCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGCTTGGCTGACTTCGCGGc                                  >  1:486287/1‑100 (MQ=255)                                                               tCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTAGCTGACTTCGCGGc                                  <  2:2573008/100‑1 (MQ=255)                                                                cATCTCTCGCTCTTCCTCAATATGGTAATATATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGc                                  <  1:647830/99‑1 (MQ=255)                                                                 acctctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGc                                  <  2:871033/96‑1 (MQ=255)                                                                  tctctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGt                               <  2:1922682/100‑1 (MQ=255)                                                                     ctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGc                            >  1:1219630/1‑100 (MQ=255)                                                                     ctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGc                            <  2:1071321/100‑1 (MQ=255)                                                                       cGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAg                          >  2:1134695/1‑100 (MQ=255)                                                                       cGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAg                          >  2:2150222/1‑100 (MQ=255)                                                                       cGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAg                          >  2:1999483/1‑100 (MQ=255)                                                                       cGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTCTAGCAg                          >  2:2246383/1‑100 (MQ=255)                                                                        gCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAgn                         <  1:2441770/100‑2 (MQ=255)                                                                          tctTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCa                       >  1:2450100/1‑100 (MQ=255)                                                                             tCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCAttt                    <  1:632145/100‑1 (MQ=255)                                                                               cTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCg                  >  2:841462/1‑100 (MQ=255)                                                                               cTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCg                  >  1:924323/1‑100 (MQ=255)                                                                                tCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCgg                 >  1:178069/1‑100 (MQ=255)                                                                                  aaTATGGTAATAGATTCAGCCATTTCAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTc               <  2:1635561/100‑1 (MQ=255)                                                                                   atatGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCa              >  2:432437/1‑100 (MQ=255)                                                                                      tGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAAc           <  1:2316420/100‑1 (MQ=255)                                                                                        gTAATAGATTTAGCTATTTAAGATAGGCACGAACCGGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCaa             <  2:1367976/96‑1 (MQ=255)                                                                                         tAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGtt        >  1:840959/1‑100 (MQ=255)                                                                                          aaTAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGttt       >  1:686429/1‑100 (MQ=255)                                                                                           aTAGATTTAGCTATTTAAGATAGGCACGAACCAGTTAAATAGTGTCGTCCACGGCTTGGCTGACTTCGCGGATGTCGCAGTCATTTCGGTCAAACGTTTc      >  2:760606/1‑100 (MQ=255)                                                                                            tAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTcc     >  2:161468/1‑100 (MQ=255)                                                                                            tAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTcc     >  1:2162282/1‑100 (MQ=255)                                                                                            tAGATTTAGCTATTTAAGATAGGCACCAACCAGTTCAATAGTGTCGTGAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTcc     <  2:899693/100‑1 (MQ=255)                                                                                              gATTTAGCTATTTAAGATAGGCACGAACCCGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCg   <  2:2007144/100‑1 (MQ=255)                                                                                              gATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCg   >  1:435497/1‑100 (MQ=255)                                                                                              gATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCg   >  1:150829/1‑100 (MQ=255)                                                                                              gATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCg   >  1:778310/1‑100 (MQ=255)                                                                                              gATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCg   <  2:2007188/100‑1 (MQ=255)                                                                                               aTTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCgg  <  1:858242/100‑1 (MQ=255)                                                                                                  |                                                                                                 AACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGG  >  NC_000913/379506‑379699

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑