Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	436,580	G→C	5.4%	H236Q (CAC→CAG)	PP_0358 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	436,580	0	G	C	5.4%	245.1 / 2.8	74	H236Q (CAC→CAG)	PP_0358	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (32/38);  new base C (2/2);  total (34/40)
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.76e-01

CGACGCTGTAAATCGCCAAACCTGCCCAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTGCAGTAGCAGCCAGCCAAGGGCCAGC  >  NC_002947/436432‑436724                                                                                                                                                     |                                                                                                                                                 cGACGCTGTAAATCGCCAAACCTGCCCAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAgg                                                                                                                                                  >  1:202295/1‑149 (MQ=255)     gCTGTAAATCGCCAAACCTGCCCAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGcc                                                                                                                                              <  4:29084/149‑1 (MQ=255)        gTAAATCGCCAAACCTGCCCAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTg                                                                                                                                           >  4:11612/1‑149 (MQ=255)         tAAATCGCCAAACCTGCCCAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTg                                                                                                                                           <  5:61066/148‑1 (MQ=255)           aaTCGCCAAACCTGCCCAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTgcgg                                                                                                                                        <  8:240485/149‑1 (MQ=255)               gCCAAACCTGCCCAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTgcggcggc                                                                                                                                    <  6:365877/149‑1 (MQ=255)                     ccTGCCCAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGaa                                                                                                                              <  6:76445/149‑1 (MQ=255)                      cTGCCCAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAg                                                                                                                             <  5:343000/149‑1 (MQ=255)                         cccAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGc                                                                                                                            <  6:13507/147‑1 (MQ=255)                          ccAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATa                                                                                                                         >  2:152937/1‑149 (MQ=255)                              aTGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAg                                                                                                                     >  8:350937/1‑149 (MQ=255)                              aTGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAg                                                                                                                     >  8:15894/1‑149 (MQ=255)                              aTGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAg                                                                                                                     <  4:323172/149‑1 (MQ=255)                                  aCATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGc                                                                                                                 >  3:255444/1‑149 (MQ=255)                                   cATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCa                                                                                                                >  7:81478/1‑149 (MQ=255)                                    aTGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCAc                                                                                                               <  3:11612/149‑1 (MQ=255)                                      gAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCAc                                                                                                               >  4:397953/1‑147 (MQ=255)                                        aGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGc                                                                                                           <  2:203772/149‑1 (MQ=255)                                         ggCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGc                                                                                                           <  2:304768/148‑1 (MQ=255)                                         ggCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGc                                                                                                           >  1:304768/1‑148 (MQ=255)                                          gccgccAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGACGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGt                                                                                                         >  7:281796/1‑149 (MQ=255)                                            cgccAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAAAAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAACGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCa                                                                                                       >  1:100804/1‑149 (MQ=255)                                             gccAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCAc                                                                                                      <  8:81478/149‑1 (MQ=255)                                             gccAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGCACGGCACCTGCCGtacac                                                                                                                                        >  6:407227/1‑111 (MQ=255)                                             gccAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGCACGGCACCTGCCGtacac                                                                                                                                        <  5:407227/115‑5 (MQ=255)                                             gccAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGCACGGCACCTGCCGtacac                                                                                                                                        <  5:412717/115‑5 (MQ=255)                                             gccAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGCACGGCACCTGCCGtacac                                                                                                                                        >  6:412717/1‑111 (MQ=255)                                                        cgacgaCAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGgccg                                                                                           >  4:392967/1‑149 (MQ=255)                                                         gacgacAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGgccgc                                                                                          >  3:128253/1‑149 (MQ=255)                                                          acgacAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGgccgc                                                                                          >  4:138261/1‑148 (MQ=255)                                                              cAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCa                                                                                     <  6:310083/149‑1 (MQ=255)                                                                     tcgtcgAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAg                                                                              >  8:225209/1‑149 (MQ=255)                                                                        tcgAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAgcgc                                                                           <  4:302279/149‑1 (MQ=255)                                                                          gAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCt                                                                         <  3:397953/149‑1 (MQ=255)                                                                            aCAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAgc                                                                             <  8:12981/143‑1 (MQ=255)                                                                            aCAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAgc                                                                             >  7:12981/1‑143 (MQ=255)                                                                                cACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCa                                                                                     <  6:89802/131‑1 (MQ=255)                                                                                cACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCa                                                                                     >  5:89802/1‑131 (MQ=255)                                                                                  cTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGctgctg                                                                 <  3:138261/149‑1 (MQ=255)                                                                                  cTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGctgctg                                                                 <  1:177061/149‑1 (MQ=255)                                                                                   tACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTgg                                                                >  6:49622/1‑149 (MQ=255)                                                                                     cGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGta                                                              <  7:15894/149‑1 (MQ=255)                                                                                       gCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGtata                                                            <  1:152937/149‑1 (MQ=255)                                                                                       gCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGta                                                              >  1:38346/1‑147 (MQ=255)                                                                                        ccTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGt                                                               >  1:165881/1‑145 (MQ=255)                                                                                        ccTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGt                                                               <  2:165881/145‑1 (MQ=255)                                                                                                    ccAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCa                                                                                                             >  5:273297/1‑87 (MQ=255)                                                                                                    ccAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCa                                                                                                             <  6:273297/87‑1 (MQ=255)                                                                                                    ccAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTgg                                               <  5:107942/149‑1 (MQ=255)                                                                                                     cAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTggg                                              <  2:202295/149‑1 (MQ=255)                                                                                                     cAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTggg                                              <  1:290718/149‑1 (MQ=255)                                                                                                       agtgcGGCGCCAGGTAATTCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGAGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCa                                            <  8:345158/145‑1 (MQ=255)                                                                                                          tCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTggg                                              >  3:32590/1‑144 (MQ=255)                                                                                                          tCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTggg                                              <  4:32590/144‑1 (MQ=255)                                                                                                          tCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTc                                         <  2:38346/149‑1 (MQ=255)                                                                                                             gcgcCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATc                                      <  4:255444/149‑1 (MQ=255)                                                                                                                ccAGGTACTGCAGGAAGCCCAACGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCt                                   <  2:100804/149‑1 (MQ=255)                                                                                                                     tACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCgg                              <  3:392967/149‑1 (MQ=255)                                                                                                                        tGCAGGAAGCCCAGCGTGGTTTACTGCAGGTGCAGTGCGGCGGCGTTGa                                                                                                                               >  6:326796/1‑49 (MQ=255)                                                                                                                        tGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGa                                                                                                                               <  5:326796/49‑1 (MQ=255)                                                                                                                        tGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTg                           >  5:70257/1‑149 (MQ=255)                                                                                                                        tGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTg                           >  7:182823/1‑149 (MQ=255)                                                                                                                          cAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTggg                                              >  1:151260/1‑128 (MQ=255)                                                                                                                          cAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTggg                                              <  2:151260/128‑1 (MQ=255)                                                                                                                           aGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTGCAg                        <  3:203112/149‑1 (MQ=255)                                                                                                                             gAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCTGGTGCAg                        >  4:158324/1‑147 (MQ=255)                                                                                                                                   cAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTGCAGTagcagc                 >  2:356210/1‑148 (MQ=255)                                                                                                                                    aGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTGCAGTAGcagcca               <  8:281796/149‑1 (MQ=255)                                                                                                                                      cGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTGCAGTAGcagccagc             >  6:147261/1‑149 (MQ=255)                                                                                                                                           tgtACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTGCAGTAGCAGCCAGCCAAgg        >  3:27195/1‑149 (MQ=255)                                                                                                                                             tACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCa                                        >  6:21965/1‑115 (MQ=255)                                                                                                                                             tACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTGCAGTAGCAGCCAGCCAAggg       <  5:49622/148‑1 (MQ=255)                                                                                                                                               cGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGTCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTGCAGTAGCAGCCAGCCAAGGGCCa    >  8:365917/1‑149 (MQ=255)                                                                                                                                                ggCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTGCAGTAGCAGCCAGCCAAGGGCCAg   >  3:53000/1‑149 (MQ=255)                                                                                                                                                 gCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTGCAGTAGCAGCCAGCCAAGGGCCAGc  <  7:350937/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CGACGCTGTAAATCGCCAAACCTGCCCAGATGAACATGAAGGCCGCCAGGGTGCTCGACGACAAATGCTCGTCGAACAGCACTACGGCCTGCAACAGCACCAAGGTCGGCGCCAGGTACTGCAGGAAGCCCAGCGTGGTGTACGGCAGGTGCCGTGCGGCGGCGTTGAAGCATACCAGCGGCACCAGCGTCACCGGCCCGGCCGCCATCAGCCACAGCGCCTCGCTGCTGGTATAGAACGCGCCCTGGGCACTCATCGCTGTCGGGTGCAGTAGCAGCCAGCCAAGGGCCAGC  >  NC_002947/436432‑436724

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