Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	939,037	T→A	7.1%	A348A (GCT→GCA)	rarA →	recombination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	939,037	0	T	A	7.1%	214.6 / 5.1	70	A348A (GCT→GCA)	rarA	recombination factor
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base A (2/3);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CGTCTGAAGACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACG  >  NC_000913/938889‑939185                                                                                                                                                     |                                                                                                                                                     cGTCTGAAGACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTg                                                                                                                                                     <  1:1049776/150‑1 (MQ=255) cGTCTGAAGACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTg                                                                                                                                                     <  1:534922/150‑1 (MQ=255)    cTGAAGACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCt                                                                                                                                                  <  1:239510/150‑1 (MQ=255)     tGAAGACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTa                                                                                                                                                 <  1:6302/150‑1 (MQ=255)      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                                                                                                                                      <  1:1207225/138‑3 (MQ=255)                            cACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCAGTGTAGACAGCGttgc                                                                                                                                      >  2:1207225/1‑136 (MQ=255)                              cGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGcc                                                                                                                        <  2:344095/150‑1 (MQ=255)                                ggCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGa                                                                                                                      >  1:1601699/1‑150 (MQ=255)                                   gATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGc                                                                                                                   >  1:582622/1‑150 (MQ=255)                                       cAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGc                                                                                                                         <  1:319918/140‑1 (MQ=255)                                       cAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGc                                                                                                                         >  2:319918/1‑140 (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:591497/1‑150 (MQ=255)                                                      gCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTAt                                                                                                <  1:920002/150‑1 (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                                                                             <  2:1593107/133‑1 (MQ=255)                                                                                          ggTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGc                                                                             >  1:1593107/1‑133 (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:423580/1‑150 (MQ=255)                                                                                                                                 tGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTAc                     >  2:640385/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)                                                                                                                                                     |                                                                                                                                                     CGTCTGAAGACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACG  >  NC_000913/938889‑939185

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