Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,834,765	C→A	6.1%	G389* (GGA→TGA)	yicJ ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,834,765	0	C	A	6.1%	262.8 / 4.4	75	G389* (GGA→TGA)	yicJ	putative transporter
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base A (3/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGTCGACGGTATCGGACATCAT  >  NC_000913/3834619‑3834906                                                                                                                                                   |                                                                                                                                              ggTTTTCAGATTGTGCGTCGTGAGAGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCg                                                                                                                                              <  4:86416/148‑1 (MQ=255)       cAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATa                                                                                                                                         <  2:382960/147‑1 (MQ=255)           ttGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGc                                                                                                                                      >  1:21957/1‑146 (MQ=255)                cgtcgtGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATc                                                                                                                               <  1:169133/148‑1 (MQ=255)                       agtgagtAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCa                                                                                                                                 >  4:159781/1‑139 (MQ=255)                       agtgagtAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCa                                                                                                                                 <  3:159781/139‑1 (MQ=255)                         tgagtAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCAGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATa                                                                                                                      <  4:255396/148‑1 (MQ=255)                          gagtAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGc                                                                                                                                      <  1:153452/131‑1 (MQ=255)                          gagtAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGc                                                                                                                                      >  2:153452/1‑131 (MQ=255)                              agtagCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAgc                                                                                                                 >  8:237028/1‑148 (MQ=255)                              agtagCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAgc                                                                                                                 >  6:47033/1‑148 (MQ=255)                               gtagCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAagag                                                                                                                  >  1:39814/1‑146 (MQ=255)                               gtagCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAagag                                                                                                                  <  2:39814/146‑1 (MQ=255)                                tagCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAgcgc                                                                                                               >  2:148186/1‑148 (MQ=255)                                 agCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGAt                                                                                                                       >  8:75516/1‑139 (MQ=255)                                 agCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGAt                                                                                                                       <  7:75516/139‑1 (MQ=255)                                 agCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATa                                                                                                                      >  1:395608/1‑140 (MQ=255)                                 agCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATa                                                                                                                      <  2:395608/140‑1 (MQ=255)                                 agCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGATcgcc                                                                                                              >  2:394625/1‑148 (MQ=255)                                   cGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGcgccgc                                                                                                            <  7:237028/148‑1 (MQ=255)                                       tAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGccgccg                                                                                                          >  3:10813/1‑146 (MQ=255)                                            aTAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCa                                                                                                   <  4:10813/148‑1 (MQ=255)                                             tAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATATCGCAATAATGATGCTTATCGTGGCGCTGTGCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGccga                                                                                                  >  6:211207/1‑146 (MQ=255)                                             tAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCaa                                                                                                  >  1:351553/1‑148 (MQ=255)                                                 cgcgCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGccgccg                                                                                                          >  2:355147/1‑136 (MQ=255)                                                 cgcgCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGccgccg                                                                                                          <  1:355147/136‑1 (MQ=255)                                                  gcgcTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAAcccc                                                                                              <  1:148186/147‑1 (MQ=255)                                                       cAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGttt                                                                                         >  7:262696/1‑147 (MQ=255)                                                         gCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAg                                                                                      >  3:236309/1‑148 (MQ=255)                                                            aaTAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGca                                                                                    >  8:125146/1‑147 (MQ=255)                                                              tAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACaa                                                                                 <  2:351553/148‑1 (MQ=255)                                                              tAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACaa                                                                                 <  5:47033/148‑1 (MQ=255)                                                               aaCAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGAt                                                                                                                       <  6:10662/109‑1 (MQ=255)                                                               aaCAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGAt                                                                                                                       >  5:10662/1‑109 (MQ=255)                                                               aaCAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCa                                                                                                   >  7:103882/1‑129 (MQ=255)                                                               aaCAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCa                                                                                                   <  8:103882/129‑1 (MQ=255)                                                                     tCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCAGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGt                                                                          <  1:394625/148‑1 (MQ=255)                                                                        ccGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGcc                                                                       <  4:40019/148‑1 (MQ=255)                                                                             aCAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCaa                                                                  <  7:376310/148‑1 (MQ=255)                                                                              cAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCaaa                                                                 <  2:21957/148‑1 (MQ=255)                                                                                aTCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAAc                                                               <  6:197641/148‑1 (MQ=255)                                                                                    tGAATAGCGCAATAATGATGCTAATAGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCAGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGAt                                                           >  7:62221/1‑148 (MQ=255)                                                                                    tGAATAGCGCAATAATGATGCTAATAGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCAGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGAt                                                           >  7:55907/1‑148 (MQ=255)                                                                                         aGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTg                                                                                        >  4:212173/1‑114 (MQ=255)                                                                                         aGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTg                                                                                        <  3:212173/114‑1 (MQ=255)                                                                                           cgcAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGcc                                                                       >  8:369744/1‑129 (MQ=255)                                                                                           cgcAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGcc                                                                       <  7:369744/129‑1 (MQ=255)                                                                                            gcAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCAGCCATAAGCCAGCATCCAGCAGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGAACAAACAGCGTGCCAGCAAAAATGATCCCGGt                                                     >  4:346026/1‑146 (MQ=255)                                                                                                taatGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTc                                                    >  8:319315/1‑143 (MQ=255)                                                                                                taatGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTc                                                    <  7:319315/143‑1 (MQ=255)                                                                                                taatGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCg                                               >  5:316149/1‑148 (MQ=255)                                                                                                 aatGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTc                                                    >  7:164881/1‑142 (MQ=255)                                                                                                 aatGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTc                                                    <  8:164881/142‑1 (MQ=255)                                                                                                 aatGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGt                                              >  3:233050/1‑148 (MQ=255)                                                                                                      tgCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTAcc                                         <  4:236309/148‑1 (MQ=255)                                                                                                           aTCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGc                                    >  8:205631/1‑148 (MQ=255)                                                                                                           aTCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGAACCCGGTCAGGCGTTTACCATTGc                                    >  2:386816/1‑148 (MQ=255)                                                                                                             cGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCaaa                                                                 <  1:27455/117‑1 (MQ=255)                                                                                                             cGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCaaa                                                                 >  2:27455/1‑117 (MQ=255)                                                                                                              gTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGcac                                  >  6:356595/1‑147 (MQ=255)                                                                                                                           gcgcTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGCGCCAGCAAAACTGATCCCgg                                                      >  7:389857/1‑114 (MQ=255)                                                                                                                           gcgcTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGCGCCAGCAAAACTGATCCCgg                                                      <  8:389857/114‑1 (MQ=255)                                                                                                                            cgcTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGTCg                   <  7:125146/148‑1 (MQ=255)                                                                                                                            cgcTTTTTCTGCCGCATCATATCAGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGt                     <  8:55907/146‑1 (MQ=255)                                                                                                                            cgcTTTTTCTGCCGCATCATATCAGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGt                     <  8:62221/146‑1 (MQ=255)                                                                                                                                 tttCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGTCGACg                <  6:316149/146‑1 (MQ=255)                                                                                                                                  ttCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGTCGACGGTa             >  4:80629/1‑148 (MQ=255)                                                                                                                                    cTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGTCGACGGTAt            <  1:386816/147‑1 (MQ=255)                                                                                                                                       ccGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGTCGACGGTATCgg         <  8:262696/147‑1 (MQ=255)                                                                                                                                           atcatATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGTCGACGGTATCGGAcatc    <  3:80629/148‑1 (MQ=255)                                                                                                                                            tcatATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGTCGACGGTATCGGAcatca   <  7:205631/148‑1 (MQ=255)                                                                                                                                             catATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGTCGACGGTATCGGAcatcat  >  5:356281/1‑148 (MQ=255)                                                                                                                                             catATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGTCGACGGTATCGGAcatcat  >  2:259753/1‑148 (MQ=255)                                                                                                                                                atCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGa                                                                                       >  7:245655/1‑60 (MQ=255)                                                                                                                                                atCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGa                                                                                       <  8:245655/60‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              GGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCCAACCCCAGTTTGAGCACAAACAGCGTGCCAGCAAAACTGATCCCGGTCAGGCGTTTACCATTGCACCACTCGCCGTAGTCGACGGTATCGGACATCAT  >  NC_000913/3834619‑3834906

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