Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,734,244	G→A	6.3%	P474S (CCG→TCG)	hutH ←	histidine ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,734,244	0	G	A	6.3%	188.5 / 2.1	63	P474S (CCG→TCG)	hutH	histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base A (2/2);  total (33/30)
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.81e-01

AAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGTTTTCGGCCATTTCCCACAGGCGCTTGCCG  >  NC_002947/5734104‑5734392                                                                                                                                             |                                                                                                                                                     aaGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCAc                                                                                                                                              >  1:155591/1‑149 (MQ=255)  aGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACt                                                                                                                                             >  5:214132/1‑149 (MQ=255)      ccgccCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTgcg                                                                                                                                          >  6:66297/1‑148 (MQ=255)          ccAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAgcg                                                                                                                                     >  7:144212/1‑149 (MQ=255)          ccAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAgcg                                                                                                                                     >  6:159694/1‑149 (MQ=255)                  cccAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGgcgc                                                                                                                             >  5:118318/1‑149 (MQ=255)                   ccAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGgcgcg                                                                                                                            >  1:167738/1‑149 (MQ=255)                     gcgggCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGcc                                                                                                                          <  6:266281/147‑1 (MQ=255)                      ggggCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCt                                                                                                                         >  2:219218/1‑149 (MQ=255)                              aCAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAg                                                                                                                  <  8:144212/148‑1 (MQ=255)                                aGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTc                                                                                                                     >  4:136305/1‑143 (MQ=255)                                aGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTc                                                                                                                     <  3:136305/143‑1 (MQ=255)                                  gCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTgg                                                                                                             >  6:154673/1‑149 (MQ=255)                                  gCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGACTGGCGCGCCTGCTCCAGCTTgg                                                                                                             >  8:154781/1‑149 (MQ=255)                                   cTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTgg                                                                                                             <  5:154673/148‑1 (MQ=255)                                       gCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGaa                                                                                                        <  8:110091/149‑1 (MQ=255)                                       gCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGACACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCg                                                                                                          >  4:41379/1‑147 (MQ=255)                                       gCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGACACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCg                                                                                                          >  4:45553/1‑147 (MQ=255)                                             cacCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCtt                                                                                                   <  5:161822/148‑1 (MQ=255)                                              acCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGACACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCa                                                                                                 <  3:41379/149‑1 (MQ=255)                                              acCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGACACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCa                                                                                                 <  3:45553/149‑1 (MQ=255)                                               cccGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAg                                                                                                <  8:212703/149‑1 (MQ=255)                                                ccgccgGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAgg                                                                                               >  7:154474/1‑149 (MQ=255)                                                         gcaAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTgcg                                                                                       >  1:35176/1‑148 (MQ=255)                                                           aaGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAg                                                                                    <  7:261865/149‑1 (MQ=255)                                                             gCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAg                                                                                    <  6:233367/147‑1 (MQ=255)                                                             gCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAg                                                                                    >  5:233367/1‑147 (MQ=255)                                                               cGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTcc                                                                                <  5:159694/149‑1 (MQ=255)                                                                ggTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTc                                                                                 >  2:61917/1‑147 (MQ=255)                                                                 gTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAg                                                                              <  6:35299/149‑1 (MQ=255)                                                                    aGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGc                                                                            <  2:266784/148‑1 (MQ=255)                                                                     gACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGccc                                                                          <  6:194388/149‑1 (MQ=255)                                                                       cTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTg                                                                        <  5:243777/149‑1 (MQ=255)                                                                          cccTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTgg                                                                       <  2:155591/147‑1 (MQ=255)                                                                           ccTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTTAGGCCTTTGCGCAGGTCCAGGCCCTGGCAc                                                                    >  2:178022/1‑149 (MQ=255)                                                                           ccTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCAc                                                                    >  2:18181/1‑149 (MQ=255)                                                                           ccTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCAc                                                                    <  1:219218/149‑1 (MQ=255)                                                                                gCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACgcgc                                                                <  2:167738/148‑1 (MQ=255)                                                                                  cAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGcc                                                               >  1:242419/1‑147 (MQ=255)                                                                                  cAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCa                                                             <  5:87693/149‑1 (MQ=255)                                                                                   aaCAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAg                                                            >  6:108832/1‑149 (MQ=255)                                                                                        gTCCACGGACTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGAGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTACAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGTCCCTGGCACGCGCCCAGCCAt                                                        >  8:120312/1‑148 (MQ=255)                                                                                           cACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCACTGGCACGCGCCCAGCCATTCGa                                                    >  6:275881/1‑149 (MQ=255)                                                                                           cACGGCCTTTTCGATATCAGGAGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGa                                                    >  6:99941/1‑149 (MQ=255)                                                                                            aCGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGAt                                                   <  8:255002/149‑1 (MQ=255)                                                                                            aCGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGAt                                                   <  5:108832/149‑1 (MQ=255)                                                                                             cGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGACACGCGCCCAGCCATTCGATg                                                  <  1:61917/149‑1 (MQ=255)                                                                                                   tttCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGCGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAg                                             >  2:163863/1‑148 (MQ=255)                                                                                                    ttCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGc                                            <  1:286355/148‑1 (MQ=255)                                                                                                       gATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCCGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGc                                        >  4:23524/1‑149 (MQ=255)                                                                                                             cGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACggg                                  >  1:122863/1‑149 (MQ=255)                                                                                                                   gaagaaACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTTAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGtttt                            <  1:178022/149‑1 (MQ=255)                                                                                                                   gaagaaACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGtttt                            >  4:179934/1‑149 (MQ=255)                                                                                                                   gaagaaACCGTCACGGTCGAAGTGCGGCACCTCACTGCGCAGCGCATGGCGCGCCTGCTCCAGAATGGCCGAAGTCTTCAGCCCTTTCCGCAGGTCCACGCCCTGGCACGATCCCAGCCATTCGATGGCCAGCACTCCACTGGTAtttt                            >  8:288811/1‑149 (MQ=255)                                                                                                                        aaCGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGTTTTCGGc                        <  2:242419/148‑1 (MQ=255)                                                                                                                          cggtcacggtcGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGTTTTCGGCCAt                     >  3:265820/1‑149 (MQ=255)                                                                                                                           ggtcacggtcGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCCGCCATTCGATGGCCAGCACGCCACGGGTGTTTTCGGCCAtt                    >  5:31609/1‑149 (MQ=255)                                                                                                                             tcacggtcGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGTTTTCGGCCATTTc                  <  7:146377/149‑1 (MQ=255)                                                                                                                             tcacggtcGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGTTTTCGGCCATTTc                  <  4:265820/149‑1 (MQ=255)                                                                                                                               acggtcGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGTTTTCGGCCATTTcc                 >  6:62411/1‑148 (MQ=255)                                                                                                                                 ggtcGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGTTTTCGGCCATTTCCca               >  5:184147/1‑148 (MQ=255)                                                                                                                                   tcGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCCGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGTTTTCGGCCATTTCCCACAg            <  3:23524/149‑1 (MQ=255)                                                                                                                                        gTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGTGTGTTTTCGGCCATTTCCCACAGGCGCt       <  5:62411/149‑1 (MQ=255)                                                                                                                                             gCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGTTTTCGGCCATTTCCCACAGGCGCTTGCcg  >  4:174495/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                     AAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCACGGTCGTAGTGCGGCACCTCACTGCGCAGCGCCTGGCGCGCCTGCTCCAGCTTGGCCGAAGTCTTCAGGCCTTTGCGCAGGTCCAGGCCCTGGCACGCGCCCAGCCATTCGATGGCCAGCACGCCACGGGTGTTTTCGGCCATTTCCCACAGGCGCTTGCCG  >  NC_002947/5734104‑5734392

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