Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	344,032	A→G	10.8%	intergenic (+121/+68)	PP_5440 → / ← PP_t72	hypothetical protein/tRNA‑Phe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	344,032	0	A	G	10.8%	166.7 / 10.0	65	intergenic (+121/+68)	PP_5440/PP_t72	hypothetical protein/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (25/33);  new base G (4/3);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

CCATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGAC  >  NC_002947/343895‑344178                                                                                                                                          |                                                                                                                                                     ccATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑c                                                                                                                                           <  5:45277/148‑1 (MQ=255)  cATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑c                                                                                                                                           >  6:111438/1‑147 (MQ=255)      tCGGACGTGTTAGGTTTGGCCCGCAGTATGCTCATGGCTGAGGGTTTCTTGATTGCATGGCTATCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGTCGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑Gc                                                                                                                                     <  3:252192/148‑1 (MQ=255)       cGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACaaaa                                                                                                                                                     <  5:53609/133‑1 (MQ=255)       cGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACaaaa                                                                                                                                                     >  6:53609/1‑133 (MQ=255)        ggACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCggg                                                                                                                                  >  7:234283/1‑149 (MQ=255)        ggACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCgg                                                                                                                                   <  3:191861/148‑1 (MQ=255)            gAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTc                                                                                                                              >  4:372050/1‑149 (MQ=255)            gAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTc                                                                                                                              <  7:89830/149‑1 (MQ=255)              gTAAGGTTTGGCCCGCAGTCCGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCtt                                                                                                                            <  7:145220/149‑1 (MQ=255)               tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCtt                                                                                                                            <  8:319139/148‑1 (MQ=255)                aaGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCttgt                                                                                                                          >  4:266707/1‑149 (MQ=255)                  ggTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑‑CTAGTGCGGGTTTTttgtt                                                                                                                         <  3:111325/148‑1 (MQ=39)                  ggTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑‑CTAGTGCGGGTTTTttgtt                                                                                                                         >  4:111325/1‑148 (MQ=39)                    tttGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTc                                                                                                                      <  8:234283/149‑1 (MQ=255)                         cccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTc                                                                                                                 <  4:139715/149‑1 (MQ=255)                           cGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑‑CTAGTGCGGGTTTTttgtt                                                                                                                         >  3:104485/1‑139 (MQ=39)                           cGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑‑CTAGTGCGGGTTTTttgtt                                                                                                                         <  4:104485/139‑1 (MQ=39)                           cGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCCGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCgg                                                                                                               >  3:216637/1‑149 (MQ=255)                                    cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTg                                                                                                      >  5:185741/1‑149 (MQ=255)                                    cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTAGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTg                                                                                                      <  7:239151/149‑1 (MQ=255)                                       aTGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCaa                                                                                                   >  5:325033/1‑149 (MQ=255)                                       aTGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCCGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGc                                                                                                     <  1:368381/147‑1 (MQ=255)                                                   ttCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCt                                                                                       >  3:269149/1‑149 (MQ=255)                                                    tctcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCt                                                                                       <  3:266707/148‑1 (MQ=255)                                                      tcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑‑CTAGTGCGGGTTTTttg                                                                                                                           <  4:153474/110‑1 (MQ=38)                                                      tcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑‑CTAGTGCGGGTTTTttg                                                                                                                           >  3:153474/1‑110 (MQ=38)                                                        gATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTg                                                                                     <  1:334482/146‑1 (MQ=255)                                                            gCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAggtg                                                                               <  7:114582/148‑1 (MQ=255)                                                                 gCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGcc                                                                          <  7:84517/148‑1 (MQ=255)                                                                  cTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCa                                                                        <  4:269149/149‑1 (MQ=255)                                                                    aGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCa                                                                        >  1:225902/1‑147 (MQ=255)                                                                       cTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGAc                                                                   >  4:20630/1‑149 (MQ=255)                                                                       cTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCCGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGAc                                                                   <  4:216637/149‑1 (MQ=255)                                                                           aGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGa                                                                <  5:182384/148‑1 (MQ=255)                                                                            ggTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGt                                                              >  4:238891/1‑149 (MQ=255)                                                                             gTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTc                                                             <  3:372050/149‑1 (MQ=255)                                                                              ttGCGTTGTGATGTCTCGTGTTCTTCGACTCCGACTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCg                                                            <  6:325033/149‑1 (MQ=255)                                                                               tGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCg                                                            >  8:102623/1‑148 (MQ=255)                                                                                gCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGa                                                           <  1:148893/148‑1 (MQ=255)                                                                                  gTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAcc                                                        >  8:185542/1‑149 (MQ=255)                                                                                      tgAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTTG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAccgcc                                                     >  4:2402/1‑148 (MQ=255)                                                                                        aGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGa                                                   <  6:77455/148‑1 (MQ=255)                                                                                          gTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGACCAGAGACGGAGTCGAACCGCCGacac                                                >  6:227440/1‑149 (MQ=255)                                                                                             tcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCAGCTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGGGTCGAACCGCCGacac                                                >  1:293081/1‑147 (MQ=255)                                                                                             tcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAg                                             >  8:330985/1‑149 (MQ=255)                                                                                             tcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAg                                             >  4:179702/1‑149 (MQ=255)                                                                                               gtgtTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGCCACGAGGa                                           >  8:156773/1‑149 (MQ=255)                                                                                               gtgtTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGa                                              <  6:375460/146‑1 (MQ=255)                                                                                                      tcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCa                                     <  2:185123/148‑1 (MQ=255)                                                                                                        gACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTTG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATc                                  <  3:2402/149‑1 (MQ=255)                                                                                                        gACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAAt                                   <  3:238891/148‑1 (MQ=255)                                                                                                        gACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATc                                  >  8:109065/1‑149 (MQ=255)                                                                                                        gACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATc                                  >  3:270372/1‑149 (MQ=255)                                                                                                          cgccgcCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCct                                <  7:102623/149‑1 (MQ=255)                                                                                                           gccgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCctc                               <  5:111438/149‑1 (MQ=255)                                                                                                                cTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGctc                          >  1:69336/1‑149 (MQ=255)                                                                                                                    aGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAcc                      <  6:185741/149‑1 (MQ=255)                                                                                                                      ttAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGa                    >  7:3127/1‑149 (MQ=255)                                                                                                                          aaaTCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                <  2:36844/149‑1 (MQ=255)                                                                                                                             tCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCt             <  3:179702/149‑1 (MQ=255)                                                                                                                                   acaacaAGAAACCCGC‑‑CTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTgg       >  5:378818/1‑149 (MQ=39)                                                                                                                                   acaacaAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTgg       >  5:367500/1‑149 (MQ=255)                                                                                                                                    caacaaAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTgg       <  7:330985/148‑1 (MQ=255)                                                                                                                                     aacaaAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                >  5:263257/1‑138 (MQ=255)                                                                                                                                     aacaaAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                <  6:263257/138‑1 (MQ=255)                                                                                                                                     aacaaAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc     <  3:20630/149‑1 (MQ=255)                                                                                                                                        aaaaaaCCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac  >  7:35371/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                                     CCATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGAC  >  NC_002947/343895‑344178

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