Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	939,043	C→G	7.2%	Y350* (TAC→TAG)	rarA →	recombination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	939,043	0	C	G	7.2%	225.7 / 7.5	69	Y350* (TAC→TAG)	rarA	recombination factor
Reads supporting (aligned to +/- strand):  ref base C (35/29);  new base G (2/3);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GAAGACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGC  >  NC_000913/938894‑939192                                                                                                                                                      |                                                                                                                                                      gAAGACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTac                                                                                                                                                       >  1:423920/1‑150 (MQ=255)    gACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACt                                                                                                                                                    >  1:1540119/1‑150 (MQ=255)      cGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGTCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGc                                                                                                                                                  >  1:1629257/1‑150 (MQ=255)            tAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTaa                                                                                                                                            >  2:33735/1‑150 (MQ=255)             aaTGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCAGTGTAGACAGCGttgctt                                                                                                                                           <  1:1140241/150‑5 (MQ=255)               tGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGc                                                                                                                                         >  1:1178330/1‑150 (MQ=255)               tGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCAGTGTAGACAGCGttgctttt                                                                                                                                         <  1:201786/150‑7 (MQ=255)               tGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGACTGGGATTGCTTTACTCGCGTTGGCCCGGAGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGAGTTTAAAGc                                                                                                                                         >  2:578723/1‑150 (MQ=255)                 ccGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCcg                                                                                                                                       >  2:1280737/1‑150 (MQ=255)                  cGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCTCGTTGGCCCGGCGGAAGGTTAACGCGCCATTTCTCAGTCGATTGTTTACCTGGCCTGCGAGCAAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCcgc                                                                                                                                      <  1:578723/150‑1 (MQ=255)                   gATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCcgcg                                                                                                                                     >  2:102117/1‑150 (MQ=255)                     tCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCt                                                                                                                                   >  1:517165/1‑150 (MQ=255)                      ccACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGATGAACGCGCCAATGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTg                                                                                                                                  >  2:1344631/1‑150 (MQ=255)                       cACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCAGTGTAGACAGCGttgc                                                                                                                                             >  2:1207225/1‑136 (MQ=255)                       cACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCAGTGTAGACAGCGttgc                                                                                                                                             <  1:1207225/138‑3 (MQ=255)                         cGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGcc                                                                                                                               <  2:344095/150‑1 (MQ=255)                           ggCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGa                                                                                                                             >  1:1601699/1‑150 (MQ=255)                              gATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGc                                                                                                                          >  1:582622/1‑150 (MQ=255)                                  cAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGc                                                                                                                                >  2:319918/1‑140 (MQ=255)                                  cAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGc                                                                                                                                <  1:319918/140‑1 (MQ=255)                                  cAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTcgc                                                                                                                      >  1:2099/1‑150 (MQ=255)                                         cAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGcc                                                                                                               >  2:593306/1‑150 (MQ=255)                                             tGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGa                                                                                                           <  1:1272653/150‑1 (MQ=255)                                                ggCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTa                                                                                                        >  1:491113/1‑150 (MQ=255)                                                 gCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTAt                                                                                                       <  1:920002/150‑1 (MQ=255)                                                 gCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTAt                                                                                                       >  1:591497/1‑150 (MQ=255)                                                  ccTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATg                                                                                                      <  1:20678/150‑1 (MQ=255)                                                  ccTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATg                                                                                                      <  1:944608/150‑1 (MQ=255)                                                           gatttACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCgg                                                                                             >  1:64646/3‑150 (MQ=255)                                                           gCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTCTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCgg                                                                                             >  1:999953/1‑150 (MQ=255)                                                            cTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGt                                                                                            <  2:1394988/150‑1 (MQ=255)                                                            cTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGt                                                                                            <  2:466980/150‑1 (MQ=255)                                                             tttACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGtt                                                                                           <  2:831559/150‑1 (MQ=255)                                                             tttACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGtt                                                                                           <  2:878997/150‑1 (MQ=255)                                                               tACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCa                                                                                         <  1:237541/150‑1 (MQ=255)                                                                 cTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCAtt                                                                                       >  1:12661/1‑150 (MQ=255)                                                                     cgTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCg                                                                                   <  2:1259511/150‑1 (MQ=255)                                                                          gCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATg                                                                              <  2:1318492/150‑1 (MQ=255)                                                                          gCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTATACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATg                                                                              >  2:1325194/1‑150 (MQ=255)                                                                              ggcggAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGcc                                                                          <  2:729772/150‑1 (MQ=255)                                                                               gcggAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCcg                                                                         <  2:491113/150‑1 (MQ=255)                                                                                     ggTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGc                                                                                    >  1:1593107/1‑133 (MQ=255)                                                                                     ggTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGc                                                                                    <  2:1593107/133‑1 (MQ=255)                                                                                           cgcgCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATg                                                             <  2:1171684/150‑1 (MQ=255)                                                                                                    gCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATg                                                    >  1:901098/1‑150 (MQ=255)                                                                                                       cAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGc                                                 <  2:12661/150‑1 (MQ=255)                                                                                                         ggCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTa                                               >  1:1202728/1‑150 (MQ=255)                                                                                                          gCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGATGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTTCGTAATGCGCCGACGAAATTAATGAAGGAAATTGGCTAc                                              >  2:191526/1‑150 (MQ=255)                                                                                                           cGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACg                                             >  2:856282/1‑150 (MQ=255)                                                                                                             tttGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACggg                                           <  2:394714/149‑1 (MQ=255)                                                                                                              ttGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTCGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGc                                          <  1:1348140/150‑1 (MQ=255)                                                                                                                gTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAg                                        >  2:588593/1‑150 (MQ=255)                                                                                                                       tGGCCTGCGCGCCAAAAAGAAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATATGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATc                                 >  1:1064809/1‑150 (MQ=255)                                                                                                                           cTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTa                             <  1:593306/150‑1 (MQ=255)                                                                                                                            tGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTAc                            >  2:640385/1‑150 (MQ=255)                                                                                                                            tGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTAc                            >  2:423580/1‑150 (MQ=255)                                                                                                                            tGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTAc                            >  1:1045063/1‑150 (MQ=255)                                                                                                                              cgcgcCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGc                          <  2:2099/150‑1 (MQ=255)                                                                                                                               gcgcCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCt                         >  1:162339/1‑150 (MQ=255)                                                                                                                                 gcCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCa                       <  1:1468853/150‑1 (MQ=255)                                                                                                                                 gcCAAAAAGCAACGCGGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCa                       <  2:1467454/150‑1 (MQ=255)                                                                                                                                  ccAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACTCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCat                      <  2:423920/150‑1 (MQ=255)                                                                                                                                     gctttaaaCGCAGTGTAGACAGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCatgat                   >  2:1228679/7‑150 (MQ=255)                                                                                                                                          cAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGc              <  1:247464/150‑1 (MQ=255)                                                                                                                                              gCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAAc          >  2:12004/1‑150 (MQ=255)                                                                                                                                               cTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACg         >  1:177619/1‑150 (MQ=255)                                                                                                                                                 gTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCt       >  2:497038/1‑150 (MQ=255)                                                                                                                                                   cTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTa     >  2:842978/1‑150 (MQ=255)                                                                                                                                                    tACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTAt    <  1:1384885/150‑1 (MQ=255)                                                                                                                                                    tACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGAGTAATGCGCCGACGATATTAATGAAGGAAATGGGCTACGGTCAGGAAAATCGTTTCGGTCTTTATTAATCAAACGCTTAt    <  2:719025/150‑1 (MQ=255)                                                                                                                                                     acacTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGATATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTAtg   <  2:1402476/150‑1 (MQ=255)                                                                                                                                                     acacTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTAtg   >  2:1116503/1‑150 (MQ=255)                                                                                                                                                     acacTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTAtg   <  1:640385/150‑1 (MQ=255)                                                                                                                                                      cacTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCTAACGCCCTGATTATGACGTTCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGTAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTAtgc  <  1:935047/150‑1 (MQ=255)                                                                                                                                                      cacTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTAtgc  <  1:497038/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GAAGACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGC  >  NC_000913/938894‑939192

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