Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,096,065	A→G	5.3%	T133A (ACC→GCC)	PP_5607 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,096,065	0	A	G	5.3%	224.3 / 2.5	75	T133A (ACC→GCC)	PP_5607	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (45/26);  new base G (2/2);  total (47/28)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

TGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCGGCCT  >  NC_002947/4095926‑4096213                                                                                                                                            |                                                                                                                                                     tGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGc                                                                                                                                               >  6:253808/1‑147 (MQ=255)      gCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGcc                                                                                                                                        >  5:111183/1‑149 (MQ=255)       ccTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCCCCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGccc                                                                                                                                       >  5:198656/1‑149 (MQ=255)         tGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTa                                                                                                                                     >  2:118669/1‑149 (MQ=255)                 ctgcATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTg                                                                                                                             <  1:237624/149‑1 (MQ=255)                    catggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTgg                                                                                                                            <  7:246035/147‑1 (MQ=255)                    catggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTgg                                                                                                                            <  2:222745/147‑1 (MQ=255)                    catggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGc                                                                                                                           >  2:272469/1‑148 (MQ=255)                     atggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGc                                                                                                                           >  2:141394/1‑147 (MQ=255)                     atggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTa                                                                                                                         <  7:131040/149‑1 (MQ=255)                      tggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAg                                                                                                                        <  3:49704/149‑1 (MQ=255)                         catggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTATGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGc                                                                                                                      >  6:303769/1‑148 (MQ=255)                           tggGCAGTTTGGTCACCGGCTGGTTCGCCTGCGACCAGGGGTGTGCACCGGTACAGACATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAAc                                                                                                                   <  6:198656/149‑1 (MQ=255)                                  tttGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTg                                                                                                            <  6:317605/149‑1 (MQ=255)                                      ggCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCgg                                                                                                        <  1:141394/149‑1 (MQ=255)                                               tGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGc                                                                                               >  3:343619/1‑149 (MQ=255)                                                ggTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGcc                                                                                              >  8:28726/1‑149 (MQ=255)                                                ggTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGcc                                                                                              >  3:5981/1‑149 (MQ=255)                                                     cccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTAt                                                                                                                                    >  3:60942/1‑106 (MQ=255)                                                     cccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTAt                                                                                                                                    <  4:60942/106‑1 (MQ=255)                                                     cccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCACGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGgc                                                                                         <  2:165439/149‑1 (MQ=255)                                                      ccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcgagag                                                                                                                                      <  3:72448/103‑5 (MQ=255)                                                      ccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcgagag                                                                                                                                      >  4:72448/1‑99 (MQ=255)                                                      ccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGgcg                                                                                        >  1:291767/1‑149 (MQ=255)                                                       cTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGgcgc                                                                                       >  7:87723/1‑149 (MQ=255)                                                       cTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGgcgc                                                                                       >  6:62756/1‑149 (MQ=255)                                                           gACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGtt                                                                                   >  7:291382/1‑149 (MQ=255)                                                                   gtgtGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTgg                                                                           <  1:272469/149‑1 (MQ=255)                                                                              tACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCgtggtg                                                                <  2:291767/149‑1 (MQ=255)                                                                                 aGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTg                                                             <  7:281885/149‑1 (MQ=255)                                                                                 aGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTg                                                             >  8:343476/1‑149 (MQ=255)                                                                                  gCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTgg                                                            <  7:28726/149‑1 (MQ=255)                                                                                     aTCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTACCCCGCCCTGGCGCTGTTTTc                                                                                >  7:183988/1‑126 (MQ=255)                                                                                         aCCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGatcat                                                     >  8:330177/1‑149 (MQ=255)                                                                                         aCCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGatcat                                                     >  6:107893/1‑149 (MQ=255)                                                                                           cTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGAtcatca                                                   >  4:369410/1‑149 (MQ=255)                                                                                           cTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGAtcatca                                                   >  7:93917/1‑149 (MQ=255)                                                                                           cTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGAtcatca                                                   >  6:152161/1‑149 (MQ=255)                                                                                            tCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcgaga                                                                                                                                       >  2:8659/1‑61 (MQ=37)                                                                                            tCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcgaga                                                                                                                                       <  1:8659/64‑4 (MQ=37)                                                                                              gcagcaACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCgtg                                                                   >  4:6694/1‑130 (MQ=255)                                                                                              gcagcaACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCgtg                                                                   <  3:6694/130‑1 (MQ=255)                                                                                              gcagcaACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCAc                                                  <  5:303769/147‑1 (MQ=255)                                                                                              gcagcaACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCg                                                >  2:26289/1‑149 (MQ=255)                                                                                               cagcaACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGt                                               >  2:276254/1‑149 (MQ=255)                                                                                                agcaACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTc                                              <  1:26289/149‑1 (MQ=255)                                                                                                 gCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCG                                             >  8:77223‑M1/1‑149 (MQ=255)                                                                                                  cAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGC                                            <  3:173767‑M1/149‑1 (MQ=255)                                                                                                        cATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGAT                                      >  6:174527‑M1/1‑149 (MQ=255)                                                                                                        cATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTG                                        <  7:330177‑M1/147‑1 (MQ=255)                                                                                                         aTAACCTGTTTGGTCTGCTGATGCTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATg                                     >  6:8579/1‑149 (MQ=255)                                                                                                          tAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATgg                                    <  5:62756/149‑1 (MQ=255)                                                                                                              cTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGcc                                 >  4:81726/1‑148 (MQ=255)                                                                                                                gTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAg                              >  2:181983/1‑149 (MQ=255)                                                                                                                 tttGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgc                             >  1:354932/1‑149 (MQ=255)                                                                                                                 tttGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgc                             >  8:317803/1‑149 (MQ=255)                                                                                                                     gTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcggcg                         >  6:187247/1‑149 (MQ=255)                                                                                                                     gTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcggcg                         >  2:109247/1‑149 (MQ=255)                                                                                                                     gTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcggcg                         >  7:166419/1‑149 (MQ=255)                                                                                                                     gTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcggc                          <  7:77223/148‑1 (MQ=255)                                                                                                                       ctgctgATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcggc                          <  4:195771/146‑1 (MQ=255)                                                                                                                       ctgctgATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcggc                          >  3:195771/1‑146 (MQ=255)                                                                                                                          ctgATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCagag                    >  8:206999/1‑149 (MQ=255)                                                                                                                          ctgATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCa                       <  4:343619/146‑1 (MQ=255)                                                                                                                           tgATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCACATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGt                   >  8:245538/1‑149 (MQ=255)                                                                                                                              tGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGc                <  3:81726/149‑1 (MQ=255)                                                                                                                                  ccTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGCCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAg             <  1:118669/148‑1 (MQ=255)                                                                                                                                    tctcGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTg          >  3:205203/1‑149 (MQ=255)                                                                                                                                    tctcGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTg          >  8:191682/1‑149 (MQ=255)                                                                                                                                    tctcGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTg          >  1:132510/1‑149 (MQ=255)                                                                                                                                        gCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATgg                                    >  5:153801/1‑119 (MQ=255)                                                                                                                                        gCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATgg                                    <  6:153801/119‑1 (MQ=255)                                                                                                                                        gCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCg      >  3:110302/1‑149 (MQ=255)                                                                                                                                        gCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCg      >  7:100713/1‑149 (MQ=255)                                                                                                                                         cTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCgg     >  6:200851/1‑149 (MQ=255)                                                                                                                                            aCCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCGGCCt  <  7:191682/149‑1 (MQ=255)                                                                                                                                            aCCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCGGCCt  <  5:152161/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                     TGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCGGCCT  >  NC_002947/4095926‑4096213

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