Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	379,603	C→T	62.0%	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	62.0%	62.8 / 93.4	92	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 (20/15);  new base T (29/28);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

ACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGG  >  NC_000913/379507‑379699                                                                                                 |                                                                                                 aCGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGc                                                                                               <  2:1988424/100‑1 (MQ=255) aCGTCAATTTAAACGATTTACAGCGGTTTACCGGGAGCAAATGCAACGGCATAACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGc                                                                                               <  2:559352/100‑1 (MQ=255)   gTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTa                                                                                             >  1:1830298/1‑100 (MQ=255)      aaTTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTAttt                                                                                          <  1:1813367/100‑1 (MQ=255)       aTTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCACTATGGTAATAGATTTAGCTAttt                                                                                          >  1:771761/1‑99 (MQ=255)        tttCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTaa                                                                                        >  1:1729317/1‑100 (MQ=255)         ttCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAg                                                                                       >  1:507453/1‑100 (MQ=255)           cAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGAt                                                                                     >  1:1118218/1‑100 (MQ=255)            aaCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATa                                                                                    >  1:290584/1‑100 (MQ=255)             aCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAg                                                                                   >  2:1205726/1‑100 (MQ=255)             aCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAg                                                                                   >  1:173829/1‑100 (MQ=255)             aCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAg                                                                                   <  1:2134960/100‑1 (MQ=255)               gATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGc                                                                                 >  1:1135897/1‑100 (MQ=255)                aTTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCa                                                                                <  2:1980842/100‑1 (MQ=255)                aTTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCACCGCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCa                                                                                >  2:1377805/1‑100 (MQ=255)                 tttCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCAc                                                                               <  1:1046854/100‑1 (MQ=255)                 tttCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTATCGCTCTTCCGCAATATGGTAAGAGATTTAGCTATTTAAGATAgg                                                                                  >  2:1337268/1‑97 (MQ=255)                  ttCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACg                                                                              >  1:1737290/1‑100 (MQ=255)                     cAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAAc                                                                           >  1:1388264/1‑100 (MQ=255)                     cAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAACAGATTTAGCTATTTAAGATAGGCACGAAc                                                                           <  1:853460/100‑1 (MQ=255)                      aGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAAcc                                                                          >  1:139727/1‑100 (MQ=255)                      aGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAAcc                                                                          <  1:129207/100‑1 (MQ=255)                      aGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAAcc                                                                          <  1:1330305/100‑1 (MQ=255)                      aGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAAcc                                                                          >  1:372905/1‑100 (MQ=255)                       gCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCa                                                                         >  1:1944151/1‑100 (MQ=255)                        cGGTTTACCGGGAGCAAATGCAACGGCAGCAGGTTATTTCATGTCTCGCTCTTCCTCAAGATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAg                                                                        <  2:1545411/100‑1 (MQ=255)                         ggTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGt                                                                       >  2:890093/1‑100 (MQ=255)                          gTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGtt                                                                      >  2:1556315/1‑100 (MQ=255)                           tttACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTc                                                                     <  1:2012285/100‑1 (MQ=255)                           tttACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTc                                                                     <  1:2106062/100‑1 (MQ=255)                           tttACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTc                                                                     <  1:1615933/100‑1 (MQ=255)                             tACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCaa                                                                   <  2:861706/100‑1 (MQ=255)                               ccGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATa                                                                 >  2:1485128/1‑100 (MQ=255)                                cGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAg                                                                >  2:685077/1‑100 (MQ=255)                                cGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAg                                                                >  2:574617/1‑100 (MQ=255)                                 gggAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAgt                                                               >  2:1154617/1‑100 (MQ=255)                                   gAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAgtgt                                                             <  2:139727/100‑1 (MQ=255)                                   gAGCAAATGCAACGGCAGCACGTGATTTCATCGCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAgtg                                                              <  2:21309/99‑1 (MQ=255)                                     gCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcg                                                           <  1:999945/100‑1 (MQ=255)                                     gCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcg                                                           >  1:971293/1‑100 (MQ=255)                                      cAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcgt                                                          <  1:706350/100‑1 (MQ=255)                                      cAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcgt                                                          >  1:757680/1‑100 (MQ=255)                                       aaaTGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcgtc                                                         <  1:2007949/100‑1 (MQ=255)                                       aaaTGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcgtc                                                         <  2:1120811/100‑1 (MQ=255)                                       aaaTGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcgtc                                                         >  2:250459/1‑100 (MQ=255)                                          tGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAAc                                                      >  1:798680/1‑100 (MQ=255)                                           gCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACg                                                     >  1:594762/1‑100 (MQ=255)                                            cAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACgg                                                    <  2:1827569/100‑1 (MQ=255)                                             aaCGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACgg                                                    <  2:267708/99‑1 (MQ=255)                                               cGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGCCACGAACCAGTTCAATAGTGTCGTCAACGGAtt                                                 >  2:2149296/1‑100 (MQ=255)                                                ggCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGAtt                                                 <  1:755928/99‑1 (MQ=255)                                                ggCAGCACGTGATTTCATCTCCCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTg                                                <  1:485519/100‑1 (MQ=255)                                                  cagcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACTAACCAGTTCAATAGTGTCGTCAACGGATTGGc                                              <  2:300577/100‑1 (MQ=255)                                                    gcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGc                                              <  2:979268/98‑1 (MQ=255)                                                    gcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTg                                            <  1:643646/100‑1 (MQ=255)                                                    gcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGAtt                                                 >  1:262460/1‑95 (MQ=255)                                                       cGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATCGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGAc                                          >  1:1428286/1‑99 (MQ=255)                                                        gTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACtt                                        <  1:1488691/100‑1 (MQ=255)                                                         tGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTc                                       >  2:1539715/1‑100 (MQ=255)                                                          gATTTAATCTCTCGCTCTTCCTCAATATGGTAATAGAGTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTCAc                                          <  2:1712354/96‑1 (MQ=255)                                                           aTTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTcgc                                     >  2:1432695/1‑100 (MQ=255)                                                           aTTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTcgc                                     >  2:1750375/1‑100 (MQ=255)                                                           aTTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTc                                       <  2:1004810/98‑1 (MQ=255)                                                            tttCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTcgcg                                    <  2:466548/100‑1 (MQ=255)                                                                 tctctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGAGGCTGt                               >  1:812908/1‑100 (MQ=255)                                                                    ctcGCTCTTTCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGc                            <  1:1886728/100‑1 (MQ=255)                                                                    ctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGc                            <  2:751372/100‑1 (MQ=255)                                                                    ctcGCTCTTCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGc                            <  2:751238/100‑1 (MQ=255)                                                                      cGCTCTTCCTCAATATTGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAg                          <  1:634785/100‑1 (MQ=255)                                                                       gCTCTTCCTCAATATGGTAATAGATTTAGCTACTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAg                          >  2:1296087/1‑99 (MQ=255)                                                                       gCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGt                         >  2:177487/1‑100 (MQ=255)                                                                       gCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTAAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGt                         >  1:1608166/1‑100 (MQ=255)                                                                         tctTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCa                       >  1:1263666/1‑100 (MQ=255)                                                                            tCCTCAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCAttt                    >  1:151746/1‑100 (MQ=255)                                                                             ccTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTc                   <  1:205828/100‑1 (MQ=255)                                                                               tCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCgg                 <  1:2122626/100‑1 (MQ=255)                                                                               tCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCgg                 >  1:1044184/1‑100 (MQ=255)                                                                               tCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCgg                 <  2:2158356/100‑1 (MQ=255)                                                                                cAATATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGt                >  1:1011898/1‑100 (MQ=255)                                                                                cAATATGGTAATAGATTTAGCTATTTAAGATAGACACAAACAAGTTCAATAGTGTCGTCAACGGATTGGCTGACATCGCGGCTGTAGCAGTCATTTCGGt                >  2:529072/1‑100 (MQ=255)                                                                                cAATAGGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGATGTAGCAGTCATTTCgg                 <  1:2218982/99‑1 (MQ=255)                                                                                 aaTATGGTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGt                <  1:685077/99‑1 (MQ=255)                                                                                      ggTAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACg          <  2:335797/100‑1 (MQ=255)                                                                                        tAATAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGtt        >  2:509270/1‑100 (MQ=255)                                                                                        tAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGtt        >  1:1345638/1‑100 (MQ=255)                                                                                        tAATAGATTCAGCTATTTAAGATAGGCACGAACCAGATCAATAGTGTCGTCAACGGATTGGCTGACTTAGCGGCTGTAGCAGTCATGTCGGTCCAACGt         >  2:1671194/1‑99 (MQ=255)                                                                                           tAGATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTcc     >  1:2166659/1‑100 (MQ=255)                                                                                           tAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCATACGTTTc      <  1:1174414/99‑1 (MQ=255)                                                                                             gATTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCAGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCg   >  2:949611/1‑100 (MQ=255)                                                                                             gATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCg   >  1:669357/1‑100 (MQ=255)                                                                                              aTTTAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCgt  <  2:1669285/100‑2 (MQ=255)                                                                                              aTTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCgg  >  2:510055/1‑100 (MQ=255)                                                                                                 |                                                                                                 ACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGG  >  NC_000913/379507‑379699

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