Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	996,022	T→A	8.5%	Q525L (CAG→CTG)	PP_0860 ←	sulfite reductase flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	996,022	0	T	A	8.5%	172.7 / 5.9	59	Q525L (CAG→CTG)	PP_0860	sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (30/24);  new base A (3/2);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

TCTGCAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGG  >  NC_002947/995875‑996168                                                                                                                                                     |                                                                                                                                                   tCTGCAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCctgctg                                                                                                                                                   >  4:203394/1‑149 (MQ=255)     cAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAAt                                                                                                                                               >  5:255468/1‑149 (MQ=255)     cAGCCTGCAGGCCAAGCAGGTACACCGGGTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGCCCAGTTGGCAAAGGGTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCtgctgc                                                                                                                                                  <  5:43255/146‑1 (MQ=255)       gCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCctgct                                                                                                                                                    <  6:89440/142‑1 (MQ=255)       gCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCctgct                                                                                                                                                    >  5:89440/1‑142 (MQ=255)          tGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATtgctg                                                                                                                                          >  5:222182/1‑149 (MQ=255)            cAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATtgctg                                                                                                                                          <  3:16021/147‑1 (MQ=255)               gCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAggg                                                                                                                                     >  1:233601/1‑149 (MQ=255)               gCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAATGGCAGCAgga                                                                                                                                     <  5:82742/149‑2 (MQ=255)               gCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAATGGCAGCAgga                                                                                                                                     <  6:255468/149‑2 (MQ=255)                 cAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGggcg                                                                                                                                   <  6:222182/149‑1 (MQ=255)                  aaGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGggcgg                                                                                                                                  >  6:176363/1‑149 (MQ=255)                  aaGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGggcgg                                                                                                                                  >  3:86312/1‑149 (MQ=255)                          acacCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCg                                                                                                                          >  3:20336/1‑149 (MQ=255)                            acCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGc                                                                                                                        >  5:343948/1‑149 (MQ=255)                            acCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGc                                                                                                                        <  7:124273/149‑1 (MQ=255)                               ggTTGGCCTTGGCTGCCGTGGTTGAGCAGTTCGCGGCGCGCAAGAGACCAGTGGGCAAAGGTTGGCGGCTCGC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGGCATTGATGCAGGGCGGCCTGGTCGGAGtt                                                                                                                     >  6:86980/1‑149 (MQ=255)                                gTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTg                                                                                                                    <  2:233601/149‑1 (MQ=255)                                gTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTg                                                                                                                    >  3:171626/1‑149 (MQ=255)                                    gCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTcc                                                                                                                 <  2:308603/148‑1 (MQ=255)                                     ccTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCAc                                                                                                               >  3:116223/1‑149 (MQ=255)                                      cTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACt                                                                                                              >  4:39143/1‑149 (MQ=255)                                                   gTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGa                                                                                                 <  8:8904/149‑1 (MQ=255)                                                       aGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATg                                                                                              <  6:343948/148‑1 (MQ=255)                                                           gTTCGCGGCCCACCAGTGCCCAGTTGGCAAAGGTTGGCGGCTGCC‑AAGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGCATGCActg                                                                                         >  6:282471/1‑149 (MQ=255)                                                            ttCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCGAGGCGCCCACTGGCTGGACACCGGTAATTTGCCCCAGATCCTCCTCCCAt                                                                                                                                               >  1:76069/1‑95 (MQ=255)                                                                   cgcACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCAcc                                                                                 <  3:203394/149‑1 (MQ=255)                                                                                  ttGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCa                                                                  >  5:230928/1‑149 (MQ=255)                                                                                  ttGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCa                                                                  >  8:328333/1‑149 (MQ=255)                                                                                     gCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATg                                                                <  6:230928/148‑1 (MQ=255)                                                                                     gCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGc                                                               >  3:248502/1‑149 (MQ=255)                                                                                       aaaGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCt                                                             >  7:62725/1‑149 (MQ=255)                                                                                       aaaGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCt                                                             <  3:39143/149‑1 (MQ=255)                                                                                          ggTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGaa                                                          <  8:62725/149‑1 (MQ=255)                                                                                                cggcTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGcgac                                                    <  5:176363/149‑1 (MQ=255)                                                                                                 ggcTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGcgac                                                    >  8:95857/1‑148 (MQ=255)                                                                                                  gcTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGt                                                                                                                   >  8:235070/1‑84 (MQ=255)                                                                                                  gcTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGt                                                                                                                   <  7:235070/84‑1 (MQ=255)                                                                                                     gCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGc                                                >  5:68913/1‑148 (MQ=255)                                                                                                      cc‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGa                                              <  7:328333/149‑1 (MQ=255)                                                                                                         aGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGaaa                                            <  4:248502/149‑1 (MQ=255)                                                                                                            cggcCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGc                                                >  5:261506/1‑142 (MQ=255)                                                                                                            cggcCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGc                                                <  6:261506/142‑1 (MQ=255)                                                                                                                cACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAg                                      >  1:206737/1‑148 (MQ=255)                                                                                                                     gctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTg                                 <  7:263443/148‑1 (MQ=255)                                                                                                                         gACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCg                                         <  4:172168/136‑1 (MQ=255)                                                                                                                         gACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCg                                         >  3:172168/1‑136 (MQ=255)                                                                                                                         gACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCggg                             >  3:99665/1‑148 (MQ=255)                                                                                                                         gACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCgg                              <  2:206737/147‑1 (MQ=255)                                                                                                                            acCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACt                         >  2:223565/1‑149 (MQ=255)                                                                                                                              cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCgg                              <  4:51161/142‑1 (MQ=255)                                                                                                                              cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCgg                              >  3:51161/1‑142 (MQ=255)                                                                                                                              cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGa                       <  6:68913/149‑1 (MQ=255)                                                                                                                              cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGa                       <  5:282471/149‑1 (MQ=255)                                                                                                                              cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTAGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGa                       >  8:12555/1‑149 (MQ=255)                                                                                                                                     ttGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGc                                                                                             <  1:222235/72‑1 (MQ=255)                                                                                                                                     ttGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGc                                                                                             >  2:222235/1‑72 (MQ=255)                                                                                                                                     ttGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCAc                >  7:228805/1‑149 (MQ=255)                                                                                                                                        cgcCAGCTCCTTCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCATGAAGGCGACTAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGa             <  5:86980/149‑1 (MQ=255)                                                                                                                                           cAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAggg          >  4:223893/1‑149 (MQ=255)                                                                                                                                           cAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAggg          >  3:350535/1‑149 (MQ=255)                                                                                                                                              gcccTGCAGCAATGGCAGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGc         >  5:62193/3‑147 (MQ=255)                                                                                                                                                ccTGCAGCAATGGCAGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAg      >  6:241068/1‑148 (MQ=255)                                                                                                                                                 ctgctgCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGc         <  7:305401/144‑1 (MQ=255)                                                                                                                                                 ctgctgCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGc         >  8:305401/1‑144 (MQ=255)                                                                                                                                                 ctgctgCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAg      >  6:82758/1‑147 (MQ=255)                                                                                                                                                    ctgcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAgg  <  8:228805/148‑1 (MQ=255)                                                                                                                                                    ctgcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAg   <  4:99665/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TCTGCAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGG  >  NC_002947/995875‑996168

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