Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,575,203	T→C	9.2%	intergenic (‑189/+21)	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,203	0	T	C	9.2%	174.3 / 7.6	65	intergenic (‑189/+21)	pcaP/galP‑II	porin/porin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base C (4/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

TAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGCACCA  >  NC_002947/1575056‑1575351                                                                                                                                                    |                                                                                                                                                     tAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTAt                                                                                                                                                      <  2:371491/148‑1 (MQ=255)     gATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCcgc                                                                                                                                                  >  2:170934/1‑148 (MQ=255)      aTCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCcgcg                                                                                                                                                 >  1:76808/1‑148 (MQ=255)       tCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCcgcg                                                                                                                                                 >  7:224671/1‑147 (MQ=255)           gATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAATAgg                                                                                                                                           <  5:364260/149‑1 (MQ=255)            aTAGCCTGCAGTAACCACCATCTTCAGCGCATCGAATGGCGGCAGCTCGCACTTTGCAATTCATGGCAGTGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGTGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGc                                                                                                                                          <  4:292205/149‑1 (MQ=255)               gCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCaa                                                                                                                                       <  6:221253/149‑1 (MQ=255)                 cTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAAtt                                                                                                                                     <  2:76808/149‑1 (MQ=255)                  tGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGGCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAAtt                                                                                                                                     >  8:53168/1‑148 (MQ=255)                   cgatctaccaccATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGCGCAATTCATGGCAGGGTTTTCCTCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                                   <  7:123917/143‑1 (MQ=255)                     aGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTGTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGt                                                                                                                                 <  3:70864/149‑1 (MQ=255)                       tCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTc                                                                                                                                >  7:87399/1‑148 (MQ=255)                         accaccATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATTCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGATGACAATTGAGTca                                                                                                                               >  6:219245/1‑147 (MQ=255)                         accaccATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTca                                                                                                                               >  4:164784/1‑147 (MQ=255)                         accaccATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTca                                                                                                                               >  4:367769/1‑147 (MQ=255)                          ccaccaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACag                                                                                                                            <  5:57432/149‑1 (MQ=255)                          ccaacaACATCAGCTCAGCGAAAGGCGCCAGCACCAACTCGGCAAATCATCCCAGGGTTAACATCACGATGATCACAAAGATCCCATAAAAATCGAGCGCCACTGGCCAATTCCCGCGAATATCCGCGAAGAGGCCAATTCAgagaacg                                                                                                                            >  8:123917/1‑143 (MQ=255)                           caccaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACaga                                                                                                                           <  1:95635/149‑1 (MQ=255)                            accaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACagaa                                                                                                                          <  2:393378/149‑1 (MQ=255)                                        tcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGa                                                                                                               >  7:227695/1‑148 (MQ=255)                                          agcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATa                                                                                                             >  2:269882/1‑148 (MQ=255)                                          agcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGATAAACCCGCGAatag                                                                                                                                            >  8:81547/1‑117 (MQ=255)                                          agcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGATAAACCCGCGAatag                                                                                                                                            <  7:81547/117‑1 (MQ=255)                                              aaTGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTca                                                                                                        >  5:330009/1‑149 (MQ=255)                                                   cggcAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATc                                                                                                   <  1:269882/149‑1 (MQ=255)                                                          tCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGt                                                                                            >  4:249593/1‑149 (MQ=255)                                                          tCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGt                                                                                            >  8:330561/1‑149 (MQ=255)                                                            gCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGtt                                                                                           <  8:224671/148‑1 (MQ=255)                                                                tGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCAAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCTGTTCTCAt                                                                                      <  5:219245/149‑1 (MQ=255)                                                                  ggCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGGCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATc                                                                                     <  7:53168/148‑1 (MQ=255)                                                                  ggCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATc                                                                                     >  2:385269/1‑148 (MQ=255)                                                                  ggCAATTCATGGCAGGGTATTCATCACGCTGTTTCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGCTTTATCCGCGAAGAGGCCAATAGAGTCACAGAACATTCAGCGGATAACTCACAATCAAGCGGTACTCATc                                                                                     >  4:321572/1‑148 (MQ=255)                                                                     aaTTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAAc                                                                                 >  3:81993/1‑149 (MQ=255)                                                                         cATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGc                                                                                                                  >  4:156231/1‑112 (MQ=255)                                                                         cATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGc                                                                                                                  <  3:156231/112‑1 (MQ=255)                                                                          aTGGCATGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGt                                                                             <  3:161900/148‑1 (MQ=255)                                                                           tGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTg                                                                           <  6:330009/149‑1 (MQ=255)                                                                                      tcatcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGt                                                                 <  4:43918/148‑1 (MQ=255)                                                                                        atcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTCGTGCTGAAGTcgc                                                              >  2:199867/1‑149 (MQ=255)                                                                                           aCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGa                                                            <  7:330561/148‑1 (MQ=255)                                                                                           aCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGAc                                                           <  7:160022/149‑1 (MQ=255)                                                                                             gCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACg                                                          >  6:308237/1‑148 (MQ=255)                                                                                                  tGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGc                                                                                                                  >  2:302326/1‑87 (MQ=255)                                                                                                  tGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGc                                                                                                                  <  1:302326/87‑1 (MQ=255)                                                                                                   gCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGt                                                    <  5:118345/148‑1 (MQ=255)                                                                                                    ccAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTc                                                   >  4:379852/1‑148 (MQ=255)                                                                                                     cAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCCCAATACAGCGGTCCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATggac                                                   >  8:353490/1‑145 (MQ=255)                                                                                                               ttCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCGGAGAAGAGGCCAAGTGAGTCACAGAAGATTCAGCGGAGAGCTCACAATCAAGCGGGGCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGGTGCGCCAt                                       >  2:364817/1‑149 (MQ=255)                                                                                                                tCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgc                                            <  5:305771/143‑1 (MQ=255)                                                                                                                tCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgc                                            >  6:305771/1‑143 (MQ=255)                                                                                                                         agggCGTCTGGCCTATTCGGGGGGTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCTGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGGCGAGTTGCGCCATTTCa                                   <  5:258111/143‑1 (MQ=255)                                                                                                                         agggCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACGCACAATCAAGCGGTTCTCATCGAACGCGATGGTGCTGAAGTCGCCACGCATGGTCGAGTTGCGCCAttgaa                                   >  6:258111/1‑140 (MQ=255)                                                                                                                             caattGGCCTCTTCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTc                         >  3:325474/5‑149 (MQ=38)                                                                                                                              gTCTGTCCTATTCTCGGTTTTATCCTCGAAGAGGCCAATTGAGTCACAGAAGATTCATCGGATAACTCACAATCAAGCGTTTCTCTTCGAACGCGTTGGTGATGATGTCGCGAAGAATGGTCGAGTTGCGCCATTTCAGCGACAGGTc                         <  4:81993/148‑1 (MQ=255)                                                                                                                              gTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTcc                        <  3:367769/149‑1 (MQ=255)                                                                                                                                 tGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCtt                      >  1:225537/1‑148 (MQ=255)                                                                                                                                  ggCCTCTTCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTg                     >  7:169694/1‑148 (MQ=255)                                                                                                                                   gCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGaa                   <  1:170934/149‑1 (MQ=255)                                                                                                                                   gCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGa                    <  2:225537/148‑1 (MQ=255)                                                                                                                                     ctcttCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGc                 >  3:123030/4‑149 (MQ=255)                                                                                                                                     cTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCt                       <  5:129801/143‑1 (MQ=255)                                                                                                                                     cTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGAGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCt                       >  6:129801/1‑143 (MQ=255)                                                                                                                                        ttCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTg               <  8:227695/148‑1 (MQ=255)                                                                                                                                            cgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGa          <  5:308237/149‑1 (MQ=255)                                                                                                                                             gggTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGAc         >  2:183448/1‑149 (MQ=255)                                                                                                                                              ggTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCGGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACt        <  8:87399/149‑1 (MQ=255)                                                                                                                                                    tCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGcacca  >  3:279646/1‑149 (MQ=255)                                                                                                                                                    tCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGcacca  >  1:56391/1‑149 (MQ=255)                                                                                                                                                    tCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGcacc   <  3:249593/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGCACCA  >  NC_002947/1575056‑1575351

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