Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,575,198	G→A	7.8%	intergenic (‑184/+26)	pcaP ← / ← galP‑II	porin/porin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,575,198	0	G	A	7.8%	223.8 / 6.5	77	intergenic (‑184/+26)	pcaP/galP‑II	porin/porin‑like protein
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base A (4/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

ATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTG  >  NC_002947/1575061‑1575346                                                                                                                                           |                                                                                                                                                     aTCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCcgcg                                                                                                                                            >  4:223022/1‑148 (MQ=255) aTCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCcgcg                                                                                                                                            >  2:137407/1‑148 (MQ=255) aTCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTAGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCcgcg                                                                                                                                            >  1:153002/1‑148 (MQ=255)  tCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGaa                                                                                                                                          <  6:332256/149‑1 (MQ=255)  tCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGaa                                                                                                                                          <  4:275346/149‑1 (MQ=255)     ggATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAagag                                                                                                                                       <  3:223022/149‑1 (MQ=255)           ccTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAAt                                                                                                                                 >  4:310927/1‑149 (MQ=255)            cTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAAtt                                                                                                                                >  8:321198/1‑149 (MQ=255)               cAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                              <  5:285412/148‑1 (MQ=255)                    accaccATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCAGCGAAGAGGCCAATTGAGTcaca                                                                                                                        >  3:73335/1‑149 (MQ=255)                     ccaccaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTC‑TTCGCGGATAAACCCGCGAATAGGCCaga                                                                                                                                 <  5:340974/139‑3 (MQ=38)                     ccaccaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTC‑TTCGCGGATAAACCCGCGAATAGGCCaga                                                                                                                                 >  6:340974/1‑137 (MQ=38)                     ccaccaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTca                                                                                                                          >  5:365680/1‑146 (MQ=255)                     ccaccaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACag                                                                                                                       <  3:334881/149‑1 (MQ=255)                       accaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACagaa                                                                                                                     <  6:35197/149‑1 (MQ=255)                                agctcagcGAATGGCGGCAGCTCTCACTGGGGAATTCATGGCAGGGTTTTCATCACTCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCTTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGATGCCAATTGAGTCACAGAAGATTCAGc                                                                                                             <  2:129519/148‑1 (MQ=255)                                agctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGCGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGTGGCCAATTGAGTCACAGAAGATTCAg                                                                                                              <  6:219273/147‑1 (MQ=255)                                agctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGCGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGTGGCCAATTGAGTCACAGAAGATTCAg                                                                                                              >  5:219273/1‑147 (MQ=255)                                agctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCg                                                                                                            <  1:35897/149‑1 (MQ=255)                                 gctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAg                                                                                                                             >  7:173230/1‑131 (MQ=255)                                 gctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAg                                                                                                                             <  8:173230/131‑1 (MQ=255)                                     agcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATaa                                                                                                       >  7:162430/1‑149 (MQ=255)                                     agcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATaa                                                                                                       <  7:321198/149‑1 (MQ=255)                                        gAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACt                                                                                                     >  4:98281/1‑148 (MQ=255)                                         aaTGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTca                                                                                                   >  5:96973/1‑149 (MQ=255)                                          aTGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTcac                                                                                                  <  2:150890/149‑1 (MQ=255)                                                gcAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATc                                                                                              >  7:16804/1‑147 (MQ=255)                                                    cTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCgg                                                                                        <  8:162430/149‑1 (MQ=255)                                                       gCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTtc                                                                                     >  5:39232/1‑149 (MQ=255)                                                         aCTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGc                                                                                                             >  2:203559/1‑123 (MQ=255)                                                         aCTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGc                                                                                                             <  1:203559/123‑1 (MQ=255)                                                         aCTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTtct                                                                                    >  5:96365/1‑148 (MQ=255)                                                           tGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCAt                                                                                 >  4:72757/1‑149 (MQ=255)                                                             tgCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATc                                                                                <  4:110815/147‑1 (MQ=255)                                                             ggCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCAt                                                                                 <  2:153002/147‑1 (MQ=255)                                                                 aTTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAAcg                                                                           >  2:95439/1‑149 (MQ=255)                                                                        gCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGt                                                                    <  3:310927/149‑1 (MQ=255)                                                                          aGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAg                                                                                                              >  7:249956/1‑105 (MQ=255)                                                                          aGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAg                                                                                                              <  8:249956/105‑1 (MQ=255)                                                                           gggTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCt                                                                 <  6:96973/149‑1 (MQ=255)                                                                              ttttCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGaa                                                              <  6:365680/149‑1 (MQ=255)                                                                                ttCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGt                                                            >  7:44545/1‑149 (MQ=255)                                                                                  catcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTc                                                           <  2:30840/148‑1 (MQ=255)                                                                                    tcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGaa                                                              >  7:193376/1‑143 (MQ=255)                                                                                    tcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGaa                                                              <  8:193376/143‑1 (MQ=255)                                                                                    tcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTcgc                                                         >  2:180206/1‑148 (MQ=255)                                                                                        gCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACg                                                     >  4:124934/1‑148 (MQ=255)                                                                                              gCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTg                                                                      >  5:104165/1‑125 (MQ=255)                                                                                              gCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTg                                                                      <  6:104165/125‑1 (MQ=255)                                                                                              gCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGt                                               <  3:98281/148‑1 (MQ=255)                                                                                               ccAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTc                                              <  6:96365/148‑1 (MQ=255)                                                                                                cAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTc                                              >  1:241788/1‑148 (MQ=255)                                                                                                cAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTc                                              >  5:191795/1‑147 (MQ=255)                                                                                                     atgCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgc                                       <  8:16804/149‑1 (MQ=255)                                                                                                        ccTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGtt                                         <  2:3542/144‑1 (MQ=255)                                                                                                        ccTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGtt                                         >  1:3542/1‑144 (MQ=255)                                                                                                         cTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACagaa                                                                                                                     >  6:88782/1‑67 (MQ=255)                                                                                                         cTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAg                                                                                                              >  5:238884/1‑74 (MQ=255)                                                                                                         cTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAg                                                                                                              <  6:238884/74‑1 (MQ=255)                                                                                                         cTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTcgc                                                         >  7:327470/1‑127 (MQ=255)                                                                                                         cTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTcgc                                                         <  8:327470/127‑1 (MQ=255)                                                                                                             aaaaGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAttt                                >  6:287152/1‑148 (MQ=255)                                                                                                                      ggCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAg                       <  1:180206/148‑1 (MQ=255)                                                                                                                         gTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTcc                   <  1:137407/149‑1 (MQ=255)                                                                                                                              gCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCtt                 <  8:44545/146‑1 (MQ=255)                                                                                                                               ccTC‑TTCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGc                            >  7:264685/1‑134 (MQ=255)                                                                                                                               ccTC‑TTCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGc                            <  8:264685/134‑1 (MQ=255)                                                                                                                               ccTC‑TTCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAg             >  2:181859/1‑149 (MQ=255)                                                                                                                                 ctcttCGCGGATAAACCCGCGATTAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGc            >  3:368256/4‑149 (MQ=255)                                                                                                                                    ttCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTg          <  6:39232/148‑1 (MQ=255)                                                                                                                                    ttCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGc         >  6:272650/1‑149 (MQ=255)                                                                                                                                      cgcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGccc       >  6:293221/1‑149 (MQ=255)                                                                                                                                      cgcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGcc        <  7:29810/148‑1 (MQ=255)                                                                                                                                      cgcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGcc        <  5:287152/148‑1 (MQ=255)                                                                                                                                       gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGccc       <  1:329126/148‑1 (MQ=255)                                                                                                                                       gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGccc       >  6:271141/1‑148 (MQ=255)                                                                                                                                       gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGAGCTGAAGTCGCGACGCATGGTCGa                                            >  8:145857/1‑111 (MQ=255)                                                                                                                                         gggTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGAc    <  3:72757/149‑1 (MQ=255)                                                                                                                                           gTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTg  >  1:186459/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                     ATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTG  >  NC_002947/1575061‑1575346

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