Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,575,190	A→C	7.1%	intergenic (‑176/+34)	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,190	0	A	C	7.1%	196.5 / 6.8	71	intergenic (‑176/+34)	pcaP/galP‑II	porin/porin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base C (3/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

ACCTGTACCTGATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAG  >  NC_002947/1575044‑1575335                                                                                                                                                   |                                                                                                                                                   aCCTGTACCTGATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑tt                                                                                                                                                 >  8:29810/1‑149 (MQ=255)  ccTGTACCTGATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTc                                                                                                                                                >  4:334881/1‑149 (MQ=255)    tGTACCTGATAGTGATCAGGATAGCCTGCAGTCACAACCATCATCAGCTCAGCGAAAGGAGGCAGCTCGCAATGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCGTTCAAAAGGGAGGGCGTCTGGCCTA‑TTagc                                                                                                                                              >  1:129519/1‑149 (MQ=255)                  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                                                                                                                      >  6:340974/1‑137 (MQ=38)                                      ccaccaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTC‑TTCGCGGATAAACCCGCGAATAGGCCaga                                                                                                                      <  5:340974/139‑3 (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                                                                                                   >  5:219273/1‑147 (MQ=255)                                                 agctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGCGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGTGGCCAATTGAGTCACAGAAGATTCAg                                                                                                   <  6:219273/147‑1 (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                                                                                                  <  1:203559/123‑1 (MQ=255)                                                                          aCTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGc                                                                                                  >  2:203559/1‑123 (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                                                   <  8:193376/143‑1 (MQ=255)                                                                                                     tcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGaa                                                   >  7:193376/1‑143 (MQ=255)                                                                                                     tcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTcgc                                              >  2:180206/1‑148 (MQ=255)                                                                                                         gCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACg                                          >  4:124934/1‑148 (MQ=255)                                                                                                               gCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTg                                                           <  6:104165/125‑1 (MQ=255)                                                                                                               gCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTg                                                           >  5:104165/1‑125 (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                 <  8:264685/134‑1 (MQ=255)                                                                                                                                                ccTC‑TTCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGc                 >  7:264685/1‑134 (MQ=255)                                                                                                                                                ccTC‑TTCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAg  >  2:181859/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   ACCTGTACCTGATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTA‑TTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAG  >  NC_002947/1575044‑1575335

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