Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,561,293	G→C	6.2%	R54G (CGT→GGT)	glpG ←	rhomboid protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,561,293	0	G	C	6.2%	217.7 / 4.0	65	R54G (CGT→GGT)	glpG	rhomboid protease
Reads supporting (aligned to +/- strand):  ref base G (26/35);  new base C (2/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GGCGATCATCATCACCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATCAACAAACGCCTGCGCCACGCGGGGGTTAGCAAAAGAGGTA  >  NC_000913/3561144‑3561441                                                                                                                                                      |                                                                                                                                                     ggCGATCATCATCACCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAAcg                                                                                                                                                      >  1:1013715/1‑150 (MQ=255)   cGATCATCATCACCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAAcgcg                                                                                                                                                    >  1:649508/1‑150 (MQ=255)   cGATCATCATCACCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAAcgcg                                                                                                                                                    >  1:639086/1‑150 (MQ=255)   cGATCATCATCACCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAAcgcg                                                                                                                                                    >  1:1249790/1‑150 (MQ=255)            tcaCCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCg                                                                                                                                           <  2:669489/150‑1 (MQ=255)             caCCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGc                                                                                                                                          <  2:729829/150‑1 (MQ=255)             caCCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGc                                                                                                                                          >  2:524575/1‑150 (MQ=255)             caCCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGc                                                                                                                                          >  2:845993/1‑150 (MQ=255)             caCCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGc                                                                                                                                          <  2:96225/150‑1 (MQ=255)             caCCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGc                                                                                                                                          >  1:148674/1‑150 (MQ=255)              aCCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGcc                                                                                                                                         <  1:1245256/150‑1 (MQ=255)              aCCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGcc                                                                                                                                         <  1:643986/150‑1 (MQ=255)              aCCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGcc                                                                                                                                         <  2:1273278/150‑1 (MQ=255)               cccAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGccc                                                                                                                                        >  1:1116603/1‑150 (MQ=255)                 cAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGt                                                                                                                                      >  2:559312/1‑150 (MQ=255)                    gTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAcg                                                                                                                                   <  2:663092/150‑1 (MQ=255)                    gTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAcg                                                                                                                                   >  1:79811/1‑150 (MQ=255)                     ttACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAcgc                                                                                                                                  <  2:593397/150‑1 (MQ=255)                             ccTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGcc                                                                                                                          <  1:1245010/150‑1 (MQ=255)                             ccTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGcc                                                                                                                          <  2:288065/150‑1 (MQ=255)                                gcgcGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGAGCTGGCGCGTTTTCTCGAAAACCCGGCAGATCCGcg                                                                                                                                         >  2:1223671/1‑131 (MQ=37)                                gcgcGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGAGCTGGCGCGTTTTCTCGAAAACCCGGCAGATCCGcg                                                                                                                                         <  1:1223671/132‑2 (MQ=37)                                 cgcgTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTggg                                                                                                                      <  1:524575/150‑1 (MQ=255)                                 cgcgTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTggg                                                                                                                      >  2:957063/1‑150 (MQ=255)                                       cACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCAGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCAt                                                                                                                >  2:72267/1‑150 (MQ=255)                                          gCAAGGCGGCAAAGAAAGGATAACGTCGATAATGCAGGCCTCTTCAGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGTATCTGCCGGGTTTTCGAGAAAACGCGCGAGCTCCGCCCGTACGCTCTCGGCCTTGGACTCATCCg                                                                                                             <  2:480047/150‑1 (MQ=255)                                          gCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCg                                                                                                             <  1:640675/150‑1 (MQ=255)                                                   cAAAGAAAGGATAAAGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGa                                                                                                    >  2:1025673/1‑150 (MQ=255)                                                    aaagaaagGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGAc                                                                                                   <  1:1005760/150‑1 (MQ=255)                                                      agaaagGATAACGGCGATAATGCAGGTCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACAt                                                                                                 >  2:976305/1‑150 (MQ=255)                                                      agaaagGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACAt                                                                                                 >  1:10894/1‑150 (MQ=255)                                                       gaaagGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATc                                                                                                >  1:450366/1‑150 (MQ=255)                                                         aagGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGc                                                                                              <  1:833435/150‑1 (MQ=255)                                                             aTAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTg                                                                                          <  2:1249790/150‑1 (MQ=255)                                                             aTAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTg                                                                                          <  2:1233695/150‑1 (MQ=255)                                                               aaCGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGt                                                                                        >  2:997252/1‑150 (MQ=255)                                                                 cGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTa                                                                                      >  1:584084/1‑150 (MQ=255)                                                                        caggCAGGCCACTGCCGGTATGGCCTGCCTGACAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATgttgtt                                                                               >  2:939398/4‑150 (MQ=255)                                                                         aTGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGttgttg                                                                              <  2:10894/150‑1 (MQ=255)                                                                         aTGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGttgttg                                                                              <  1:58890/150‑1 (MQ=255)                                                                         aTGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGttgttg                                                                              <  1:248482/150‑1 (MQ=255)                                                                          tGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGa                                                                             >  1:803720/1‑150 (MQ=255)                                                                                 cACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGa                                                                      <  2:317040/150‑1 (MQ=255)                                                                                    tGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCTTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGGTGTTGAATCGTGAGGa                                                                   <  2:478090/150‑1 (MQ=255)                                                                                        ggTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAac                                                               >  1:261477/1‑150 (MQ=255)                                                                                          tATGGCCTGCCTGCCAGCTCGCCGCCATATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAacag                                                             <  2:639086/150‑2 (MQ=255)                                                                                           aTGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACAcc                                                            >  1:358825/1‑150 (MQ=255)                                                                                           aTGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACAcc                                                            >  1:353578/1‑150 (MQ=255)                                                                                            tGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACAccc                                                           >  1:788983/1‑150 (MQ=255)                                                                                             ggCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCt                                                          >  1:766684/1‑150 (MQ=255)                                                                                                   cctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCACACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCg                                                    <  1:324015/150‑1 (MQ=255)                                                                                                      gccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCa                                                 <  2:353578/150‑1 (MQ=255)                                                                                                      gccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCa                                                 <  2:358825/150‑1 (MQ=255)                                                                                                      gccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCa                                                 <  2:1013715/150‑1 (MQ=255)                                                                                                         aGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGt                                              >  2:1138226/1‑150 (MQ=255)                                                                                                          gCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTa                                             >  2:1357415/1‑150 (MQ=255)                                                                                                           cTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTaa                                            >  2:1405641/1‑150 (MQ=255)                                                                                                           cTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTaa                                            <  2:1116603/150‑1 (MQ=255)                                                                                                             cgccgcCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATc                                          <  2:79811/150‑1 (MQ=255)                                                                                                              gccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATca                                         <  1:957063/150‑1 (MQ=255)                                                                                                               ccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACCCCCTGCGTCGCCATGTAATcaa                                        >  1:1418544/1‑150 (MQ=255)                                                                                                               ccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATcaa                                        <  1:397738/150‑1 (MQ=255)                                                                                                                          cgcgGAGCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATCAACAAACGCCTgc                             <  1:909912/150‑1 (MQ=255)                                                                                                                              gagcTGGCGCGTTTTCTCGAAAACCCGGCAGATCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATCAACAAACGCCTGCGCCa                         >  1:42674/4‑150 (MQ=38)                                                                                                                                         ttttCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGGCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATCAACAAACGCCTGCGCCACGCGGGGGTTa              <  2:134296/150‑1 (MQ=255)                                                                                                                                          tttCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATCAACAAACGCCTGCGCCACGCGGGGGTTAg             <  1:939398/150‑1 (MQ=255)                                                                                                                                           ttctcgaAAACCCGGCAGATCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATCAACAAACGCCTGCGCCACGCGGGGGTTAGc            <  2:42674/145‑1 (MQ=255)                                                                                                                                               agaAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATCAACAAACGCCTGCGCCACGCGGGGGTTAGCaaaa        >  2:1224569/1‑150 (MQ=255)                                                                                                                                                  aaaCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATCAACAAACGCCTGCGCCACGCGGGGGTTAGCAAAagag     <  1:1224569/150‑1 (MQ=255)                                                                                                                                                     cgcgCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATCAACAAACGCCTGCGCCACGCGGGGGTTAGCAAAAGAGGTa  >  2:243631/1‑150 (MQ=255)                                                                                                                                                     cgcgCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATCAACAAACGCCTGCGCCACGCGGGGGTTAGCAAAAGAGGTa  >  1:976456/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GGCGATCATCATCACCCAGGTTACCGGACCTGCGCGTTCACGCAAGGCGGCAAAGAAAGGATAACGGCGATAATGCAGGCCACTGCCGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGGATAACACCCTGCGTCGCCATGTAATCAACAAACGCCTGCGCCACGCGGGGGTTAGCAAAAGAGGTA  >  NC_000913/3561144‑3561441

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