Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,425,180	C→A	6.5%	intergenic (‑238/‑562)	PP_5647 ← / → PP_3917	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,425,180	0	C	A	6.5%	194.2 / 3.8	62	intergenic (‑238/‑562)	PP_5647/PP_3917	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base A (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

GCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAAGGGTCGGTTATCAGTGTAATCGTCGCGATCCT  >  NC_002947/4425032‑4425321                                                                                                                                                      |                                                                                                                                              gCCCTAATCCTTGCTTATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGc                                                                                                                                               <  7:200208/149‑1 (MQ=255)       aTCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTACCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCaa                                                                                                                                         <  2:367460/149‑1 (MQ=255)         ccTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGAGTTCAAt                                                                                                                                        >  2:274345/1‑148 (MQ=255)         ccTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGAGTTCAAt                                                                                                                                        >  2:269376/1‑148 (MQ=255)          cTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTa                                                                                                                                      >  1:361688/1‑149 (MQ=255)           ttGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTAc                                                                                                                                     >  6:180414/1‑149 (MQ=255)             gCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACg                                                                                                                                    <  8:343662/148‑1 (MQ=255)             gCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGc                                                                                                                                   >  2:307249/1‑149 (MQ=255)               tGATAGATATCCGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGCCTCCGCGTGCAATTACGCtt                                                                                                                                 >  8:197610/1‑149 (MQ=255)                gatagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTg                                                                                                                                >  3:185857/1‑149 (MQ=255)                gatagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGACCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTg                                                                                                                                >  7:137007/1‑149 (MQ=255)                 atagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTgg                                                                                                                               >  8:330218/1‑149 (MQ=255)                 atagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTgg                                                                                                                               >  2:165297/1‑149 (MQ=255)                          aGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACt                                                                                                                       >  3:327566/1‑148 (MQ=255)                           gCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTc                                                                                                                      >  5:68966/1‑148 (MQ=255)                              ggACAGCGCCAGTAGAGAACCTAGCCCGATCAGTGAGCCCGTCCAGCAGACAGAAGTTCAATTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGGAGCTACTCCAAt                                                                                                                  <  4:176825/149‑1 (MQ=255)                                    cgcCAGTAGAGAAACGAGCCAGATCACTGAGCCCGACCACCAGACAGAAGTTCAAGTAGCAAAAAGCCAAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGAATAGACTCAGCGTAAAATTACTCTTGGCAGCTACTCCACTGGgcgc                                                                                                            >  8:294848/1‑149 (MQ=255)                                     gcCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCCGCAGACAGCAGTTCAAGTAGCAAAAAGCCCAGCGAGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGATCAATTACGCTTGGCAGCTACTCCAATGCGCGCt                                                                                                           >  8:233840/1‑149 (MQ=255)                                      ccAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCtt                                                                                                          <  5:180414/149‑1 (MQ=255)                                       cAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTtg                                                                                                         <  2:361688/149‑1 (MQ=255)                                       cAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTtg                                                                                                         >  3:325194/1‑149 (MQ=255)                                         gTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTATTAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTtgtg                                                                                                       <  1:307249/149‑1 (MQ=255)                                                  cGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCta                                                                                                >  2:28164/1‑147 (MQ=255)                                                        cGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGCTAGCTATAGCGCTc                                                                                        >  4:109633/1‑149 (MQ=255)                                                             aCTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCAt                                                                                   >  7:22477/1‑149 (MQ=255)                                                             aCTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCAt                                                                                   <  4:327566/149‑1 (MQ=255)                                                                 aGCCCGTCCTGCAGACAGTAGTTCAAGTAGCAAAAAGGCCATCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTgg                                                                                <  2:357305/148‑1 (MQ=255)                                                                 aGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTg                                                                                 <  5:103094/147‑1 (MQ=255)                                                                 aGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCACTTGTGATAGCTATAGCGCTCGGCATTg                                                                                 <  1:165297/147‑1 (MQ=255)                                                                    ccGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATa                                                                                                    >  2:320525/1‑125 (MQ=255)                                                                     cGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGcc                                                                           <  8:12929/149‑1 (MQ=255)                                                                      gTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGccc                                                                          <  5:353353/149‑1 (MQ=255)                                                                      gTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGccc                                                                          >  6:353353/1‑149 (MQ=255)                                                                           gcagACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGTAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGAt                                                                     <  7:330218/149‑1 (MQ=255)                                                                           gcagACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGAt                                                                     <  3:169791/149‑1 (MQ=255)                                                                           gcagACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGAt                                                                     >  2:250365/1‑149 (MQ=255)                                                                               acagaAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGt                                                                 >  1:110995/1‑149 (MQ=255)                                                                                      ttCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTg                                                           >  1:269655/1‑149 (MQ=255)                                                                                      ttCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTgg                                                          <  1:28164/149‑1 (MQ=255)                                                                                        cAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTgg                                                          <  8:22477/147‑1 (MQ=255)                                                                                         aaGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTATAATCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGATCTGCATTGGATGCCCTCGATGGGTCCGGTGGGc                                                        <  7:233840/148‑1 (MQ=255)                                                                                          aGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCt                                                      >  1:140285/1‑149 (MQ=255)                                                                                            tAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACGCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTg                                                    >  7:318834/1‑149 (MQ=255)                                                                                              gCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGt                                                   <  7:288179/148‑1 (MQ=255)                                                                                               cAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGAGTTCAATTACGCTTGGCAGCTACTCCAATGGGCTCTTGTGATAGCTTTAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGttt                                                 <  1:269376/149‑1 (MQ=255)                                                                                               cAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGAGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTTTAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGttt                                                 <  1:274345/149‑1 (MQ=255)                                                                                                 aaaaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGa                                                                      <  4:357153/126‑1 (MQ=255)                                                                                                 aaaaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGa                                                                      >  3:357153/1‑126 (MQ=255)                                                                                                   aaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTgc                                             <  2:269655/149‑1 (MQ=255)                                                                                                    aGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTgcg                                            <  6:68966/149‑1 (MQ=255)                                                                                                          gcgcGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTc                                      >  7:14289/1‑149 (MQ=255)                                                                                                                  cggggctTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAAGGGt                              <  8:137007/147‑1 (MQ=255)                                                                                                                     ggctTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAAGGGTCgg                           <  7:294848/149‑1 (MQ=255)                                                                                                                           tGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAAGGGTCGGTTAt                       >  6:253673/1‑147 (MQ=255)                                                                                                                             tCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAAGGGTCGGTTATCa                     <  6:128757/147‑1 (MQ=255)                                                                                                                             tCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAAGGGTCGGTTATCa                     <  4:325194/147‑1 (MQ=255)                                                                                                                                gca‑gAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGTTCCGGTGGGCCTTGTTTCTGCGCGCCTCGaaa                                  <  1:250770/131‑1 (MQ=255)                                                                                                                                gca‑gAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAAGGGTCGGTTATCAgtgt                 >  8:186656/1‑148 (MQ=255)                                                                                                                                gca‑gAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGATCCGGTGGGCCTTGTTTCTGCGCGCCTCGaaa                                  >  2:250770/1‑131 (MQ=255)                                                                                                                                    gagCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAAGGGTCGGTTATCAGTGTAAtc             >  2:254460/1‑149 (MQ=255)                                                                                                                                          gACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCACGCCTCGAAAGGGTCGGTTATCAGTGTAATCGTcgc        <  1:250365/148‑1 (MQ=255)                                                                                                                                           aCTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAAGGGTCGGTTATCAGTGTAATCGTCGCGa      <  4:185857/149‑1 (MQ=255)                                                                                                                                               gACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAAGGGTCGGTTATCAGTGTAATCGTCGCGATCct  >  1:15356/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                              GCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCA‑GAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAAGGGTCGGTTATCAGTGTAATCGTCGCGATCCT  >  NC_002947/4425032‑4425321

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