Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,161,481	G→A	5.7%	L80L (CTG→CTA)	PP_5741 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,161,481	0	G	A	5.7%	211.4 / 2.8	70	L80L (CTG→CTA)	PP_5741	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/35);  new base A (2/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

ACGTTCGCGGTTATCTGTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCACTGGTCGTGTAGAAACGG  >  NC_002947/6161339‑6161626                                                                                                                                               |                                                                                                                                                   tcGTTCGCTGTTATCTGTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGTAGCGTTATACCTCCCACCTGCTATTc                                                                                                                                              <  1:319274/148‑1 (MQ=255)      cgcgGTTATCTGTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGcc                                                                                                                                          <  7:159057/148‑1 (MQ=255)      cgcgGTTATCTGTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGcc                                                                                                                                          <  6:115302/148‑1 (MQ=255)       gcgGTTATCTGTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCt                                                                                                                                        >  4:171992/1‑149 (MQ=255)          gTTATCTGTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCGACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTggc                                                                                                                                     >  8:176201/1‑149 (MQ=255)             aTCTGTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTggcgg                                                                                                                                   <  2:164892/148‑1 (MQ=255)               cTGTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGcc                                                                                                                                 >  7:260029/1‑148 (MQ=255)                tGTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGa                                                                                                                               <  3:24122/149‑1 (MQ=255)                 gTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGccgcca                                                                                                                                      >  4:119576/1‑137 (MQ=255)                 gTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGccgcca                                                                                                                                      <  3:119576/141‑5 (MQ=255)                 gTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGAc                                                                                                                              >  4:99297/1‑149 (MQ=255)                       agagTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTcgc                                                                                                                        >  5:216997/1‑149 (MQ=255)                        gagTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTcgcc                                                                                                                       >  2:181737/1‑149 (MQ=255)                        gagTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTcgcc                                                                                                                       >  4:60325/1‑149 (MQ=255)                        gagTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTcgcc                                                                                                                       >  6:227179/1‑149 (MQ=255)                         cgTTGGTGAAGCTGGGCCAGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGACGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTcgccg                                                                                                                      <  8:353557/148‑1 (MQ=255)                         agTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTcgccg                                                                                                                      >  3:165512/1‑149 (MQ=255)                          gTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTcgccg                                                                                                                      >  8:75501/1‑148 (MQ=255)                             ggTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑g                                                                                                                  <  5:227179/148‑1 (MQ=255)                             ggTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑g                                                                                                                  <  3:60325/148‑1 (MQ=255)                                    cTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATAcc                                                                                                           >  3:110310/1‑148 (MQ=255)                                           tGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGgc                                                                                                   >  5:108693/1‑149 (MQ=255)                                            gAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGgcg                                                                                                  >  6:138975/1‑149 (MQ=255)                                              aaaaTCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCa                                                                                                <  4:110310/149‑1 (MQ=255)                                                 aTCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCatat                                                                                             <  6:377034/149‑1 (MQ=255)                                                    gggCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGa                                                                                          >  5:53726/1‑149 (MQ=255)                                                     ggCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGAt                                                                                         <  5:256843/149‑1 (MQ=255)                                                        cacaGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCt                                                                                      >  4:336899/1‑149 (MQ=255)                                                         acaGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTc                                                                                     >  7:167240/1‑149 (MQ=255)                                                                 gCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAgg                                                                              <  3:148766/148‑1 (MQ=255)                                                                 gCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAgg                                                                              >  4:148766/1‑148 (MQ=255)                                                                 gCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGc                                                                             <  8:167240/149‑1 (MQ=255)                                                                  cTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTg                                                                                           >  7:107071/1‑134 (MQ=255)                                                                  cTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTg                                                                                           <  8:107071/134‑1 (MQ=255)                                                                      cTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAAc                                                                        >  3:343421/1‑149 (MQ=255)                                                                      cTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAAc                                                                        >  5:336280/1‑149 (MQ=255)                                                                      cTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTCGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGATTACCGAAACGGCGCATATTGATGCTCGCAAAGGCCCAAc                                                                        >  5:159400/1‑149 (MQ=255)                                                                       tCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAAct                                                                       >  6:183453/1‑149 (MQ=255)                                                                          gcgcgcCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTg                                                                    <  5:183453/149‑1 (MQ=255)                                                                           cgcgcCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTgg                                                                   <  7:75501/149‑1 (MQ=255)                                                                            gcgcCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGc                                                                  <  3:336899/149‑1 (MQ=255)                                                                               ccTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCgg                                                                                                                                            >  2:358106/1‑73 (MQ=255)                                                                               ccTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCgg                                                                                                                                            <  1:358106/73‑1 (MQ=255)                                                                                   ctacGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAAt                                                             <  6:108693/147‑1 (MQ=255)                                                                                      cGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAAtt                                                        <  3:99297/149‑1 (MQ=255)                                                                                         gggATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTc                                                     <  8:260029/149‑1 (MQ=255)                                                                                              gCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTACTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAAc                                                >  6:24134/1‑149 (MQ=255)                                                                                              gCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTACTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAAc                                                >  6:18862/1‑149 (MQ=255)                                                                                                   cagcagCAAGCCACAGCTGTTCTGGAGCGTTATACCTCCCACCTGCTCTTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTg                                            <  4:365982/148‑1 (MQ=255)                                                                                                   cAGAAGGAAGCCATAGCTTTTCTGGAGCTTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGTCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTg                                            <  6:216997/148‑1 (MQ=255)                                                                                                   cAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTg                                                                    >  8:82903/1‑124 (MQ=255)                                                                                                   cAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTg                                                                    <  7:82903/124‑1 (MQ=255)                                                                                                    aGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTgg                                           <  3:171992/148‑1 (MQ=255)                                                                                                    aGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGc                                          <  6:53726/149‑1 (MQ=255)                                                                                                     gaaggaagCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCt                                                                                      <  2:204561/104‑1 (MQ=255)                                                                                                     gaaggaagCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCt                                                                                      >  1:204561/1‑104 (MQ=255)                                                                                                           agCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCt                                   <  7:176201/149‑1 (MQ=255)                                                                                                            gCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTg                                                                                           <  2:247217/92‑1 (MQ=255)                                                                                                            gCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTg                                                                                           >  1:247217/1‑92 (MQ=255)                                                                                                            gCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCtt                                  >  6:184941/1‑149 (MQ=255)                                                                                                                      ttCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTtgatgatga                          <  5:138975/147‑1 (MQ=255)                                                                                                                             gCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCACTg                 <  1:181737/149‑1 (MQ=255)                                                                                                                              cGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCACTgg                <  4:343421/149‑1 (MQ=255)                                                                                                                               gTTATACCTCCCACCTACTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCACTGGt               <  5:24134/149‑1 (MQ=255)                                                                                                                               gTTATACCTCCCACCTACTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCACTGGt               <  5:18862/149‑1 (MQ=255)                                                                                                                                   taCCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCACTGGTCgt            <  6:336280/148‑1 (MQ=255)                                                                                                                                       tCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCACTGGTCGTGTAg        >  2:115282/1‑148 (MQ=255)                                                                                                                                       tCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCAATGGTCGTGTAg        >  2:120791/1‑148 (MQ=255)                                                                                                                                          cACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCACTGGTCGTGTAGAAAc    <  1:115282/149‑1 (MQ=255)                                                                                                                                              tGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCACTGGTCGTGTAGAAACgg  >  3:73233/1‑147 (MQ=255)                                                                                                                                               gCTATTCGGCCCTGGCGGCCGACTCTCGCCGCAGGGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCACTGGTCGTGTAGAAACgg  >  1:244641/1‑147 (MQ=255)                                                                                                                                               |                                                                                                                                                   ACGTTCGCGGTTATCTGTTATCAGAGTTGGTGAAGCTGTGCCTGAAAAATCGGGCCACAGTTCCGCTAGCTCGGCGCGCCTACTACGAGGGATGCGTTCAGAAGGAAGCCATAGCTTTTCTGGAGCGTTATACCTCCCACCTGCTATTCGGCCCTGGCGGCCGACTCTCGCCGCA‑GGAATACCGATACGGCGCATATTGATGCTCGCAAAGGCCCAACTCTGGCGTAATCAATTGTCCCAACGTTGGCAGAGGCTTGATGATGAAAATCACTGGTCGTGTAGAAACGG  >  NC_002947/6161339‑6161626

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