Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	233,834	A→T	9.8%	intergenic (+22/+41)	PP_0180 → / ← PP_0181	cytochrome c family protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	233,834	0	A	T	9.8%	147.1 / 6.0	59	intergenic (+22/+41)	PP_0180/PP_0181	cytochrome c family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/29);  new base T (3/3);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

GTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTC  >  NC_002947/233686‑233977                                                                                                                                                     |                                                                                                                                                 gTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCcaca                                                                                                                                                  <  5:376540/149‑1 (MQ=255)     gaatggctggGGATTCATGCAGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTTGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAATCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAgtg                                                                                                                                               <  2:249974/145‑3 (MQ=255)       ctggctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGAcc                                                                                                                                            <  4:256717/149‑1 (MQ=255)              gggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCaa                                                                                                                                      <  4:342971/148‑1 (MQ=255)                 aTTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtc                                                                                                                                  >  4:175850/1‑149 (MQ=255)                  ttCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtct                                                                                                                                 <  1:186813/149‑1 (MQ=255)                    cATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtctc                                                                                                                                <  7:151033/148‑1 (MQ=255)                     aTGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAgg                                                                                                                                                <  4:116/131‑1 (MQ=255)                     aTGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAgg                                                                                                                                                >  3:116/1‑131 (MQ=255)                        ccGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATc                                                                                                                           <  1:73899/149‑1 (MQ=255)                        ccGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATc                                                                                                                           <  6:292482/149‑1 (MQ=255)                         cGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCt                                                                                                                          >  7:249305/1‑149 (MQ=255)                         cGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGATCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGCGCTCAGATCt                                                                                                                          >  6:281069/1‑149 (MQ=255)                                tACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCATGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCCCCCCAGGTCTCAGATCTTGCGCaa                                                                                                                   >  2:108539/1‑149 (MQ=255)                                  cTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGc                                                                                                                 <  4:242555/149‑1 (MQ=255)                                  cTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGc                                                                                                                 <  3:175850/149‑1 (MQ=255)                                   tCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCCCACGGGCCGCGCAAGGGGTCAGATCATTCGCaaac                                                                                                                 >  2:346162/1‑146 (MQ=255)                                     gCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCt                                                                                                               <  8:249305/148‑1 (MQ=255)                                      cTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGa                                                                                                             <  2:264117/149‑1 (MQ=255)                                       tGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGAt                                                                                                            <  8:142964/149‑1 (MQ=255)                                                   cGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGAGCATCAggc‑ttgc                                                                                               <  4:224202/149‑1 (MQ=255)                                                      tGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCCGTACCACAGGGACCCCGCAAGGGCTCAGATCTTGCGCAAGCCTGAAGCTCAGGC‑TGGCggt                                                                                            >  4:312133/1‑149 (MQ=255)                                                          cTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtctc                                                                                                                                <  3:205765/110‑1 (MQ=255)                                                                ggCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCtt                                                                                    >  1:232093/1‑147 (MQ=255)                                                                 gCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTga                                                                                  <  3:312133/148‑1 (MQ=255)                                                                  ccATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggc‑tgg                                                                                                >  1:230254/1‑133 (MQ=255)                                                                  ccATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggc‑tgg                                                                                                <  2:230254/133‑1 (MQ=255)                                                                     tCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATg                                                                             >  1:241176/1‑149 (MQ=255)                                                                       gtgtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGtt                                                                           <  1:362220/149‑1 (MQ=255)                                                                         gtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCt                                                                         <  7:224028/149‑1 (MQ=255)                                                                         gtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAACGGCTCAGAAATTGCGCAAGCCAGATGCTAAGGC‑TGGGTGTTGCGCCTTGAGAATGTTca                                                                         >  2:313513/1‑148 (MQ=255)                                                                          ttCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTg                                                                        >  6:308214/1‑149 (MQ=255)                                                                             cTTTACGGCCGCTCGCGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTa                                                                      <  5:73048/148‑1 (MQ=255)                                                                                  cGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAAtt                                                                >  2:58537/1‑149 (MQ=255)                                                                                       gCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACCGGGACCGCGCAAGGTGTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGaca                                                           >  2:60191/1‑149 (MQ=255)                                                                                           gTGGCTGGCCGAAAAGCGCGGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCCCAGATCTTGCGCAAGCCGGATGCTAAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACtt                                                        >  6:109574/1‑148 (MQ=255)                                                                                            tggctggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTc                                                      <  1:346162/149‑1 (MQ=255)                                                                                              gctggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCa                                                     <  1:329158/148‑1 (MQ=255)                                                                                              gctggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCa                                                     <  1:333434/148‑1 (MQ=255)                                                                                                tggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGGCCGCGCCAGGGGTCAGCTCTTTCGCAAGCCGGATGCTCAGGC‑TGGCGGGGGGGGCTTGAGAATGTTGTGTAGCCATTTGACACGTTGAACg                                                  >  8:18391/1‑149 (MQ=255)                                                                                                 ggcCGAAAAGCGCAGTGCTGCGGCCAACAGACTGGAGCCTGCGCCGTCCCCCAGGGGCCGCGCAAGGTGTGAGATCTTTCGCAAGCCGGATGCTCAGGC‑TGGGGGTGGCGCCGTGAGAATGTTCTGTAGCAATTGGACACTTTCAACgg                                                 >  1:208285/1‑149 (MQ=255)                                                                                                  tccGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGc                                                <  6:330117/148‑1 (MQ=255)                                                                                                       aaaGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCAt                                           <  5:308214/149‑1 (MQ=255)                                                                                                            gcAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGAccc                                      <  2:241176/149‑1 (MQ=255)                                                                                                             cAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGAcc                                       >  5:181515/1‑147 (MQ=255)                                                                                                              aGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACTGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCa                                     <  6:181515/148‑1 (MQ=255)                                                                                                                     gggCCAACTGACTGGAGGCTGCGCGGTCCCACAGTGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTc                             <  1:167738/148‑1 (MQ=255)                                                                                                                                ctgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAgccgc                   >  3:46181/1‑148 (MQ=255)                                                                                                                                  agaccttGCGCGGTCCCTGTGGGACCGCGCAGGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCAt                >  6:169036/7‑149 (MQ=255)                                                                                                                                     gCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCg             >  3:295383/1‑149 (MQ=255)                                                                                                                                      cttGCGCGGTCCCTGTGGGACCGCGCAGGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCgg            >  7:126087/3‑149 (MQ=255)                                                                                                                                      ccTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCgg            >  8:54469/1‑149 (MQ=255)                                                                                                                                       ttGCGCGGTCCCTGTGGGACCGCGCAGGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACtt                                                        >  2:151607/2‑104 (MQ=39)                                                                                                                                       ttGCGCGGTCCCTGTGGGACCGCGCAGGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACtt                                                        <  1:151607/103‑1 (MQ=39)                                                                                                                                         ccgttctCCCACAGGGACCGCTCAAGGTCTCAGCTCTTGCGCAAGCCTGCTGCTCAGGC‑TGGCTGTGGCGCCTTGAGTATGTTCTGTTGCATTTGGACACTTTCACCGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCt          <  4:225047/142‑1 (MQ=255)                                                                                                                                           gcgGTCCCTGTGGGACCGCGCAGGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGAt       <  5:169036/149‑1 (MQ=255)                                                                                                                                                 ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCa                 >  6:115810/1‑133 (MQ=255)                                                                                                                                                 ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCa                 <  5:115810/133‑1 (MQ=255)                                                                                                                                                 ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGACCgat   >  6:228848/1‑145 (MQ=255)                                                                                                                                                 ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCACGCTGGGCAGGTACCCGCCATCCGGCTGagaacgg  >  3:225047/1‑142 (MQ=255)                                                                                                                                                  cacaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCGTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGt   >  1:276884/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGC‑TGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTC  >  NC_002947/233686‑233977

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