Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	367,830	G→A	18.9%	Y197Y (TAC→TAT)	mhpR ←	mhp operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	367,830	0	G	A	18.9%	177.0 / 32.0	84	Y197Y (TAC→TAT)	mhpR	mhp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (37/30);  new base A (10/6);  total (47/37)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

CACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATT  >  NC_000913/367731‑367929                                                                                                    |                                                                                                    cacCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCg                                                                                                     >  2:1366035/1‑100 (MQ=255) cacCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCg                                                                                                     <  1:1723376/100‑1 (MQ=255) cacCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGACGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCg                                                                                                     >  1:1594949/1‑100 (MQ=255)   ccAGATTCAGACAGCCAATTACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTa                                                                                                   <  1:1913583/100‑1 (MQ=255)   ccAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATa                                                                                                   >  1:1294187/1‑100 (MQ=255)   ccAGATTCAGACAGCCAATCAAACGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGACGCGGTAGTTCTGTCCATa                                                                                                   >  2:210251/1‑100 (MQ=255)    cAGATTCAGACAGCCAATCACCCTTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTaa                                                                                                  <  1:1144651/100‑1 (MQ=255)    cAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATTTTCTCCTCCTGATCCCAGCCGCGGTAGTTCGTTCCGTa                                                                                                   >  1:921607/1‑99 (MQ=255)     aGATTCAGACAGCCAATCACCCGTTGCTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAAc                                                                                                 <  1:839228/100‑1 (MQ=255)      gATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAAcc                                                                                                >  2:115235/1‑100 (MQ=255)       aTTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAAccc                                                                                               <  1:1365386/100‑1 (MQ=255)            gACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCtt                                                                                           >  2:2109952/1‑99 (MQ=255)             aCAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTg                                                                                         <  2:582270/100‑1 (MQ=255)                gCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTCGTTCTGTCCGTAACCCTCTTTgcg                                                                                       >  2:1342245/1‑99 (MQ=255)                  cAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcg                                                                                     >  2:303258/1‑99 (MQ=255)                   aaTCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcg                                                                                   >  1:1294063/1‑100 (MQ=255)                      cACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgcgc                                                                                >  1:213977/1‑100 (MQ=255)                       aCCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTgcgcgcgcgcg                                                                               >  2:1223456/1‑100 (MQ=255)                         ccGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGcc                                                                             >  2:677209/1‑100 (MQ=255)                          cGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCa                                                                            >  2:1332892/1‑100 (MQ=255)                             tgttCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGaa                                                                         >  2:1850088/1‑100 (MQ=255)                               ttCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATg                                                                       <  2:1122127/100‑1 (MQ=255)                               ttCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATg                                                                       >  2:1370854/1‑100 (MQ=255)                                 cACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTTCGCGCGCGCGCCAGAATGGc                                                                     >  1:2354104/1‑100 (MQ=255)                                 cACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGc                                                                     >  1:146692/1‑100 (MQ=255)                                 cACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGc                                                                     >  2:1564483/1‑100 (MQ=255)                                  aCTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCt                                                                    >  1:847389/1‑100 (MQ=255)                                     gcgcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTcc                                                                 <  2:342355/100‑1 (MQ=255)                                     gcgcAGCGGTACGGCGATAAAGGCGATCTTCTCCTCCTGATCCCAGCCGCGCTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTcc                                                                 >  2:1046431/1‑100 (MQ=255)                                       gcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAGCCCTCTTTGCGCGCGCGCGCCAGAATGGCTTcc                                                                 <  2:1585412/98‑1 (MQ=255)                                       gcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAg                                                               <  1:1630891/100‑1 (MQ=255)                                         aGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCTCGCGCCAGAATGGCTTCCAg                                                               <  2:1457594/98‑1 (MQ=255)                                         aGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCt                                                             <  1:1959382/100‑1 (MQ=255)                                           cGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCttt                                                           <  2:2338833/100‑1 (MQ=255)                                            ggTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTa                                                          <  2:370573/100‑1 (MQ=255)                                             gTATGGCGATAGAGGCGAACTTCTCCTCATGATCACCGCCGCGGTAGTTCTGTCCGTAAGCCTCTTTGTGCGCGCACGCCAGAATGTCTTCCAGCTTTa                                                          >  2:293929/1‑99 (MQ=255)                                              tACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAAc                                                        >  1:100259/1‑100 (MQ=255)                                               aCGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACg                                                       <  2:887534/100‑1 (MQ=255)                                               aCGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACg                                                       >  1:2356538/1‑100 (MQ=255)                                                cGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACgg                                                      <  1:346353/100‑1 (MQ=255)                                                 ggCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGt                                                     >  1:223666/1‑100 (MQ=255)                                                   cGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGTGCGCCAGAATGGCTTCCAGCTTTAACGGTTc                                                   <  2:471227/100‑1 (MQ=255)                                                    gATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTcc                                                  <  1:872390/100‑1 (MQ=255)                                                      tAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCg                                                <  2:431452/100‑1 (MQ=255)                                                        gagGCGATCTTCGCCACCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTg                                              <  1:690797/100‑1 (MQ=255)                                                         agGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGc                                             >  2:2324739/1‑100 (MQ=255)                                                         agGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGc                                             >  2:430076/1‑100 (MQ=255)                                                          ggCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGcc                                            <  1:779384/100‑1 (MQ=255)                                                          ggCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGcc                                            >  1:1412095/1‑100 (MQ=255)                                                          ggCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGcc                                            >  2:2021175/1‑100 (MQ=255)                                                          ggCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGc                                             >  1:2125253/1‑99 (MQ=255)                                                            cGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAg                                          >  1:2108211/1‑100 (MQ=255)                                                             gATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGt                                         >  2:607858/1‑100 (MQ=255)                                                               tcttgtCCTCCTGCTCCCAGCCGCGGTAGTTCTGTCCGTAACCCCCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGtt                                        <  1:1342259/99‑1 (MQ=255)                                                               tcttctCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTg                                       <  2:1978171/100‑1 (MQ=255)                                                                cttctCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGa                                      <  1:1900114/100‑1 (MQ=255)                                                                  tctcCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCTCGCGCGCTCCAGAATGGCTTCCAGCTTTAACGCTTCCCGTGCCAGtt                                        >  2:837285/1‑96 (MQ=255)                                                                  tctcCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATa                                    >  1:696418/1‑100 (MQ=255)                                                                    tcctccTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAg                                   <  2:1144487/99‑1 (MQ=255)                                                                     cctcctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGtc                                 >  2:1862304/1‑100 (MQ=255)                                                                      ctcctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTCACGGTTCCCGTGCCAGTTGATAGtca                                >  2:432508/1‑100 (MQ=255)                                                                        cctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCCCTTTGCGCCCCCGCGCAAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGacatc                              >  2:1638549/1‑100 (MQ=255)                                                                         ctGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGtcatca                             >  2:666026/1‑100 (MQ=255)                                                                          tGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCAc                            <  2:1234854/100‑1 (MQ=255)                                                                           gATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGtcatca                             >  2:349136/1‑98 (MQ=255)                                                                              cccAGCCGCGGTAGTTCTGTCCGTATCCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCggg                        <  2:1764670/100‑1 (MQ=255)                                                                              cccAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTTCCGTGCCAGTTGATAGTCATCACCggg                        <  2:1074528/100‑1 (MQ=255)                                                                              cccAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCggg                        >  1:1324129/1‑100 (MQ=255)                                                                              cccAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCggg                        >  2:339390/1‑100 (MQ=255)                                                                               ccAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCAACggg                        >  1:1641003/1‑99 (MQ=255)                                                                               ccAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCgggg                       >  2:1992766/1‑100 (MQ=255)                                                                                   ccGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCggg                   <  2:146689/100‑1 (MQ=255)                                                                                    cgcgGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCggg                   <  2:899677/99‑1 (MQ=255)                                                                                     gcgGTAGTTCTGTCCTTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAg                 <  1:2364230/100‑1 (MQ=255)                                                                                     gcgGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAg                 <  1:178050/100‑1 (MQ=255)                                                                                     gcgGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAg                 <  1:1563964/100‑1 (MQ=255)                                                                                      cgGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAgg                >  2:756782/1‑100 (MQ=255)                                                                                       ggTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGc               <  1:1568044/100‑1 (MQ=255)                                                                                           gTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAAc           <  1:756791/100‑1 (MQ=255)                                                                                           gTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAAc           >  1:52990/1‑100 (MQ=255)                                                                                            ttCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACa          <  1:2308836/100‑1 (MQ=255)                                                                                              cTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACAtt        >  2:1666295/1‑100 (MQ=255)                                                                                                gTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTc      <  1:671351/100‑1 (MQ=255)                                                                                                 tCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCg     >  2:1933953/1‑100 (MQ=255)                                                                                                 tCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCg     >  1:1836696/1‑100 (MQ=255)                                                                                                  ccGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGa    <  2:872495/100‑1 (MQ=255)                                                                                                   cGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGAt   >  1:1460302/1‑100 (MQ=255)                                                                                                    gTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGAtt  >  2:1185431/1‑100 (MQ=255)                                                                                                    |                                                                                                    CACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATT  >  NC_000913/367731‑367929

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: ‑