Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,156,931	G→A	19.4%	P176P (CCC→CCT)	PP_2772 ←	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,156,931	0	G	A	19.4%	162.3 / 30.9	77	P176P (CCC→CCT)	PP_2772	monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base A (7/8);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTACCAACCATCGAAGCGTGTGAGTACGGTGC  >  NC_002947/3156783‑3157073                                                                                                                                                      |                                                                                                                                               tGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTg                                                                                                                                                <  1:357373/149‑1 (MQ=255)        gggggTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATgg                                                                                                                                         <  6:142160/149‑1 (MQ=255)          gggTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTAGGGATGGCa                                                                                                                                       <  1:59306/149‑1 (MQ=255)                tCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGAt                                                                                                                                  >  7:226845/1‑148 (MQ=255)                  aGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTAGGGATGGCAGCGATca                                                                                                                                <  5:352091/148‑1 (MQ=255)                   gCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATcagc                                                                                                                              <  8:226845/149‑1 (MQ=255)                    cGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATcagca                                                                                                                             <  1:109632/149‑1 (MQ=255)                     gggCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATcagcag                                                                                                                            >  2:123363/1‑149 (MQ=255)                       gCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTCTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCAGCTAGGTATGGCACCGATaagaagaa                                                                                                                          >  7:166724/1‑147 (MQ=255)                             gTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTgtcgc                                                                                                                    >  5:117402/1‑149 (MQ=255)                               cGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCa                                                                                                                  <  7:72353/149‑1 (MQ=255)                                  gaAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCa                                                                                                               >  8:170215/1‑149 (MQ=255)                                      cATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACtt                                                                                                            >  7:132397/1‑148 (MQ=255)                                        tctcGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCa                                                                                                         >  6:154683/1‑149 (MQ=255)                                          tcGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCAt                                                                                                        >  6:361687/1‑148 (MQ=255)                                            gTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCAt                                                                                                        <  5:361687/146‑1 (MQ=255)                                            gTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCAt                                                                                                        <  5:225583/146‑1 (MQ=255)                                             ttCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATc                                                                                                       >  6:164937/1‑146 (MQ=255)                                             ttCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATc                                                                                                       <  5:164937/146‑1 (MQ=255)                                                tctTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTg                                                                                                  >  4:148568/1‑148 (MQ=255)                                                  ttCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTg                                                                                                  >  2:96270/1‑146 (MQ=255)                                                  ttCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTg                                                                                                  >  1:55919/1‑146 (MQ=255)                                                  ttCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTg                                                                                                  >  7:214858/1‑146 (MQ=255)                                                     gACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCagag                                                                                             <  8:214858/148‑1 (MQ=255)                                                      acgcgcacgcgcTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAgg                                                                                            >  6:336335/1‑148 (MQ=255)                                                        gcgcacgcgcTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAggtcg                                                                                         >  5:16600/1‑146 (MQ=255)                                                         cgcacgcgcTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAAt                                                                                        <  5:154683/149‑1 (MQ=255)                                                          gcacgcgcTCGAAAGTCGCGCTGCCGCTGTCGGTCTAGGCGCTGGCCCATGTTGTTCCAGGCATCATGCAGTGTGATCCCGGTCCGGCCGCGGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCGGAGGCAAt                                                                                        >  1:241251/1‑149 (MQ=255)                                                              gcgcTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGTGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCATTCAGc                                                                                    <  4:82187/148‑1 (MQ=255)                                                              gcgcTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGTGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCATTCAGc                                                                                    >  3:82187/1‑148 (MQ=255)                                                              gcgcTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCa                                                                                   <  7:170215/149‑1 (MQ=255)                                                                  tCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATgg                                                                                                                                         <  2:279256/91‑1 (MQ=255)                                                                  tCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATgg                                                                                                                                         >  1:279256/1‑91 (MQ=255)                                                                  tCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTc                                                                               >  6:298982/1‑149 (MQ=255)                                                                     aaaGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCATTCAGCATTTCCGa                                                                            <  1:96270/149‑1 (MQ=255)                                                                     aaaGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGa                                                                            <  2:241251/149‑1 (MQ=255)                                                                       aGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCag                                                                                               <  7:137793/128‑1 (MQ=255)                                                                       aGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCag                                                                                               >  8:137793/1‑128 (MQ=255)                                                                                    cTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCt                                                               >  2:29189/1‑147 (MQ=255)                                                                                    cTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGACGCAATCAGCATTTCCGAGTCGCTGGCGCCt                                                               <  1:29189/147‑1 (MQ=255)                                                                                       tCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAg                                                            >  7:180250/1‑147 (MQ=255)                                                                                                  gcTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACtcttc                                                >  7:306394/1‑149 (MQ=255)                                                                                                  gcTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACtcttc                                                >  2:51215/1‑149 (MQ=255)                                                                                                  gcTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACtcttc                                                <  2:55919/149‑1 (MQ=255)                                                                                                    tGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCGGAAACTcttctt                                              <  8:132397/149‑1 (MQ=255)                                                                                                      gCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGTCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGc                                            <  6:117402/149‑1 (MQ=255)                                                                                                         cATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACtt                                                                                                            <  3:169131/82‑1 (MQ=255)                                                                                                         cATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACtt                                                                                                            >  4:169131/1‑82 (MQ=255)                                                                                                             ttgttCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGaaaa                                     >  5:128081/1‑149 (MQ=255)                                                                                                             ttgttCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCAGaaaa                                     >  5:6400/1‑149 (MQ=255)                                                                                                             ttgttCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAgg                                                                                            <  3:70037/94‑1 (MQ=255)                                                                                                             ttgttCCAGTCATCATGCAGTCTGATACCGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAgg                                                                                            >  4:70037/1‑94 (MQ=255)                                                                                                                ttCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTc                                  <  1:123363/149‑1 (MQ=255)                                                                                                                ttCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTc                                  >  4:329657/1‑149 (MQ=255)                                                                                                                  ccAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCg                                 >  7:266056/1‑148 (MQ=255)                                                                                                                  ccAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCg                                 >  8:93755/1‑148 (MQ=255)                                                                                                                      tcatcaTGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGaa                                                       >  5:84749/1‑122 (MQ=255)                                                                                                                      tcatcaTGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGaa                                                       <  6:84749/122‑1 (MQ=255)                                                                                                                           atGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTaccaacca                        <  1:51215/148‑1 (MQ=255)                                                                                                                            tGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGc                                                           >  2:350245/1‑112 (MQ=255)                                                                                                                            tGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGc                                                           <  1:350245/112‑1 (MQ=255)                                                                                                                              cAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTc                                  >  4:124500/1‑135 (MQ=255)                                                                                                                              cAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTc                                  <  3:124500/135‑1 (MQ=255)                                                                                                                               aGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACtcttc                                                <  7:103486/120‑1 (MQ=255)                                                                                                                               aGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACtcttc                                                >  8:103486/1‑120 (MQ=255)                                                                                                                                 tgtgATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTACCAACCATCGaa                   <  4:202294/147‑1 (MQ=255)                                                                                                                                     aTACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGGGCCTGAGCAGAAACTCTTCTTGCCGGAAAAc                                    >  2:39689/1‑126 (MQ=255)                                                                                                                                     aTACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGGGCCTGAGCAGAAACTCTTCTTGCCGGAAAAc                                    <  1:39689/126‑1 (MQ=255)                                                                                                                                      tACCGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTACCAACCATCGAAGCgtgt             >  3:20016/1‑148 (MQ=255)                                                                                                                                      tACCGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTACCAACCATCGAAGCgtgt             <  5:336335/148‑1 (MQ=255)                                                                                                                                       aCCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTacc                             <  2:199377/131‑1 (MQ=255)                                                                                                                                       aCCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTacc                             >  1:199377/1‑131 (MQ=255)                                                                                                                                        ccGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTACCAACCATCGAAGCGTGTGa           <  8:306394/148‑1 (MQ=255)                                                                                                                                        ccGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTACCAACCATCGAAGCGTGTGAg          >  7:202614/1‑149 (MQ=255)                                                                                                                                        ccGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTACCAACCATCGAAGCGTGTGa           <  4:20016/148‑1 (MQ=255)                                                                                                                                        ccGGTCCGGCCGCTAGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTACCAACCATCGAAGCGTGTGa           <  5:298982/148‑1 (MQ=255)                                                                                                                                            tCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTACCAACCATCGAAGCGTGTGAGTACg      <  8:266056/149‑1 (MQ=255)                                                                                                                                                gCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTACCAACCATCGAAGCGTGTGAGTACGGTGc  >  8:161858/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                               TGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCT‑GGCGCTGGCCCATGTTGTTCCAGTCATCATGCAGTGTGATACCGGTCCGGCCGCTGGGGATGGCAGCGATCAGCAGCTTGCCGCCAGCACTTTCATCGACTGCAGAGGCAATCAGCATTTCCGAGTCGCTGGCGCCTGAGCAGAAACTCTTCTTGCCGGAAAACTCGTACCAACCATCGAAGCGTGTGAGTACGGTGC  >  NC_002947/3156783‑3157073

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