Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,951,401	Δ1 bp	5.9%	coding (496/666 nt)	PP_3487 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,951,401	0	G	.	5.9%	222.1 / 7.0	68	coding (496/666 nt)	PP_3487	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/29);  new base . (2/2);  total (37/31)
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.78e-01

AACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAGGTCACCTGGGCGGCCAGCGCTGG  >  NC_002947/3951256‑3951548                                                                                                                                                  |                                                                                                                                                    aaCAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCTACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTc                                                                                                                                                  <  1:117036/149‑1 (MQ=255)    aGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGa                                                                                                                                                >  7:65572/1‑148 (MQ=255)      cTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACAt                                                                                                                                             >  7:49008/1‑149 (MQ=255)         cATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTgg                                                                                                                                          <  7:298460/149‑1 (MQ=255)         cATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGCGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTg                                                                                                                                           <  1:267955/148‑1 (MQ=255)          aTGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGc                                                                                                                                         <  5:33565/149‑1 (MQ=255)                    cAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCg                                                                                                                               >  3:54746/1‑149 (MQ=255)                       tCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGt                                                                                                                             >  4:3687/1‑148 (MQ=255)                       tCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTc                                                                                                                            <  8:65572/149‑1 (MQ=255)                        cgccgcTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGt                                                                                                                             >  1:154373/1‑147 (MQ=255)                        cgccgcTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCA‑GTCGACATTGGCCTTGTCGCCGGTCAg                                                                                                                          >  1:367225/1‑149 (MQ=255)                         gccgcTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAg                                                                                                                          >  4:170552/1‑149 (MQ=255)                         gccgcTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAg                                                                                                                          <  4:54746/149‑1 (MQ=255)                          ccgcTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAgg                                                                                                                         <  3:170552/149‑1 (MQ=255)                                       aTCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCg                                                                                                            >  6:39619/1‑149 (MQ=255)                                           gggCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGcc                                                                                                         >  1:255872/1‑148 (MQ=255)                                               gAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATc                                                                                                     >  7:104496/1‑148 (MQ=255)                                                      gCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCA‑GTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCt                                                                                             >  1:138180/1‑148 (MQ=255)                                                       cTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCA‑GTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGc                                                                                           <  2:367225/149‑1 (MQ=255)                                                          cAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCtt                                                                                         >  7:274134/1‑149 (MQ=255)                                                          cAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCtt                                                                                         >  2:41328/1‑149 (MQ=255)                                                          cAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCt                                                                                          <  1:363628/148‑1 (MQ=255)                                                                  aCCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCCGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGgcggcg                                                                                  >  7:350776/1‑148 (MQ=255)                                                                    cTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAAc                                                                               >  6:361192/1‑149 (MQ=255)                                                                      gCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGc                                                                             >  5:370546/1‑149 (MQ=255)                                                                      gCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTCTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGCCGAACGc                                                                             >  8:75724/1‑149 (MQ=255)                                                                      gCAACATAACCGGCAATCAGCTGCGTCTAAAGCTGCTCGGGACTCGGCAGGAACAGCCCTATACCGCCAAGCACCAGGTCGCCATTGGCCTTTTCGCCGGTCAGGAACAGACAGTTCGGCCCATCCACCCCCTGCTTGGAGGCGAACGc                                                                             >  1:228251/1‑149 (MQ=255)                                                                       cAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCgg                                                                                                           >  3:121046/1‑118 (MQ=255)                                                                       cAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCgg                                                                                                           <  4:121046/118‑1 (MQ=255)                                                                              aCCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAgc                                                                      <  1:41328/148‑1 (MQ=255)                                                                                  gCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCtt                                                                          <  5:15355/140‑1 (MQ=255)                                                                                  gCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGCTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCtt                                                                          >  6:15355/1‑140 (MQ=255)                                                                                   cAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCg                                                                <  2:184603/149‑1 (MQ=255)                                                                                      tCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCgggg                                                             <  5:361192/149‑1 (MQ=255)                                                                                      tCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCgggg                                                             <  4:318896/149‑1 (MQ=255)                                                                                         gCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGgtg                                                          <  8:49008/149‑1 (MQ=255)                                                                                              gTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGtt                                                     >  2:245863/1‑149 (MQ=255)                                                                                                 gAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAg                                                  >  5:349764/1‑149 (MQ=255)                                                                                                   aTGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGc                                                 <  2:255872/148‑1 (MQ=255)                                                                                                   aTGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCg                                                <  8:104496/149‑1 (MQ=255)                                                                                                    tgctgcTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCgg                                               >  3:150006/1‑149 (MQ=255)                                                                                                       tgctCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGt                                             >  1:28247/1‑148 (MQ=255)                                                                                                           cGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGccc                                                                                                        <  3:154343/85‑1 (MQ=255)                                                                                                           cGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGccc                                                                                                        >  4:154343/1‑85 (MQ=255)                                                                                                            gggACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCg                                                       >  4:329596/1‑133 (MQ=255)                                                                                                            gggACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCg                                                       <  3:329596/133‑1 (MQ=255)                                                                                                            gggACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCActgc                                       >  1:255127/1‑149 (MQ=255)                                                                                                              gACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCggg                                                              >  4:295610/1‑124 (MQ=255)                                                                                                              gACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCggg                                                              <  3:295610/124‑1 (MQ=255)                                                                                                                cTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGAt                                   <  8:274134/149‑1 (MQ=255)                                                                                                                 tCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATc                                  >  1:181041/1‑149 (MQ=255)                                                                                                                   ggCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATc                                  <  6:349764/147‑1 (MQ=255)                                                                                                                      aGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGtt                                                     <  1:317303/125‑1 (MQ=255)                                                                                                                      aGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGtt                                                     >  2:317303/1‑125 (MQ=255)                                                                                                                      aGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGaca                              <  1:368026/148‑1 (MQ=255)                                                                                                                      aGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGaca                              >  2:368026/1‑148 (MQ=255)                                                                                                                        gAACTGCCCTATACGGCCAAGCACCA‑GTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGaa                          <  2:138180/149‑1 (MQ=255)                                                                                                                          aCTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAg                         <  6:370546/149‑1 (MQ=255)                                                                                                                          aCTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAg                         <  2:255127/149‑1 (MQ=255)                                                                                                                              cccTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAGGTCa                     >  6:274590/1‑149 (MQ=255)                                                                                                                              cccTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAGGTCa                     >  6:230478/1‑149 (MQ=255)                                                                                                                               ccTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAGGTCAc                    >  6:305046/1‑149 (MQ=255)                                                                                                                                cTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATc                                  <  8:6696/134‑1 (MQ=255)                                                                                                                                      ggCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAGGTCACCTGggcg             >  2:197173/1‑149 (MQ=255)                                                                                                                                          aaGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAg                         >  4:243444/1‑133 (MQ=255)                                                                                                                                          aaGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAg                         <  3:243444/133‑1 (MQ=255)                                                                                                                                              aCCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAGGTCACCTGGGCGGCCAgcgc     >  6:40779/1‑149 (MQ=255)                                                                                                                                                 aGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAGGTCACCTGGGCGGCCAGCGCTgg  <  3:3687/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCATCCACCCCCTGCTTGGCGGCGAACGCCTTCAGCGCCTCGGGGGTGTCGTTCAGCGGGTCACTGCTGATCGACACGAAGGTCACCTGGGCGGCCAGCGCTGG  >  NC_002947/3951256‑3951548

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