Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	939,046	T→A	6.7%	T351T (ACT→ACA)	rarA →	recombination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	939,046	0	T	A	6.7%	227.7 / 5.5	74	T351T (ACT→ACA)	rarA	recombination factor
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base A (2/3);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

GACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTG  >  NC_000913/938897‑939194                                                                                                                                                      |                                                                                                                                                     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                                                                                                                                  <  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: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)                                                                                                                                                     cTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTAtgctg  >  2:1288460/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGATTGTTTACCTGGCCTGCGCGCCAAAAAGCAACGCTGTCTACACTGCGTTTAAAGCCGCGCTGGCCGATGCTCGCGAACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTG  >  NC_000913/938897‑939194

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