Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,155,091	A→T	7.8%	E87V (GAA→GTA)	PP_1910 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,155,091	0	A	T	7.8%	188.5 / 5.1	64	E87V (GAA→GTA)	PP_1910	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/25);  new base T (3/2);  total (37/27)
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

GACAATGTCGGTACGGTGCAGGCGAAGTTCGATTTTGAACGAGATGAACAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGCCAATCGTGCCTGCTCA  >  NC_002947/2154949‑2155238                                                                                                                                               |                                                                                                                                                    gACAATGTCGGTACGGTGCAGGCGAAGTTCGATTTTGAACGAGATGAACAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTa                                                                                                                                                >  1:96151/1‑148 (MQ=255)   cattGTCGGTACGGTGCAGGCGAAGTTCGATTTTGAACGAGATGAACAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACAt                                                                                                                                             <  7:122500/146‑1 (MQ=255)     aTGTCGGTACGGTGCAGGCGAAGTTCGATTTTGAACGAGATGAACAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCAGCCTACCGATCCACAGCGAATGTACATAc                                                                                                                                           >  6:349097/1‑149 (MQ=255)          ggTACGGTGCAGGCGAAGTTCGATTTTGAACGAGATGAACAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACGGCGTATGTACATTCGCTgt                                                                                                                                      <  8:151874/149‑1 (MQ=255)          ggTACGGTGCAGGCGAAGTTCGATTTTGAACGAGATGAACAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACGGCGTATGTACATTCGCTgt                                                                                                                                      >  7:151874/1‑149 (MQ=255)          ggTACGGTGCAGGCGAAGTTCGATTTTGAACGAGATGAACAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCgt                                                                                                                                      >  5:111137/1‑149 (MQ=255)                      gCGAAGTTCGATTTTGAACGAGATGAACAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAaggagtg                                                                                                                          <  7:328150/149‑3 (MQ=255)                      gCGAAGTTCGATTTTGAACGAGATGAACAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAggg                                                                                                                          >  2:344205/1‑149 (MQ=255)                             tCGATTTTGAACGAGATGAACAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGaa                                                                                                                    >  7:248457/1‑148 (MQ=255)                                                aCAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAg                                                                                                >  8:191391/1‑149 (MQ=255)                                                  aGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGc                                                                                              <  2:96151/149‑1 (MQ=255)                                                       gtggtgGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGa                                                                                         >  6:165778/1‑149 (MQ=255)                                                       gtggtgGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCATACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGa                                                                                         >  2:194770/1‑149 (MQ=255)                                                         ggtggtTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAg                                                                                                             <  8:114980/127‑1 (MQ=255)                                                         ggtggtTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAg                                                                                                             >  7:114980/1‑127 (MQ=255)                                                          gtggtTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGa                                                                                         >  5:288585/1‑146 (MQ=255)                                                          gtggtTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGa                                                                                         <  6:288585/146‑1 (MQ=255)                                                             gtTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCAATACTTTGCCGAAAGGCTATGACGTgat                                                                                   >  2:89236/1‑147 (MQ=255)                                                                  ccACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAAc                                                                              >  8:22611/1‑149 (MQ=255)                                                                   cacaGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAAc                                                                              <  7:281467/148‑1 (MQ=255)                                                                   cacaGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACt                                                                             >  6:90911/1‑149 (MQ=255)                                                                    acaGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTg                                                                            >  8:288404/1‑149 (MQ=255)                                                                    acaGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTg                                                                            >  6:263506/1‑149 (MQ=255)                                                                      aGCGAGCTGGACGTCGAAGTCGAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATACACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCGGAAAGGCTATGACGTGCTGGAACTggg                                                                          >  3:32217/1‑149 (MQ=255)                                                                      aGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTggg                                                                          <  6:111137/149‑1 (MQ=255)                                                                          aGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGaa                                                                      <  5:263506/149‑1 (MQ=255)                                                                                aCGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCt                                                                >  5:358510/1‑149 (MQ=255)                                                                                  gTCGAAGTCAAGATGGTTTGCCAGCGTTGTGTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCttt                                                              >  1:30133/1‑149 (MQ=255)                                                                                     gAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGAc                                                                                        >  5:289135/1‑120 (MQ=255)                                                                                      aaGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCTAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGa                                                           >  5:82540/1‑148 (MQ=255)                                                                                         tCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTg                                                       >  4:359633/1‑149 (MQ=255)                                                                                           aaGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTg                                                       >  8:244027/1‑147 (MQ=255)                                                                                           aaGATGGTTTGCCAGCGTTGTCTAGAGCTGGTCACCCTACCGATCCACAGCGAATGTAGATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTg                                                       >  7:264437/1‑147 (MQ=255)                                                                                              aTGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGc                                                  <  6:358510/149‑1 (MQ=255)                                                                                                 gTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCAt                                                <  7:288404/148‑1 (MQ=255)                                                                                                 gTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCAt                                                <  7:191391/148‑1 (MQ=255)                                                                                                       cAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTc                                           >  6:286861/1‑147 (MQ=255)                                                                                                        aGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCgag                                        <  5:165778/149‑1 (MQ=255)                                                                                                              gTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCgaggagga                                   <  1:344205/148‑1 (MQ=255)                                                                                                               tCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGc                                 >  8:298794/1‑149 (MQ=255)                                                                                                                 ttGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCtt                               >  3:176241/1‑149 (MQ=255)                                                                                                                   gAGCTGGTCACCCTACCGATCCACATCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCttt                              <  8:248457/148‑1 (MQ=255)                                                                                                                     gCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCAt                                                <  4:230951/128‑1 (MQ=255)                                                                                                                     gCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCAt                                                >  3:230951/1‑128 (MQ=255)                                                                                                                     gCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTg                            <  7:298794/148‑1 (MQ=255)                                                                                                                     gCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGc                           <  7:244027/149‑1 (MQ=255)                                                                                                                         gTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGc                       >  8:297862/1‑149 (MQ=255)                                                                                                                            aCCCTACCGATCCACAGCGAATGTACATACGCCGTGTTGATGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGATGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTc                         <  6:319922/144‑1 (MQ=255)                                                                                                                            aCCCTACCGATCCACAGCGAATGTACATACGCCGTGGAGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATAGGACGAGGAGGAGCTTTTGCTc                         >  5:319922/1‑144 (MQ=255)                                                                                                                              ccTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGc                  <  3:359633/149‑1 (MQ=255)                                                                                                                               cTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGc                  >  7:271399/1‑148 (MQ=255)                                                                                                                                  cttcaccACGGCGTATGTACATTCGCTGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGCCAAt              >  5:165021/6‑149 (MQ=255)                                                                                                                                  ccGATCCACAGCGAATGTTCATACGCCGTGGTTAAGGAGGGTGGGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGATCTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGTAGTAGCTTTTGCTCGCCTTGCCAAt              <  6:82540/149‑1 (MQ=255)                                                                                                                                  ccGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGCCAAt              <  8:271399/149‑1 (MQ=255)                                                                                                                                  ccGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGCCAAt              <  7:22611/149‑1 (MQ=255)                                                                                                                                  ccGATCCACAGCGAATGTACATACGCCGTCGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGCCAAt              >  2:71811/1‑149 (MQ=255)                                                                                                                                   ttcaccACGGCGTATGTACATTCGCTGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCga                                         <  8:190776/117‑1 (MQ=39)                                                                                                                                   ttcaccACGGCGTATGTACATTCGCTGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCga                                         >  7:190776/5‑121 (MQ=39)                                                                                                                                   cGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGCCAATc             <  1:194770/149‑1 (MQ=255)                                                                                                                                     aTCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAg                                                    <  2:57238/108‑1 (MQ=255)                                                                                                                                     aTCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAg                                                    >  1:57238/1‑108 (MQ=255)                                                                                                                                       ccACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGCCAATCGTGc         >  7:322692/1‑149 (MQ=255)                                                                                                                                         acaGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGCCAATCGTGCCt       >  7:106534/1‑149 (MQ=255)                                                                                                                                              gAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGCCAATCGTGCCTGCtca  <  8:264437/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    GACAATGTCGGTACGGTGCAGGCGAAGTTCGATTTTGAACGAGATGAACAGCACGTGGTGGTTATCCACAGCGAGCTGGACGTCGAAGTCAAGATGGTTTGCCAGCGTTGTCTTGAGCTGGTCACCCTACCGATCCACAGCGAATGTACATACGCCGTGGTGAAGGAGGGTGCGAATACCCAGTCGTTGCCGAAAGGCTATGACGTGCTGGAACTGGGCGAAGATCCTTTGGATCTGCAGGCATTGGTCGAGGAGGAGCTTTTGCTCGCCTTGCCAATCGTGCCTGCTCA  >  NC_002947/2154949‑2155238

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑