Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,209,838	Δ1 bp	5.3%	coding (297/1134 nt)	rnd ←	ribonuclease D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,209,838	0	A	.	5.3%	236.4 / 7.1	75	coding (297/1134 nt)	rnd	ribonuclease D
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base . (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TTTCCGAGAGCGGGCGTTGCAGCCAGTCGGAACGGGTTTCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAGGAAGGCGCGTTGGCC  >  NC_002947/5209690‑5209984                                                                                                                                                     |                                                                                                                                                   tttCCGAGAGCGGGCGTTGCAGCCAGTCGGAACGGGTTTCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGa                                                                                                                                                    <  2:56416/149‑1 (MQ=255)     cGAGAGCGGGCGTTGCAGCCAGTCGGAACGGGTTTCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAgg                                                                                                                                                  >  6:224556/1‑147 (MQ=255)     cGAGAGCGGGCGTTGCAGCCAGTCGGAACGGGTTTCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAgg                                                                                                                                                  >  1:68960/1‑147 (MQ=255)      gagaGCGGGCGTTGCAGCCAGTCGGAACGGGTTTCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAg                                                                                                                                               >  4:97069/1‑149 (MQ=255)                 ttGCAGCCAGTCGGAACGGGTTTCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGt                                                                                                                                      >  4:351771/1‑147 (MQ=255)                       ccAGTGGGTAGGGGTTTCGCCTTTGGGCTATTCGATGCGCAGCACTTCCTGGACCAGGCGCGAATAGCCCATCGAGAAGCCCTGGTTGAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCgca                                                                                                                              <  8:212214/149‑1 (MQ=255)                       ccAGTCGGAACGGGTTTCGCCTTTGGGCAATTCGATGACCAGCACTTCCTGCACCAGGCGCGAATAGCCCATAGAGAAGCCCAGGTTCAGGTAAACGGCGGCCAGCTGAGTGTCTAACAGCGGCTGAGGCAGCTTTCCGGTCAGGCgca                                                                                                                              <  8:206742/149‑1 (MQ=255)                               aaCGGGTTTCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTG‑GGCAGCTTGCCGGTCAGGCGCAGCAGCACtt                                                                                                                     >  1:182664/1‑149 (MQ=255)                                 cGGGTTTCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTc                                                                                                                    <  8:287949/149‑1 (MQ=255)                                    gTTTCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAg                                                                                                                 <  5:56454/149‑1 (MQ=255)                                    gTTTCGCCATTGTGAAATTCCATTCCCAGAACTTCCTGCACCAGGCGCGAAAAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCTAg                                                                                                                 <  4:261563/149‑1 (MQ=255)                                     tttCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAgg                                                                                                                >  8:299255/1‑149 (MQ=255)                                       tCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCg                                                                                                                                        <  8:70119/123‑1 (MQ=255)                                       tCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCAc                                                                                                                       >  4:230600/1‑140 (MQ=255)                                       tCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCAc                                                                                                                       <  3:230600/140‑1 (MQ=255)                                        cGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTc                                                                                                              <  8:273720/148‑1 (MQ=255)                                         gCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTG‑GGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTc                                                                                                           <  2:182664/149‑1 (MQ=255)                                          ccTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTc                                                                                                           >  8:54163/1‑149 (MQ=255)                                           cTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCg                                                                                                          <  7:293272/149‑1 (MQ=255)                                           cTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCg                                                                                                          >  3:155182/1‑149 (MQ=255)                                                          tGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTgcagc                                                                                           <  3:97069/149‑1 (MQ=255)                                                             ccAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCAcc                                                                                        >  7:188857/1‑149 (MQ=255)                                                              cAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCt                                                                                       >  2:296870/1‑149 (MQ=255)                                                               aGCACTTGCTGCACCAGGCGCGAATATCCCATCGAGAAGCCCAGGTTCAGGTAACCGTCGGCCTGCTTAGTGTCGATCAGCGGCTTAGGCAGCTTTACGTTCAGGCGCAGAAGCACTTGTAGGTCTTCGCTTCAGGCGTGCAGCACAtt                                                                                      <  4:203995/149‑1 (MQ=255)                                                               aGCACGTCCTGCACCAGGCGCGAATCGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTTCAGGCGTGCAGCACCtt                                                                                      <  2:68960/149‑1 (MQ=255)                                                                    ttCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCgcagca                                                                                                                           <  1:231017/107‑1 (MQ=255)                                                                    ttCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCgcagca                                                                                                                           >  2:231017/1‑107 (MQ=255)                                                                           aCCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCt                                                                          >  2:226971/1‑149 (MQ=255)                                                                               ggCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGAGTCGAACAGCGGCTGCGGCCGCTTGCCGGTCAGGCGCAGCAGCCCTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACCGCGCTGGCt                                                                      >  6:37302/1‑149 (MQ=255)                                                                                  gcgAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTc                                                                    >  4:353943/1‑148 (MQ=255)                                                                                    gAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCcag                                                                 >  4:54714/1‑149 (MQ=255)                                                                                     aaTAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCcagc                                                                <  5:355798/149‑1 (MQ=255)                                                                                      aTAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCcagca                                                               <  1:226971/149‑1 (MQ=255)                                                                                          cccATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGgcg                                                                                                                                >  2:219191/1‑80 (MQ=255)                                                                                          cccATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGgcg                                                                                                                                <  1:219191/80‑1 (MQ=255)                                                                                          cccATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTc                                                           >  4:93974/1‑149 (MQ=255)                                                                                          cccATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTc                                                           >  2:357433/1‑149 (MQ=255)                                                                                           ccATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCCGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCg                                                          >  6:276157/1‑149 (MQ=255)                                                                                             aTCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGc                                                        >  3:55737/1‑149 (MQ=255)                                                                                                gagaAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCAGTCTTCCAGCAGGTCGGCCaa                                                     <  7:299255/149‑1 (MQ=255)                                                                                                 agaAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAg                                                    >  7:258774/1‑149 (MQ=255)                                                                                                  gaAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAgg                                                   <  4:155182/149‑1 (MQ=255)                                                                                                  gaAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAgg                                                   <  1:296870/149‑1 (MQ=255)                                                                                                     gCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGtt                                                >  6:189777/1‑149 (MQ=255)                                                                                                       ccAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGc                                              >  7:101401/1‑149 (MQ=255)                                                                                                               aGGTAACCGTCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGc                                       <  8:188857/148‑1 (MQ=255)                                                                                                               aGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGcc                                      >  3:322264/1‑149 (MQ=255)                                                                                                               aGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGcc                                      >  8:113459/1‑149 (MQ=255)                                                                                                                ggTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGAc                                                                                   <  4:287148/103‑1 (MQ=255)                                                                                                                ggTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGAc                                                                                   >  3:287148/1‑103 (MQ=255)                                                                                                                ggTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGtt                                                <  5:218199/138‑1 (MQ=255)                                                                                                                ggTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGtt                                                >  6:218199/1‑138 (MQ=255)                                                                                                                ggTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCa                                            >  5:2794/1‑142 (MQ=255)                                                                                                                ggTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCa                                            <  6:2794/142‑1 (MQ=255)                                                                                                                 gTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGa                                    >  2:250496/1‑149 (MQ=255)                                                                                                                      cggcggCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGtt                                         <  5:186867/139‑1 (MQ=255)                                                                                                                      cggcggCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGtt                                         >  6:186867/1‑139 (MQ=255)                                                                                                                      cggcggCCAGCTGAGTGTCGAACAGCGGCT‑‑GGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGtt                                         >  7:48930/1‑137 (MQ=255)                                                                                                                      cggcggCCAGCTGAGTGTCGAACAGCGGCT‑‑GGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGtt                                         <  8:48930/137‑1 (MQ=255)                                                                                                                           gcCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCagcagc                            <  6:183849/147‑1 (MQ=255)                                                                                                                           gcCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCagcagc                            >  5:183849/1‑147 (MQ=255)                                                                                                                             cAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGt                        <  8:101401/149‑1 (MQ=255)                                                                                                                             cAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGt                        <  3:93974/149‑1 (MQ=255)                                                                                                                             cAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGt                        <  1:250496/149‑1 (MQ=255)                                                                                                                                   aGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATc                   <  1:357433/148‑1 (MQ=255)                                                                                                                                   aGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCa                  <  3:351771/149‑1 (MQ=255)                                                                                                                                       tCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCGAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAgg                >  2:341155/1‑147 (MQ=255)                                                                                                                                       tCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAgg                >  3:250520/1‑147 (MQ=255)                                                                                                                                       tCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAGGaa              <  7:113459/149‑1 (MQ=255)                                                                                                                                       tCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAGGa               <  3:54714/148‑1 (MQ=255)                                                                                                                                        cGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAGGAAg             >  3:27677/1‑149 (MQ=255)                                                                                                                                         gAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAGGAAgg            >  3:101381/1‑149 (MQ=255)                                                                                                                                         gAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAGGAAg             <  8:258774/148‑1 (MQ=255)                                                                                                                                           aCAGAGGCTGAGGCAGCGTGACTTTCATTCGCAGCATCACTTCGATGTCTTCTCTGCAGGCGTGAAGAACCTTGACAACAGCGCTGTCTTCCAGAAGGTCGGACAAGGGTTTCCATTTGACGATCAGCAGCGGGTCGATCAGGAAGgc           <  5:37302/148‑1 (MQ=255)                                                                                                                                                ggCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAGGAAGGCGCGTTg     >  2:340299/1‑149 (MQ=255)                                                                                                                                                 gCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAGGAAGGCGCGTTgg    >  1:179833/1‑149 (MQ=255)                                                                                                                                                     aGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAGGAAGGCGCGTTGGcc  <  5:224556/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TTTCCGAGAGCGGGCGTTGCAGCCAGTCGGAACGGGTTTCGCCTTTGGGCAATTCGATGCCCAGCACTTCCTGCACCAGGCGCGAATAGCCCATCGAGAAGCCCAGGTTCAGGTAACCGGCGGCCAGCTGAGTGTCGAACAGCGGCTGAGGCAGCTTGCCGGTCAGGCGCAGCAGCACTTCGAGGTCTTCGCTGCAGGCGTGCAGCACCTTGACCACAGCGCTGTCTTCCAGCAGGTCGGCCAAGGGTTGCCAGTTGCCGATCAGCAGCGGGTCGATCAGGAAGGCGCGTTGGCC  >  NC_002947/5209690‑5209984

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