Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,491,493	T→C	6.5%	K309R (AAG→AGG)	PP_3096 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,491,493	0	T	C	6.5%	176.2 / 2.6	61	K309R (AAG→AGG)	PP_3096	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base C (2/2);  total (29/32)
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.30e-01

CACAAGGCGACCCCATCGTAAGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACT  >  NC_002947/3491349‑3491640                                                                                                                                                 |                                                                                                                                                     cacaAGGCGACCCCATCGTAAGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGa                                                                                                                                                   <  7:62723/148‑1 (MQ=255)       gCGACCCCATCGTAAGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCgg                                                                                                                                            <  3:232982/149‑1 (MQ=255)         gACCCCATCGTAAGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCCTCAAgcgcaagg                                                                                                                                          <  6:313853/149‑5 (MQ=255)           ccccATCGTAAGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑Ca                                                                                                                                       <  3:29821/149‑1 (MQ=255)            cccATCGTAAGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCACGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCCTCAAGcgcga                                                                                                                                            <  6:306260/144‑2 (MQ=255)            cccATCGTAAGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCACGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCCTCAAGcgcga                                                                                                                                            >  5:306260/1‑143 (MQ=255)              cATCGTAAGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAAtt                                                                                                                                    >  2:200440/1‑149 (MQ=255)                   tAAGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGaa                                                                                                                               <  8:8683/149‑1 (MQ=255)                    aaGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTAACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCCACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGaaa                                                                                                                              >  6:331575/1‑149 (MQ=255)                     aGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGaaa                                                                                                                              <  3:67987/148‑1 (MQ=255)                     aGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGaaa                                                                                                                              <  7:49776/148‑1 (MQ=255)                      ggAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATc                                                                                                                            <  3:162450/149‑1 (MQ=255)                       gAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATcc                                                                                                                           <  6:369123/149‑1 (MQ=255)                              gAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCgg                                                                                                                    >  5:135115/1‑149 (MQ=255)                              gAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCgg                                                                                                                    >  1:220063/1‑149 (MQ=255)                                aGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCa                                                                                                                  <  2:129950/149‑1 (MQ=255)                                      cgcAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGc                                                                                                                   <  3:236407/142‑1 (MQ=255)                                      cgcAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGc                                                                                                                   >  4:236407/1‑142 (MQ=255)                                       gcAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGTCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAAt                                                                                                                                     <  1:84162/123‑1 (MQ=255)                                       gcAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGTCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAAt                                                                                                                                     >  2:84162/1‑123 (MQ=255)                                           cgcgAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCAc                                                                                                                 >  5:67800/1‑139 (MQ=255)                                           cgcgAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCAc                                                                                                                 <  6:67800/139‑1 (MQ=255)                                           cgcgAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGaa                                                                                                       <  3:67190/149‑1 (MQ=255)                                           cgcgAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGaa                                                                                                       <  2:35724/149‑1 (MQ=255)                                             cgAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAgc                                                                                                     >  2:171463/1‑149 (MQ=255)                                              gAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAgcg                                                                                                    >  2:101317/1‑149 (MQ=255)                                                                  aGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCCACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGc                                                                                <  5:331575/149‑1 (MQ=255)                                                                   gggCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGACAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGc                                                                                <  1:101317/149‑1 (MQ=255)                                                                   gggCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCg                                                                               >  3:86761/1‑149 (MQ=255)                                                                   gcgCGTTAGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGATGGGCGAGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAAGTTCGAAATCCAGAACTGCACCTATATAGAAGCCCACAAACTCGACCACCTTGCg                                                                               >  8:161339/3‑149 (MQ=255)                                                                      cGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTa                                                                            <  1:171463/149‑1 (MQ=255)                                                                             gCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGc                                                                                    <  7:38726/134‑1 (MQ=255)                                                                             gCCAACCCAGCCAGCCAAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTAGGCAAATATATAGAAGCGCACCAACTCGACCAGc                                                                                    >  8:38726/1‑134 (MQ=255)                                                                                 aCCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGc                                                                 >  5:358972/1‑149 (MQ=255)                                                                                   ccagccagCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGccc                                                               >  3:134086/1‑149 (MQ=255)                                                                                   ccagccagCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGccc                                                               <  6:358972/149‑1 (MQ=255)                                                                                      gccagccCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCg                                                                               <  4:255101/130‑1 (MQ=255)                                                                                      gccagccCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCg                                                                               >  3:255101/1‑130 (MQ=255)                                                                                      gccagccCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGc                                                            >  2:25332/1‑149 (MQ=255)                                                                                       ccagccCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGccc                                                               >  4:336504/1‑145 (MQ=255)                                                                                       ccagccCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGccc                                                               <  3:336504/145‑1 (MQ=255)                                                                                       ccagccCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCt                                                           <  5:364533/149‑1 (MQ=255)                                                                                          cccaaGGGCTAAATTGCGGCTCGGGCCATGTTGGGGGGGTGGAACCGCTTCGCGCTGTAGGGGCGGGG‑AAATTTCTAAATCACGTTCGGCACCAAAAAAGAACCCCCCCAACTCTA‑CAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGaa                                                       <  2:311576/145‑1 (MQ=255)                                                                                              aGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGaa                                                       <  3:179814/146‑1 (MQ=255)                                                                                              aGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTAGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACGAACACGACCCGCATGCGGTAGCAGGCGGTGCCCGGCTGGaa                                                       >  4:179814/1‑146 (MQ=255)                                                                                               gggCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAgcgc                                                   <  1:200440/149‑1 (MQ=255)                                                                                                ggCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGt                                                                             >  4:22571/1‑122 (MQ=255)                                                                                                ggCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGt                                                                             <  3:22571/122‑1 (MQ=255)                                                                                                ggCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCt                                                  >  4:69967/1‑149 (MQ=255)                                                                                                            gCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCa                                      <  2:220063/149‑1 (MQ=255)                                                                                                             cTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTgg                                                >  1:4268/1‑138 (MQ=255)                                                                                                             cTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTgg                                                <  2:4268/138‑1 (MQ=255)                                                                                                                 ggCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCt                                  >  2:125021/1‑148 (MQ=255)                                                                                                                   cccAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAg                               >  8:13811/1‑149 (MQ=255)                                                                                                                       gTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCggg                           <  1:125021/149‑1 (MQ=255)                                                                                                                             cgcgTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGc                     >  6:288007/1‑149 (MQ=255)                                                                                                                                   gTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATg               >  8:75778/1‑149 (MQ=255)                                                                                                                                      ttgCCCCGCGCCTCAAGCGCGAGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCg             >  7:269127/3‑148 (MQ=255)                                                                                                                                      ccGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCg             >  2:8911/1‑148 (MQ=255)                                                                                                                                             gcgcTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCAc      >  5:236437/1‑148 (MQ=255)                                                                                                                                                cTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACt  <  1:25332/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CACAAGGCGACCCCATCGTAAGGAATGTGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGG‑CAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACT  >  NC_002947/3491349‑3491640

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