Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,918,256	A→G	6.6%	N219D (AAC→GAC)	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,918,256	0	A	G	6.6%	172.0 / 2.1	60	N219D (AAC→GAC)	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base G (2/2);  total (30/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 = 9.63e-01

GTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCT  >  NC_002947/4918112‑4918404                                                                                                                                                  |                                                                                                                                                     gTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGc                                                                                                                                                  >  4:126881/1‑149 (MQ=255)       ggCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGc                                                                                                                                             <  6:158956/148‑1 (MQ=255)        gCAGCCTGCTCAGCGCCCTGGTGCAGGGTGAGTCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCGACGCCCTGc                                                                                                                                             >  1:219365/1‑148 (MQ=255)        gCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCGACGCCCTGCaa                                                                                                                                           >  7:218316/1‑149 (MQ=255)                cTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCAc                                                                                                                                   <  4:193542/149‑1 (MQ=255)                    gcgcCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCGACGCCCTGCAACAGCTCACCcgc                                                                                                                               <  2:219365/149‑1 (MQ=255)                    gcgcCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCGACGCCCTGCAACAGCTCACCcgc                                                                                                                               <  8:218316/149‑1 (MQ=255)                         cTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGAccc                                                                                                                          <  7:102636/149‑1 (MQ=255)                                  ggTGA‑TCAGTCATTGGGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGa                                                                                                                 <  6:16572/149‑1 (MQ=255)                                    tGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATg                                                                                                               <  3:126881/149‑1 (MQ=255)                                      a‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATgg                                                                                                              <  5:42318/148‑1 (MQ=255)                                        tCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGc                                                                                                             <  3:89220/148‑1 (MQ=255)                                          aGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGAc                                                                                                          <  6:291556/149‑1 (MQ=255)                                              aTTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTCCAACAGCTCACCCGCCACCCAGGCGTCGATGCCGACCTgc                                                                                                      >  6:12942/1‑149 (MQ=255)                                                 gcgcTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTgcgc                                                                                                    <  5:12942/148‑1 (MQ=255)                                                 gcgcTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTgcgc                                                                                                    <  2:270888/148‑1 (MQ=255)                                                      gcctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGa                                                                                                                 <  7:274115/127‑1 (MQ=255)                                                      gcctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGa                                                                                                                 >  8:274115/4‑130 (MQ=255)                                                      gcctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAgc                                                                                              >  1:163317/4‑149 (MQ=255)                                                       cctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAgcg                                                                                             >  6:17814/3‑149 (MQ=255)                                                        gTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAgcg                                                                                             <  3:346796/148‑1 (MQ=255)                                                         tGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGcc                                                                                           >  7:338914/1‑149 (MQ=255)                                                          gCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGAc                                                                                                          >  1:207736/1‑133 (MQ=255)                                                          gCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGAc                                                                                                          <  2:207736/133‑1 (MQ=255)                                                               tgaGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGc                                                                                     <  5:17814/147‑1 (MQ=255)                                                                         gcaggACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGc                                                                             >  8:289319/1‑147 (MQ=255)                                                                                  aCTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCg                                                                   >  4:175713/1‑148 (MQ=255)                                                                                    tCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGAc                                                                >  4:294553/1‑149 (MQ=255)                                                                                    tCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGAc                                                                >  6:179518/1‑149 (MQ=255)                                                                                      aGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGATACCCGGCCTGCTCAACGCCCAGCAACAGCTCACACGCGACCCCGGCGTCGATGGCGACCTGCCCGCCAGCCCCGAACGCCTGCTGGCAAGCCTGCCGGACGc                                                              >  3:106092/1‑149 (MQ=255)                                                                                        cATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGccc                                                            <  8:338914/149‑1 (MQ=255)                                                                                        cATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGccc                                                            <  2:163317/149‑1 (MQ=255)                                                                                        cATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGccc                                                            >  1:51896/1‑149 (MQ=255)                                                                                         aTTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCg                                                           >  3:300078/1‑149 (MQ=255)                                                                                            tcccAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTGGCTACCGGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGAc                                                         <  6:68267/147‑1 (MQ=255)                                                                                              ccAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAAc                                                      >  7:171265/1‑149 (MQ=255)                                                                                              ccAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAAc                                                      >  1:293277/1‑149 (MQ=255)                                                                                               cAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACt                                                     >  6:355440/1‑149 (MQ=255)                                                                                                 aCAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTgg                                                   <  2:51896/149‑1 (MQ=255)                                                                                                   aGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGCCAACTGGGc                                                 >  6:242369/1‑149 (MQ=255)                                                                                                   aGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGc                                                 >  6:6499/1‑149 (MQ=255)                                                                                                     ctgctgACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGAAACAGCTCACCAGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGa                                               >  8:198354/1‑149 (MQ=255)                                                                                                      tgctgACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGAc                                              >  7:46735/1‑149 (MQ=255)                                                                                                            cccAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAg                                        <  8:171265/149‑1 (MQ=255)                                                                                                             ccAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAgg                                       >  2:21581/1‑149 (MQ=255)                                                                                                             ccAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGAGGGCGCCCAGAGCGCCCCCGCCGAACCCCAGCTCGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAgg                                       >  1:289793/1‑149 (MQ=255)                                                                                                               aGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCa                                     <  2:293277/149‑1 (MQ=255)                                                                                                                  aaaaCCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTg                                  <  5:179217/149‑1 (MQ=255)                                                                                                                     aCCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGcc                               <  5:179518/149‑1 (MQ=255)                                                                                                                      ccGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGccc                              <  5:355440/149‑1 (MQ=255)                                                                                                                       cGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGcc                               >  6:217853/1‑147 (MQ=255)                                                                                                                       cGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGcc                               >  7:211393/1‑147 (MQ=255)                                                                                                                       cGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCa                             <  8:211393/149‑1 (MQ=255)                                                                                                                           aaGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGcc                         <  8:46735/149‑1 (MQ=255)                                                                                                                             gCCTCGCTACCCGGGCTGCTGAACGCCCTGCAAGAGCTCACCCGGGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCt                       <  4:106092/149‑1 (MQ=255)                                                                                                                              ccTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTg                      <  7:289319/149‑1 (MQ=255)                                                                                                                                tCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGaa                    >  8:282422/1‑149 (MQ=255)                                                                                                                                       ccGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGc             >  6:310402/1‑149 (MQ=255)                                                                                                                                           ccTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCa         <  1:173573/149‑1 (MQ=255)                                                                                                                                           ccTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTCGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCgg           >  7:327707/1‑147 (MQ=255)                                                                                                                                                  aaCGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTtct  <  5:217853/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCT  >  NC_002947/4918112‑4918404

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