Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	262,059	A→G	5.1%	E186G (GAA→GGA)	proA →	gamma‑glutamylphosphate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	262,059	0	A	G	5.1%	250.2 / 3.1	79	E186G (GAA→GGA)	proA	gamma‑glutamylphosphate reductase
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base G (2/2);  total (42/37)
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.49e-01

GTAATGCGGTGATCCTGCGCGGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAG  >  NC_000913/261912‑262202                                                                                                                                                    |                                                                                                                                                gTAATGCGGTGATCCTGCGCGGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGa                                                                                                                                                 <  4:229742/148‑1 (MQ=255) gTAATGCGGTGATCCTGCGCGGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGa                                                                                                                                                 <  3:136390/148‑1 (MQ=255)      gCGGTGATCCTGCGCGGTGGCAAAGAAACGTGTCGGAGTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGGGGGTGCCGGGCGGGCGATTGATATTCCTGCCCGTGCGCTGGTCAGTGAAAtgc                                                                                                                                            >  2:269296/1‑148 (MQ=255)       cGGTGATCCTGCGCGGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGGGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAAtgct                                                                                                                                           <  7:296980/148‑1 (MQ=255)       cGGGGATCCTGCGCGGTGGCAAAGAAACGTGGCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAAtgct                                                                                                                                           <  3:215807/148‑1 (MQ=255)           gagccTGCGCGGTGGCACAGAAACGTGTCGCACTCACGCTGCAACGGGGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGt                                                                                                                                       <  8:383028/145‑1 (MQ=255)            aTCCTGCGCGGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGt                                                                                                                                       <  4:230119/147‑1 (MQ=255)            aTCCTGCGCGGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTa                                                                                                                                      <  2:149297/148‑1 (MQ=255)            aTCCTGCGCGGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTa                                                                                                                                      <  8:160744/148‑1 (MQ=255)                tGCGCGGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGt                                                                                                                                       >  4:32440/1‑143 (MQ=255)                tGCGCGGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGt                                                                                                                                       <  3:32440/143‑1 (MQ=255)                 gcgcGGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGAGTGATAATCCTGCCCGTGCGCTGGTCAGTGAAATGCTGCGTATGGa                                                                                                                                  >  5:190736/1‑147 (MQ=255)                  cgcgGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGGAATGCTGCGt                                                                                                                                       >  7:19794/1‑141 (MQ=255)                  cgcgGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGGAATGCTGCGt                                                                                                                                       <  8:19794/141‑1 (MQ=255)                  cgcgGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGGAATGCTGCGt                                                                                                                                       >  7:14279/1‑141 (MQ=255)                  cgcgGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGGAATGCTGCGt                                                                                                                                       <  8:14279/141‑1 (MQ=255)                      gTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATa                                                                                                                            >  8:184849/1‑148 (MQ=255)                          cAAAGAAACGTGTCGCACTAACGCTGCCACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATc                                                                                                                        >  6:308624/1‑148 (MQ=255)                            aaGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGa                                                                                                                      <  2:146090/148‑1 (MQ=255)                             aGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGa                                                                                                                      >  7:38411/1‑147 (MQ=255)                             aGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGAc                                                                                                                     <  4:173737/148‑1 (MQ=255)                                  cGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCt                                                                                                                >  7:344220/1‑148 (MQ=255)                                       cGCACTAACGCTGCAACGGTGGCGGTGATTCGGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATcc                                                                                                           >  8:198603/1‑148 (MQ=255)                                             aaCGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGaaa                                                                                                                                               <  6:106900/106‑1 (MQ=255)                                             aaCGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGaaa                                                                                                                                               >  5:106900/1‑106 (MQ=255)                                             aaCGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCgtg                                                                                                     >  7:232566/1‑148 (MQ=255)                                             aaCGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCgtg                                                                                                     >  7:361246/1‑148 (MQ=255)                                               cGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCgtggt                                                                                                   >  3:339977/1‑148 (MQ=255)                                               cGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCgtggt                                                                                                   >  1:319559/1‑148 (MQ=255)                                                 cTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGtggtgg                                                                                                 <  8:10767/148‑1 (MQ=255)                                                    cAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGgcgc                                                                                              >  1:141628/1‑148 (MQ=255)                                                      aCGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCt                                                                                             <  8:180329/147‑1 (MQ=255)                                                      aCGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCt                                                                                             >  7:180329/1‑147 (MQ=255)                                                           ggCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTg                                                                                       >  7:28423/1‑148 (MQ=255)                                                               gTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCAt                                                                                    >  1:28268/1‑147 (MQ=255)                                                                     cAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGCCCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACtgt                                                                             >  3:17628/1‑148 (MQ=255)                                                                      aGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACtgtg                                                                            <  4:339977/148‑1 (MQ=255)                                                                        gACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGcc                                                                          <  6:171619/148‑1 (MQ=255)                                                                          cGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGt                                                                        <  8:344220/148‑1 (MQ=255)                                                                           gCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGt                                                                        >  7:222274/1‑147 (MQ=255)                                                                            cccTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTg                                                                       >  5:196633/1‑147 (MQ=255)                                                                              cTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGGGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGaa                                                                     >  3:214798/1‑147 (MQ=255)                                                                                gAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAg                                                                  >  3:177230/1‑148 (MQ=255)                                                                                     ccTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGa                                                              >  8:324800/1‑147 (MQ=255)                                                                                     ccTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGa                                                              >  6:77230/1‑147 (MQ=255)                                                                                      cTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACa                                                            <  8:220064/148‑1 (MQ=255)                                                                                      cTGCGGCTTACCGGCGGGGGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGAc                                                             <  2:332180/147‑1 (MQ=255)                                                                                       tGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACaa                                                           >  2:57711/1‑148 (MQ=255)                                                                                        gCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGa                                                                      <  6:218581/136‑1 (MQ=255)                                                                                        gCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGa                                                                      >  5:218581/1‑136 (MQ=255)                                                                                         cGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACaa                                                           <  7:183966/146‑1 (MQ=255)                                                                                          ggCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATcc                                                        <  8:361246/148‑1 (MQ=255)                                                                                              tACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGt                                                    >  3:310750/1‑148 (MQ=255)                                                                                               aCCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGATCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGt                                                    >  2:131823/1‑147 (MQ=255)                                                                                                ccGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGa                                                  <  6:363941/148‑1 (MQ=255)                                                                                                       gTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAgg                                           <  7:198603/148‑1 (MQ=255)                                                                                                        tGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAggt                                          <  4:17628/148‑1 (MQ=255)                                                                                                           cGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAggtgg                                        >  3:58479/1‑147 (MQ=255)                                                                                                               cAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAg                                   <  5:308624/148‑1 (MQ=255)                                                                                                                  gCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCg                                >  1:313853/1‑148 (MQ=255)                                                                                                                   cGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGt                               <  8:232566/148‑1 (MQ=255)                                                                                                                     aTTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCg                                                                         <  7:57125/104‑1 (MQ=255)                                                                                                                     aTTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCg                                                                         >  8:57125/1‑104 (MQ=255)                                                                                                                         aTAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAgg                                  <  7:349066/139‑1 (MQ=255)                                                                                                                         aTAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAgg                                  <  7:235896/139‑1 (MQ=255)                                                                                                                         aTAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAgg                                  >  8:235896/1‑139 (MQ=255)                                                                                                                         aTAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAgg                                  >  8:349066/1‑139 (MQ=255)                                                                                                                           aaTCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCata                       >  8:199150/1‑148 (MQ=255)                                                                                                                            aTCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCatat                      >  3:151001/1‑148 (MQ=255)                                                                                                                             tCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATAtt                     >  1:360617/1‑148 (MQ=255)                                                                                                                              ccTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATAtt                     >  7:29543/1‑147 (MQ=255)                                                                                                                                  aCCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTAc                  <  4:177230/146‑1 (MQ=255)                                                                                                                                   ccGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGtt               >  3:288653/1‑148 (MQ=255)                                                                                                                                    cGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGtt               >  5:356076/1‑147 (MQ=255)                                                                                                                                         gcTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGaa         <  7:199150/148‑1 (MQ=255)                                                                                                                                             gTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAgtg     <  8:38411/148‑1 (MQ=255)                                                                                                                                             gTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAgt      <  7:324800/147‑1 (MQ=255)                                                                                                                                             gTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAgt      <  2:313853/147‑1 (MQ=255)                                                                                                                                              tCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGTTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAgtgt    <  5:77230/148‑1 (MQ=255)                                                                                                                                              tCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAgtgt    <  8:33672/148‑1 (MQ=255)                                                                                                                                              tCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAgtg     >  1:22236/1‑147 (MQ=255)                                                                                                                                                aGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTTTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTag  <  2:319559/148‑1 (MQ=255)                                                                                                                                                   gAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTa   >  5:188115/1‑144 (MQ=255)                                                                                                                                                    |                                                                                                                                                GTAATGCGGTGATCCTGCGCGGTGGCAAAGAAACGTGTCGCACTAACGCTGCAACGGTGGCGGTGATTCAGGACGCCCTGAAATCCTGCGGCTTACCGGCGGGTGCCGTGCAGGCGATTGATAATCCTGACCGTGCGCTGGTCAGTGAAATGCTGCGTATGGATAAATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAG  >  NC_000913/261912‑262202

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