Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,099,701	C→T	6.1%	G14G (GGC→GGT)	PP_4488 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,099,701	0	C	T	6.1%	189.9 / 2.3	65	G14G (GGC→GGT)	PP_4488	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/33);  new base T (2/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

CGTTAATCTCCATCTAATCCCGATAGTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTAC  >  NC_002947/5099562‑5099845                                                                                                                                            |                                                                                                                                                 cGTTAATCTCCATCTAATCCCGATAGTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATg                                                                                                                                         <  1:30411/149‑1 (MQ=255)         tcCATCTAATCCCGATAGTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCg                                                                                                                                  >  7:145651/1‑148 (MQ=255)           cATCTAATCCCGATAGTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGaa                                                                                                                                <  3:274469/148‑1 (MQ=255)                 aTCCCGATAGTTAAAACTAGCAATACTCGCCGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAg                                                                                                                         <  7:221253/149‑1 (MQ=255)                     cGATAGTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGc                                                                                                                     <  8:145651/149‑1 (MQ=255)                          gTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAg                                                                                                                 >  1:137121/1‑148 (MQ=255)                          gTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGc                                                                                                                <  7:95668/149‑1 (MQ=255)                           ttAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAg                                                                                                                 >  1:250830/1‑147 (MQ=255)                           ttAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCa                                                                                                               >  8:227915/1‑149 (MQ=255)                              aaaCTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGa                                                                                                             >  1:273307/1‑148 (MQ=255)                              aaaCTAGCAATACGCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGCTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGa                                                                                                             <  6:170263/148‑1 (MQ=255)                              aaaCTAGAATTACTCGCGGCGCAATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGa                                                                                                             <  2:273307/148‑1 (MQ=255)                                aCTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTc                                                                                                          >  5:321932/1‑149 (MQ=255)                                 cTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTcc                                                                                                         >  5:99359/1‑149 (MQ=255)                                    gCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCgg                                                                                                       <  5:153326/148‑1 (MQ=255)                                    gCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGt                                                                                                      >  4:274959/1‑149 (MQ=255)                                        tACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGaa                                                                                                  <  3:274959/149‑1 (MQ=255)                                        tACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGaa                                                                                                  <  7:227915/149‑1 (MQ=255)                                          cTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAac                                                                                                 <  6:99359/148‑1 (MQ=255)                                               ggcgCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAacagtaca                                                                                           <  7:193181/149‑1 (MQ=255)                                                gcgcCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAacagtacag                                                                                          <  7:98785/149‑1 (MQ=255)                                                  gcCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAacagtacagt                                                                                         >  7:52980/1‑148 (MQ=255)                                                   ccATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTa                                                                                       >  1:27494/1‑149 (MQ=255)                                                    cATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTa                                                                                       >  7:172778/1‑148 (MQ=255)                                                       aaGCGCCGCGTGCCCCCTCACGACCTTAGACAGGGATATAGAAAATGAAAGCTTTACTGGGATTGGTACTTGGCTGTGTTTGCGGCGCGGCATTGGCCGGCGAAGGCAAAGATGCGGCGCAGCTTCCGGTTGCACAGTACAGTTACtgg                                                                                   >  6:235756/1‑147 (MQ=255)                                                             cgcgTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCAc                                                                             >  5:294200/1‑149 (MQ=255)                                                                   cccccTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGa                                                                        >  8:191136/1‑148 (MQ=255)                                                                      ccTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCTGTGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAacagtaca                                                                                           <  6:224301/126‑1 (MQ=255)                                                                      ccTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGTGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAacagtaca                                                                                           <  2:63676/126‑1 (MQ=255)                                                                      ccTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGTGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAacagtaca                                                                                           >  1:63676/1‑126 (MQ=255)                                                                      ccTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGTGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAacagtaca                                                                                           >  5:224301/1‑126 (MQ=255)                                                                       cTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTg                                                                   >  3:222760/1‑149 (MQ=255)                                                                       cTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTg                                                                   <  8:172778/149‑1 (MQ=255)                                                                          accaCCTTAGTCAGGTAAATAGAATATGAAAGCTTTACTGGTATTGGTACTTGGGAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGc                                                                  <  6:321932/144‑1 (MQ=255)                                                                          aCGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGccc                                                                >  4:207420/1‑149 (MQ=255)                                                                          aCGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGACGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAaca                                                                                                <  7:252232/117‑1 (MQ=255)                                                                          aCGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGACGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAaca                                                                                                >  8:252232/1‑117 (MQ=255)                                                                               cTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGtt                                                           <  7:151065/149‑1 (MQ=255)                                                                                 tAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTAt                                                         >  6:171073/1‑149 (MQ=255)                                                                                   gACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTAttt                                                       <  6:294200/149‑1 (MQ=255)                                                                                      aGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCa                                                    <  3:104962/149‑1 (MQ=255)                                                                                        gTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAacagtacag                                                                                          >  8:6777/1‑109 (MQ=255)                                                                                           aaTAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGc                                               <  8:52980/149‑1 (MQ=255)                                                                                               gAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGaa                                            <  2:137121/148‑1 (MQ=255)                                                                                               gAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAg                                           >  6:311189/1‑149 (MQ=255)                                                                                               gAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGACCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAg                                           >  3:290458/1‑149 (MQ=255)                                                                                                   aTGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGcc                                       <  4:290458/149‑1 (MQ=255)                                                                                                     gAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTa                                                          >  6:278828/1‑128 (MQ=255)                                                                                                     gAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTa                                                          <  5:278828/128‑1 (MQ=255)                                                                                                            ttACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAAtgtgtg                              >  5:102018/1‑149 (MQ=255)                                                                                                             tACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAAtgtgt                               <  3:207420/147‑1 (MQ=255)                                                                                                              aCTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGc                             <  5:171073/148‑1 (MQ=255)                                                                                                              aCTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGc                             >  6:7592/1‑148 (MQ=255)                                                                                                              aCTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGc                             <  5:311189/148‑1 (MQ=255)                                                                                                               cTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGGCCAATGTGTGCg                            >  4:207763/1‑148 (MQ=255)                                                                                                                tGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCAGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGATATTTCCATGAGCGAAGTGACCAATCTGTGCg                            >  4:184504/1‑147 (MQ=255)                                                                                                                  gTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGaa                          <  2:27494/147‑1 (MQ=255)                                                                                                                     ttGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGt                     <  5:220317/149‑1 (MQ=255)                                                                                                                        gTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTg                                                                          >  5:271014/1‑93 (MQ=255)                                                                                                                        gTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTg                                                                          <  6:271014/93‑1 (MQ=255)                                                                                                                         tACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCg                 <  8:12633/149‑1 (MQ=255)                                                                                                                         tACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCg                 <  8:17311/149‑1 (MQ=255)                                                                                                                               gCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCg            >  7:50753/1‑148 (MQ=255)                                                                                                                                       tGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGAc      <  2:250830/146‑1 (MQ=255)                                                                                                                                       tGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGAc      <  5:7592/146‑1 (MQ=255)                                                                                                                                          ggcgcggcAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTAc  >  8:119546/1‑147 (MQ=255)                                                                                                                                            |                                                                                                                                                 CGTTAATCTCCATCTAATCCCGATAGTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGGTAAATAGAAAATGAAAGCTTTACTGGTATTGGTACTTGGCAGTCTTTGCGGCGCGGCAATGGCCGGCGAAGCCAAAGATGCCGAGCAGATTCCGGTTGAACAGTACAGTTACTCGCAGCACCTGGACATTGCCCGCGTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTAC  >  NC_002947/5099562‑5099845

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