Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,498,721	G→A	R298R (CGG→CGA)	cysG →	siroheme synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,498,721	0	G	A	98.9%	305.4 / ‑4.9	87	R298R (CGG→CGA)	cysG	siroheme synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (43/43);  minor base T (1/0);  total (44/43)
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.98e-01

ACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTC  >  NC_000913/3498585‑3498846                                                                                                                                         |                                                                                                                              aCGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTg                                                                                                                            >  1:1092913/1‑140 (MQ=255)   gATATTATGAATCTGGAACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGaa                                                                                                                          <  2:1123622/140‑1 (MQ=255)     tattatGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAg                                                                                                                        <  1:14881/140‑1 (MQ=255)        tatGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTg                                                                                                                     <  2:252560/140‑1 (MQ=255)             aTCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGAt                                                                                                                >  1:1033622/1‑140 (MQ=255)               cTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATcc                                                                                                              >  1:845505/1‑140 (MQ=255)               cTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATcc                                                                                                              <  2:487948/140‑1 (MQ=255)                 ggTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGt                                                                                                            >  1:568750/1‑140 (MQ=255)                  gTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGtt                                                                                                           >  2:169565/1‑140 (MQ=255)                      gccgcGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTAtt                                                                                                       >  2:504603/1‑140 (MQ=255)                         gcgATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTAttttt                                                                                                    >  1:239724/1‑140 (MQ=255)                                  aCCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCgtggt                                                                                           >  1:499308/1‑140 (MQ=255)                                     gtgtTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGc                                                                                        >  2:448770/1‑140 (MQ=255)                                          ttCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAaga                                                                                   <  2:239724/140‑1 (MQ=255)                                          ttCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAaga                                                                                   <  2:95204/140‑1 (MQ=255)                                          ttCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAaga                                                                                   <  2:1033622/140‑1 (MQ=255)                                          ttCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAaga                                                                                   <  2:1092913/140‑1 (MQ=255)                                            cgtcgGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGc                                                                                 <  1:1089093/140‑1 (MQ=255)                                                ggCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGa                                                                             >  2:340018/1‑140 (MQ=255)                                                   aaaCGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAac                                                                          <  2:1175729/140‑1 (MQ=255)                                                         gcgGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACtgtg                                                                    <  2:1409183/140‑1 (MQ=255)                                                            ggATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCaa                                                                 <  2:640739/140‑1 (MQ=255)                                                            ggATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCaa                                                                 <  2:499308/140‑1 (MQ=255)                                                             gATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAAc                                                                <  2:766464/140‑1 (MQ=255)                                                              aTACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAAcg                                                               <  1:219202/140‑1 (MQ=255)                                                                  cACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCggg                                                           <  2:29426/140‑1 (MQ=255)                                                                       cGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCggg                                                           >  2:1114855/1‑135 (MQ=255)                                                                       cGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCggg                                                           <  1:1114855/135‑1 (MQ=255)                                                                       cGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTc                                                      >  1:888911/1‑140 (MQ=255)                                                                       cGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTc                                                      >  1:76986/1‑140 (MQ=255)                                                                       cGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTc                                                      >  1:649841/1‑140 (MQ=255)                                                                       cGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTc                                                      >  1:1306722/1‑140 (MQ=255)                                                                       cGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTc                                                      >  1:1254922/1‑140 (MQ=255)                                                                       cGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTc                                                      >  2:570296/1‑140 (MQ=255)                                                                       cGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTc                                                      >  1:375644/1‑140 (MQ=255)                                                                       cGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTc                                                      >  2:550310/1‑140 (MQ=255)                                                                        gTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTcc                                                     <  2:1415417/140‑1 (MQ=255)                                                                         tACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAAcgc                                                              <  1:462262/130‑1 (MQ=255)                                                                         tACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAAcgc                                                              >  2:462262/1‑130 (MQ=255)                                                                         tACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCg                                                    >  1:1149264/1‑140 (MQ=255)                                                                           ccccAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGtt                                                  <  1:121070/140‑1 (MQ=255)                                                                                 gAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCggt                                            >  1:375907/1‑140 (MQ=255)                                                                                 gAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCggt                                            >  1:663266/1‑140 (MQ=255)                                                                                     agaTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGtt                                        <  1:448770/140‑1 (MQ=255)                                                                                     agaTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGtt                                        <  2:1226826/140‑1 (MQ=255)                                                                                     agaTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGtt                                        <  2:888911/140‑1 (MQ=255)                                                                                      gaTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGtt                                        <  2:490130/139‑1 (MQ=255)                                                                                      gaTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGtt                                        >  1:490130/1‑139 (MQ=255)                                                                                      gaTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTc                                       >  1:922892/1‑140 (MQ=255)                                                                                      gaTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTc                                       >  2:467701/1‑140 (MQ=255)                                                                                      gaTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTc                                       >  2:48600/1‑140 (MQ=255)                                                                                      gaTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTc                                       >  1:1075545/1‑140 (MQ=255)                                                                                       aTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTc                                       <  1:318973/139‑1 (MQ=255)                                                                                       aTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTc                                       >  2:318973/1‑139 (MQ=255)                                                                                        ttAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGGGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCg                                     <  1:169565/140‑1 (MQ=255)                                                                                               gATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTa                              >  2:1060731/1‑140 (MQ=255)                                                                                               gATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTa                              >  1:77565/1‑140 (MQ=255)                                                                                               gATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGGGGTTCCGGGTATTa                              >  2:808940/1‑140 (MQ=255)                                                                                                aTCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTAc                             >  2:376115/1‑140 (MQ=255)                                                                                                aTCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTAc                             >  1:678269/1‑140 (MQ=255)                                                                                                aTCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTAc                             >  2:219342/1‑140 (MQ=255)                                                                                                 tCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGTAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTAcc                            >  1:1003248/1‑140 (MQ=255)                                                                                                 tCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGtt                                        >  1:1031842/1‑128 (MQ=255)                                                                                                 tCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGtt                                        <  2:1031842/128‑1 (MQ=255)                                                                                                  ccTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGTGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCg                           >  2:1361804/1‑140 (MQ=255)                                                                                                  ccTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCg                           <  2:568750/140‑1 (MQ=255)                                                                                                  ccTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCg                           <  1:467701/140‑1 (MQ=255)                                                                                                     gctgcGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAg                        <  1:1060731/140‑1 (MQ=255)                                                                                                     gctgcGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAg                        >  2:165446/1‑140 (MQ=255)                                                                                                        gcGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCtt                     <  2:375644/140‑1 (MQ=255)                                                                                                          gggAAGCGCAAAAAGGCAAACGCGTGGTGCGTCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCt                   >  2:907884/1‑140 (MQ=255)                                                                                                          gggAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCt                   >  1:1280892/1‑140 (MQ=255)                                                                                                          gggAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGAATTACCGCAGCTTCt                   >  2:35019/1‑140 (MQ=255)                                                                                                           ggAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTg                  <  1:165446/140‑1 (MQ=255)                                                                                                            gAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGGAAGGTGGCGATCCGTTTATTTTTGGACGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTgg                 <  1:241015/140‑1 (MQ=255)                                                                                                            gAAGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTgg                 <  2:1254922/140‑1 (MQ=255)                                                                                                              aGCGCAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGtt               >  1:418481/1‑140 (MQ=255)                                                                                                               gcgcAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGTAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTg              <  2:1003248/140‑1 (MQ=255)                                                                                                               gcgcAAAAAGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTg              <  2:375907/140‑1 (MQ=255)                                                                                                                 gcAAAAGGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGct            <  1:703690/140‑1 (MQ=255)                                                                                                                   aaaacGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGctct          <  2:1306722/140‑1 (MQ=255)                                                                                                                   aaaaaGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTACTACCGCAGCTTCTGGTTGctct          <  1:480179/140‑1 (MQ=255)                                                                                                                    aaaaGGCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTg         <  1:1395394/140‑1 (MQ=255)                                                                                                                        ggCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTa     <  1:340018/140‑1 (MQ=255)                                                                                                                        ggCAAACGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTa     <  2:1075545/140‑1 (MQ=255)                                                                                                                           aaaCGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGTCTATTc  >  1:449189/1‑140 (MQ=255)                                                                                                                           aaaCGCGTGGTGCGACTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTc  <  1:808940/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                              ACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTC  >  NC_000913/3498585‑3498846

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑