Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	832,728	C→T	5.3%	A75T (GCA→ACA)	ybiA ←	DUF1768 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	832,728	0	C	T	5.3%	242.3 / 2.3	75	A75T (GCA→ACA)	ybiA	DUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base C (34/37);  new base T (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

GTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCA  >  NC_000913/832584‑832870                                                                                                                                                 |                                                                                                                                               gTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTgcggc                                                                                                                                             <  8:94423/148‑1 (MQ=255)    ttGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGa                                                                                                                                           <  8:432717/147‑1 (MQ=255)    ttGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGa                                                                                                                                           <  7:123001/147‑1 (MQ=255)         ggCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAg                                                                                                                                     >  5:49449/1‑148 (MQ=255)         ggCGCGGTTGCCAGCAAGAGCGAACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAg                                                                                                                                     >  6:144588/1‑148 (MQ=255)             cgGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGga                                                                                                                                 <  1:260402/148‑1 (MQ=255)                  gCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGgaagaa                                                                                                                             <  3:180838/147‑1 (MQ=255)                      gCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAAccc                                                                                                                         >  3:386856/1‑147 (MQ=255)                       cAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGa                                                                                                                       >  4:351217/1‑148 (MQ=255)                                    aGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTcgc                                                                                                          <  4:121884/148‑1 (MQ=255)                                            aTGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGTGGCGACCATAGGGGAAGAAACCCGACGGa                                                                                                                   >  1:115569/1‑131 (MQ=255)                                            aTGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGTGGCGACCATAGGGGAAGAAACCCGACGGa                                                                                                                   <  2:115569/131‑1 (MQ=255)                                              gCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGAtctc                                                                                                               <  2:388002/133‑1 (MQ=255)                                              gCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGAtctc                                                                                                               >  1:388002/1‑133 (MQ=255)                                              gCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTAttttt                                                                                                 >  8:304133/1‑147 (MQ=255)                                               cTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTtca                                                                                               >  6:256994/1‑148 (MQ=255)                                                 gTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTtcat                                                                                              >  4:42187/1‑147 (MQ=255)                                                   tCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTtcatc                                                                                             >  8:387301/1‑146 (MQ=255)                                                      aaTTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCTTAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGa                                                                                        <  8:405645/148‑1 (MQ=255)                                                        tttAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCaa                                                                                           <  2:316009/143‑1 (MQ=255)                                                        tttAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCaa                                                                                           >  1:316009/1‑143 (MQ=255)                                                            gcgcGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAAttttt                                                                                  <  7:304133/148‑1 (MQ=255)                                                            gcgcGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAAtttt                                                                                   <  3:351217/147‑1 (MQ=255)                                                             cgcgAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTtg                                                                                 <  5:43253/148‑1 (MQ=255)                                                                  aGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCt                                                                            >  8:406735/1‑148 (MQ=255)                                                                    tGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTg                                                                           >  8:385044/1‑147 (MQ=255)                                                                     gCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGgaagaa                                                                                                                             <  1:10218/96‑1 (MQ=255)                                                                     gCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGgaagaa                                                                                                                             >  2:10218/1‑96 (MQ=255)                                                                       ttttCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTAtcttt                                                                                                 >  6:117933/1‑122 (MQ=255)                                                                       ttttCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTAtcttt                                                                                                 <  5:117933/122‑1 (MQ=255)                                                                       ttttCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTtc                                                                                                >  6:205520/1‑123 (MQ=255)                                                                       ttttCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTtc                                                                                                <  5:205520/123‑1 (MQ=255)                                                                       ttttCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAAttttt                                                                                  >  3:1672/1‑137 (MQ=255)                                                                       ttttCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAAttttt                                                                                  <  4:1672/137‑1 (MQ=255)                                                                         ttCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGaaaa                                                                       >  3:313977/1‑146 (MQ=255)                                                                         ttCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAAt                                                                      <  3:42187/147‑1 (MQ=255)                                                                            gCATCAACTTGTCCTTGACCGACTTCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAgt                                                                   <  3:312441/147‑1 (MQ=255)                                                                                cACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCt                                                              >  1:337775/1‑148 (MQ=255)                                                                                  ccTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGa                                                            <  7:385044/148‑1 (MQ=255)                                                                                   cTGTTCTTTGACCGTCTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCTCACCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACTGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAg                                                           <  3:225051/148‑1 (MQ=255)                                                                                   cTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGGGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAg                                                           <  7:406735/148‑1 (MQ=255)                                                                                   cTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAg                                                           >  2:183307/1‑148 (MQ=255)                                                                                   cTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAg                                                           >  8:412049/1‑148 (MQ=255)                                                                                    tGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAAtttt                                                                                   >  4:387301/1‑123 (MQ=255)                                                                                    tGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAAtttt                                                                                   <  3:387301/123‑1 (MQ=255)                                                                                      ttCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGt                                                         <  4:313977/147‑1 (MQ=255)                                                                                      ttCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGAGCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTa                                                        <  6:23185/148‑1 (MQ=255)                                                                                         tttGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGc                                                     >  5:105068/1‑148 (MQ=255)                                                                                           tGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAgg                                                      <  7:387301/145‑1 (MQ=255)                                                                                            gACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCa                                                   <  5:144588/147‑1 (MQ=255)                                                                                               cGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGTGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAg                                                  <  8:293310/145‑1 (MQ=255)                                                                                               cGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGTGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAg                                                  >  7:293310/1‑145 (MQ=255)                                                                                                     ccAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGTAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCg                                         <  7:277643/148‑1 (MQ=255)                                                                                                     ccAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCg                                         >  8:446933/1‑148 (MQ=255)                                                                                                      cAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAAttt                                                                                    <  3:342936/104‑1 (MQ=255)                                                                                                      cAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAAttt                                                                                    >  4:342936/1‑104 (MQ=255)                                                                                                        gTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCa                                      <  4:386856/148‑1 (MQ=255)                                                                                                         tttttACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCaa                                     >  5:24752/1‑148 (MQ=255)                                                                                                          ttttACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAg                                                           >  6:49170/1‑125 (MQ=255)                                                                                                          ttttACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAg                                                           <  5:49170/125‑1 (MQ=255)                                                                                                               cGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTg                                <  7:446933/147‑1 (MQ=255)                                                                                                               cGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTg                                <  2:337775/147‑1 (MQ=255)                                                                                                               cGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTg                                <  6:49449/147‑1 (MQ=255)                                                                                                                 cAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAAt                                                                      <  1:391114/107‑1 (MQ=255)                                                                                                                 cAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAAt                                                                      >  2:391114/1‑107 (MQ=255)                                                                                                                    agGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGc                          >  6:111432/1‑148 (MQ=255)                                                                                                                         tAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAggc                     >  3:286575/1‑148 (MQ=255)                                                                                                                          aGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGGAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAgg                      <  1:183307/146‑1 (MQ=255)                                                                                                                           gAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAggcgg                   >  3:19141/1‑148 (MQ=255)                                                                                                                             aCGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAggcggc                  >  7:422452/1‑147 (MQ=255)                                                                                                                               gATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCa                 >  5:233740/1‑146 (MQ=255)                                                                                                                                aTCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAg              <  4:19141/148‑1 (MQ=255)                                                                                                                                aTCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTAAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAg              >  1:291866/1‑148 (MQ=255)                                                                                                                                   gcgcCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGcc                         >  6:313932/1‑134 (MQ=255)                                                                                                                                   gcgcCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGcc                         <  5:313932/134‑1 (MQ=255)                                                                                                                                     gcCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGa         <  5:119262/148‑1 (MQ=255)                                                                                                                                            gcgtgcgGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCa  >  3:220319/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                               GTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCA  >  NC_000913/832584‑832870

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