Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	379,603	C→T	84.9%	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	84.9%	228.0 / 31.1	99	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 (9/6);  new base T (32/52);  total (41/58)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

AACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGG  >  NC_000913/379506‑379699                                                                                                  |                                                                                                 aaCGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAg                                                                                                <  2:1269543/100‑1 (MQ=255)  aCGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGc                                                                                               <  1:1897265/100‑1 (MQ=255)  aCGTCAATTTCAACGATTTACAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGc                                                                                               <  1:1130186/100‑1 (MQ=255)    gTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTa                                                                                             <  1:2321045/100‑1 (MQ=255)    gTCAATTGCAACAGGTTGCAGCAGTTTACCGGGAGCAAACGCAACGGCAGCACGTGATTTCATCGCTAGCTCTGCCTCAATAGGGTAATAGATTTAGc                                                                                               <  1:314897/98‑1 (MQ=255)       aaTTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTAttt                                                                                          >  1:1036538/1‑100 (MQ=255)         tttCACCGATTCCCAGGGGTTTACCGAGAGCAAATTATACCGCAGCACGTTATTTCATTTCTCGCGCTACCTCAATATGGTAATAGATTTAGCTATTTaa                                                                                        <  1:1099160/100‑1 (MQ=255)           tCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGa                                                                                      <  2:1448112/100‑1 (MQ=255)            cAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGAt                                                                                     <  2:606455/100‑1 (MQ=255)            cAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGAt                                                                                     <  1:189072/100‑1 (MQ=255)            cAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGAt                                                                                     <  2:2852166/100‑1 (MQ=255)             aaCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATa                                                                                    >  2:784002/1‑100 (MQ=255)             aaCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATa                                                                                    >  2:415820/1‑100 (MQ=255)              aCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAg                                                                                   <  2:807025/100‑1 (MQ=255)              aCGATGTGCAGCGGTGTACCGGGCGCATATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAg                                                                                   <  1:2813372/100‑1 (MQ=255)                gATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGc                                                                                 <  1:2905758/100‑1 (MQ=255)                gATTTCCAGAGCTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAAAATGGTAATAGATTTAGCTATTTAAGATAGGc                                                                                 <  1:1417590/100‑1 (MQ=255)                 aTTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCa                                                                                >  2:338119/1‑100 (MQ=255)                 aTTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCa                                                                                >  2:2086578/1‑100 (MQ=255)                   ttCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACg                                                                              >  1:1824719/1‑100 (MQ=255)                    tCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTGCAGCTCTCTCTCTTCCTCAATATGGTAATAGATTTAGCTA‑TTAAGATAGGCACg                                                                              <  1:1772341/98‑1 (MQ=255)                     ccAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGaa                                                                            >  1:2898955/1‑100 (MQ=255)                      cAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAAc                                                                           >  2:2265745/1‑100 (MQ=255)                      cAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAAc                                                                           >  2:1951793/1‑100 (MQ=255)                       aGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAAcc                                                                          <  2:2288702/100‑1 (MQ=255)                       aGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAAcc                                                                          <  2:284328/100‑1 (MQ=255)                           gTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGtt                                                                      <  1:2256569/100‑1 (MQ=255)                           gTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGt                                                                       >  2:1469859/1‑99 (MQ=255)                            tttACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTc                                                                     <  2:1243366/100‑1 (MQ=255)                            tttACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTc                                                                     <  1:897741/100‑1 (MQ=255)                            tttACCCGGAGCAAATGCAACGGCAGCACGTTATTTCATCGCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAg                                                                        <  1:1211292/97‑1 (MQ=255)                             ttACCGGGAGCAAATGCAACTGCACCATGTGATTTCAGCTCGCGCGCTTCCTTAATATGGTAATAGAGTTAGCCATTTAAGATAGGCACTAACCAGTCCa                                                                    <  1:795343/100‑1 (MQ=255)                              tACCGGGAGCAAATGCAACGGCAGAACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCaa                                                                   <  1:1159027/100‑1 (MQ=255)                               aCCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTGCGCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGACCCAGTTCact                                                                  >  1:438119/1‑98 (MQ=255)                               aCCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAAt                                                                  <  2:1572257/100‑1 (MQ=255)                                ccGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATa                                                                 >  2:1549748/1‑100 (MQ=255)                                ccGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATa                                                                 >  2:992018/1‑100 (MQ=255)                                ccGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATa                                                                 >  2:2356692/1‑100 (MQ=255)                                 cGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAg                                                                <  1:758973/100‑1 (MQ=255)                                 cGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAg                                                                <  2:1313217/100‑1 (MQ=255)                                  gggAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAgt                                                               >  1:2113299/1‑100 (MQ=255)                                   ggAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAgtg                                                              <  1:853119/100‑1 (MQ=255)                                    gAGCAAATGCAACGGCAGCACGGGATTTCGTCTCGCGGTCTGCCTCAATATGGGAATAGATTTAGCTATTTAAGATAGGGACGAACCAGTTCAATAgt                                                               <  1:543265/98‑1 (MQ=255)                                      gCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcg                                                           >  1:1826652/1‑100 (MQ=255)                                      gCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcg                                                           >  1:2869105/1‑100 (MQ=255)                                      gCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcg                                                           <  2:1162987/100‑1 (MQ=255)                                         atctcAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCa                                                        <  2:918141/96‑1 (MQ=255)                                          aTGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCaa                                                       >  1:1157994/1‑100 (MQ=255)                                          aTGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCaa                                                       <  2:204176/100‑1 (MQ=255)                                           tGCAACGGCAGCACGTGATTTCATCTCGCGCTCTTCCTCAATATGGTATTAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAAc                                                      <  2:2313623/100‑1 (MQ=255)                                           tGCAACGGCAGCACGTGATTTCATCTATCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCaa                                                       >  1:813562/1‑99 (MQ=255)                                             cAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCaagat                                                    >  1:918138/1‑97 (MQ=255)                                                  gcagcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTgg                                               <  1:2613194/100‑1 (MQ=255)                                                   cagcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGc                                              <  2:902598/100‑1 (MQ=255)                                                    agcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCt                                             <  1:2910873/100‑1 (MQ=255)                                                    agcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCt                                             <  2:2179065/100‑1 (MQ=255)                                                      caCGTGATTTCATCTATCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGa                                           <  2:2869124/100‑1 (MQ=255)                                                       aCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGAc                                          <  2:1036541/100‑1 (MQ=255)                                                         gTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACt                                         >  1:2437729/1‑99 (MQ=255)                                                          tGATTTCATCTCTCGCTCTTCCTCTATATGGTAGTAGATGTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGat                                          <  1:1177897/97‑2 (MQ=255)                                                          tGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTCGCTGACTTc                                       >  1:2485107/1‑100 (MQ=255)                                                             tttCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTcgcg                                    <  2:1573565/100‑1 (MQ=255)                                                              ttCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCgg                                   >  2:148701/1‑100 (MQ=255)                                                              ttCATCTCTCGCACCTCCTCAATATGGTAACAGATTTAGCTATTAAAGATAGGCACGAGCCAGTTCAATAATGCCGTCAACGGATTGGCTGACTTCGCgg                                   <  1:1802566/100‑1 (MQ=255)                                                               tCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGc                                  >  1:2751416/1‑100 (MQ=255)                                                               tCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGc                                  <  1:851256/100‑1 (MQ=255)                                                               tCATCTCTCGCTCTTCCTCAATACGGTTATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGc                                  <  1:2176529/100‑1 (MQ=255)                                                                cATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCt                                 >  2:2801883/1‑100 (MQ=255)                                                                  tctctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGt                               <  2:2384162/100‑1 (MQ=255)                                                                   ctctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTa                              >  2:73522/1‑100 (MQ=255)                                                                    tctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAg                             >  2:2144643/1‑100 (MQ=255)                                                                     ctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAg                             <  2:219536/99‑1 (MQ=255)                                                                     ctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGc                            >  2:2093846/1‑100 (MQ=255)                                                                       cGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAg                          >  2:2142082/1‑100 (MQ=255)                                                                       cGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAg                          >  2:744207/1‑100 (MQ=255)                                                                          tctTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTTTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAg                          <  1:2890976/97‑1 (MQ=255)                                                                           ctTCCTCAATATGGTAATAGTTTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCTTCAACGGATTGGCTGACTTCGCGGCTTTAGCAGTCAt                      >  1:2621127/1‑100 (MQ=255)                                                                           ctTCCTCAATATGGTAATAGATTTAGCTATTTAAGTTAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCa                       <  2:750165/99‑1 (MQ=255)                                                                             tCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCAttt                    >  2:1213574/1‑100 (MQ=255)                                                                              ccTCAATATGGTAATAGATGTGGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTc                   <  1:1329723/100‑1 (MQ=255)                                                                               cTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCTTTTCg                  >  2:1442288/1‑100 (MQ=255)                                                                               cTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCg                  >  1:1374748/1‑100 (MQ=255)                                                                               cTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCg                  <  2:1824727/100‑1 (MQ=255)                                                                                  aaTATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTc               <  1:2356677/100‑1 (MQ=255)                                                                                   atatGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCa              <  2:557400/100‑1 (MQ=255)                                                                                     atGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCaaa            >  2:207347/1‑100 (MQ=255)                                                                                      tGGTAATAGATTTAGCTATTTAATATAGGCATGAACCAGTTCAATAGTGGCGTGAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAAc           <  2:765462/100‑1 (MQ=255)                                                                                      tGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAAc           <  2:859148/100‑1 (MQ=255)                                                                                       ggTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACg          >  2:2585620/1‑100 (MQ=255)                                                                                        gTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGt         >  1:420113/1‑100 (MQ=255)                                                                                          aaTAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGttt       >  1:52111/1‑100 (MQ=255)                                                                                           aTAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTc      <  2:1033672/100‑1 (MQ=255)                                                                                           aTAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTc      <  1:1213569/100‑1 (MQ=255)                                                                                           aTAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGCAGCAGTCATTTCGGTCAAACGTTTc      <  1:555436/100‑1 (MQ=255)                                                                                            tAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTcc     >  2:1930484/1‑100 (MQ=255)                                                                                            tAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTcc     >  2:2374197/1‑100 (MQ=255)                                                                                               aTTTAGCTATTTAAGATATGCACGAACCAGTTCAATAGCGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCg   <  1:48222/99‑1 (MQ=255)                                                                                               aTTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCgg  >  2:1680441/1‑100 (MQ=255)                                                                                               aTTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCgg  <  2:2899447/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: ‑