Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	344,041	Δ1 bp	11.6%	intergenic (+130/+59)	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,041	0	A	.	11.6%	204.4 / 22.3	69	intergenic (+130/+59)	PP_5440/PP_t72	hypothetical protein/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base . (5/3);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CCATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATT  >  NC_002947/343895‑344188                                                                                                                                                   |                                                                                                                                                      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)        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                                                                                                                                     >  3:153474/1‑110 (MQ=38)                                                      tcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑‑CTAGTGCGGGTTTTttg                                                                                                                                     <  4:153474/110‑1 (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)                                                                                                                                        aagaaaCCCGC‑‑CTAGTGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac            >  2:375028/4‑149 (MQ=255)                                                                                                                                        aaaaaaCCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac            >  7:35371/1‑149 (MQ=255)                                                                                                                                               ccGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTg      >  7:146772/1‑148 (MQ=255)                                                                                                                                                  cA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAAtt  <  8:146772/149‑1 (MQ=255)                                                                                                                                                   a‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAAtt  >  6:315927/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CCATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCA‑CTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATT  >  NC_002947/343895‑344188

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