Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,276,382	T→C	7.4%	E105G (GAG→GGG)	PP_2875 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,276,382	0	T	C	7.4%	152.7 / 2.4	54	E105G (GAG→GGG)	PP_2875	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base C (2/2);  total (27/27)
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.29e-01

CGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTGGGCAGGTGCTCTTCAAATT  >  NC_002947/3276236‑3276529                                                                                                                                                   |                                                                                                                                                    cGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTcg                                                                                                                                                   >  2:92279/1‑149 (MQ=255)    gCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCgg                                                                                                                                                <  4:264406/149‑1 (MQ=255)         ggCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGc                                                                                                                                           >  2:87005/1‑149 (MQ=255)          gCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCg                                                                                                                                          >  1:65314/1‑149 (MQ=255)            aaCAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCCGCTGCCCGCGAGCAGCagg                                                                                                                                        >  1:50762/1‑149 (MQ=255)              cAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCCGCTGCCCGCGAGCAGCAGGca                                                                                                                                      <  2:50762/149‑2 (MQ=255)                 gcgcCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCAtt                                                                                                                                   >  8:162523/1‑149 (MQ=255)                 gcgcCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCAt                                                                                                                                    >  6:45499/1‑148 (MQ=255)                    ccTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCtt                                                                                                                                >  4:14852/1‑149 (MQ=255)                        gCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCttgttg                                                                                                                            >  2:156021/1‑149 (MQ=255)                        gCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCttgt                                                                                                                              <  1:126183/147‑1 (MQ=255)                          ctctGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCa                                                                                                                          <  2:348313/149‑1 (MQ=255)                             tGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCCGCTGCCCGCGAGCAGCAGGca                                                                                                                                      <  7:326010/134‑2 (MQ=255)                             tGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCCGCTGCCCGCGAGCAGCAGGca                                                                                                                                      >  8:326010/1‑133 (MQ=255)                                 ggAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCttgt                                                                                                                              >  5:25450/1‑138 (MQ=255)                                 ggAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCttgt                                                                                                                              <  6:25450/138‑1 (MQ=255)                                  gAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTgc                                                                                                                  <  3:14852/149‑1 (MQ=255)                                   aaaCCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTgca                                                                                                                 <  7:162523/149‑1 (MQ=255)                                     aCCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTgcagc                                                                                                               <  5:45499/149‑1 (MQ=255)                                       ccGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCAc                                                                                                             >  6:360332/1‑149 (MQ=255)                                                      ccAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTc                                                                                              <  1:87005/149‑1 (MQ=255)                                                         ggCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAgt                                                                                           >  7:88872/1‑149 (MQ=255)                                                         ggCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAgt                                                                                           <  7:99164/149‑1 (MQ=255)                                                               ccagaGCCCGGCCTGCATAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTcggc                                                                                     >  2:329274/5‑149 (MQ=255)                                                               acaaaGCCCAGCCTATCCAGGCCGTGCACTGGTACGCCTGCCAACCCACAATTCACCGGTTCGATCACCTTGACTGCCTGCTCCTCGCGGGCAGCGGTCCATTCTTGATGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTcgc                                                                                      <  8:151206/146‑2 (MQ=255)                                                               aaaaaGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAgt                                                                                           >  1:143869/1‑143 (MQ=255)                                                                  aaGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTcggcggc                                                                                  <  5:360332/149‑1 (MQ=255)                                                                       cGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCa                                                                                             >  8:277881/1‑133 (MQ=255)                                                                       cGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCa                                                                                             <  7:277881/133‑1 (MQ=255)                                                                           cTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAgtgt                                                                                         >  6:59713/1‑133 (MQ=255)                                                                           cTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAgtgt                                                                                         <  5:59713/133‑1 (MQ=255)                                                                                            tGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTg                                                          <  1:156021/147‑1 (MQ=255)                                                                                            tGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTAAGTGTCGGCGGCGGCATTc                                                                           >  2:306136/1‑130 (MQ=255)                                                                                              gTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTcggcggcg                                                                                 <  4:354795/122‑1 (MQ=255)                                                                                              gTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTcggcggcg                                                                                 >  3:354795/1‑122 (MQ=255)                                                                                                    gTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCAcgg                                                   <  4:32317/146‑1 (MQ=255)                                                                                                    gTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCAcgg                                                   >  3:32317/1‑146 (MQ=255)                                                                                                       cGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTa                                              <  2:65314/148‑1 (MQ=255)                                                                                                                     aGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTgg                                                         >  7:191271/1‑123 (MQ=255)                                                                                                                     aGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCggg                                <  5:308815/148‑1 (MQ=255)                                                                                                                     aGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGAAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTgg                                                         <  8:191271/123‑1 (MQ=255)                                                                                                                            cGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCAt                         <  8:88872/148‑1 (MQ=255)                                                                                                                              aTCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCt                      >  6:245651/1‑149 (MQ=255)                                                                                                                               tCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCCCGCATGCTg                     >  8:252565/1‑149 (MQ=255)                                                                                                                                  ccTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCtt                                                           >  8:122608/1‑108 (MQ=255)                                                                                                                                  ccTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCtt                                                           <  7:122608/108‑1 (MQ=255)                                                                                                                                  ccTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTggg                   >  5:260182/1‑148 (MQ=255)                                                                                                                                   cTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTggg                   >  8:232361/1‑147 (MQ=255)                                                                                                                                   cTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTggg                   >  1:168926/1‑147 (MQ=255)                                                                                                                                      gACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTGGGCAg                <  2:168926/147‑1 (MQ=255)                                                                                                                                       aCTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTGGGCAGGt              <  6:260182/148‑1 (MQ=255)                                                                                                                                         tgcctgcTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGTAGCGGGCACGCATGCTGGGCAGGTg             <  4:70829/147‑1 (MQ=255)                                                                                                                                         tgcctgcTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGTAGCGGGCACGCATGCTGGGCAGGTg             >  3:70829/1‑147 (MQ=255)                                                                                                                                           cctgcTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTGGGCAGGTGctc          <  5:141535/148‑1 (MQ=255)                                                                                                                                              gctgctCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTGGGCAGGTGCTCTTc       >  2:299406/1‑148 (MQ=255)                                                                                                                                               ctgctCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTGGGCAGGTGCTCTTc       <  5:245651/147‑1 (MQ=255)                                                                                                                                               ctgctCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTGGGCAGGTGCTCTTc       <  1:329274/147‑1 (MQ=255)                                                                                                                                                   tCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTGGGCAGGTGCTCTTCa      <  2:304515/144‑1 (MQ=255)                                                                                                                                                   tCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTGGGCAGGTGCTCTTCa      >  1:304515/1‑144 (MQ=255)                                                                                                                                                   tCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTGGGCAGGTGCTCTTCAAAtt  >  3:140919/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATGCAGGCCGGGCTCTGGTAAGGGTGCGAACCCTCAATTCAGCGGTTCGATCACCTTGACTGCCTGCTGCTCGCGGGCAGCGGGCCATTCTTGTTGCAGGGTCTGCAGCACGCGGCCAACGAATTCAGTGTCGGCGGCGGCTTTCTTGCCGACGTAACCTTGGCCACGGCGGTAGACCTTGAAGCGGGCACGCATGCTGGGCAGGTGCTCTTCAAATT  >  NC_002947/3276236‑3276529

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