Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,033,557	A→G	10.9%	Q33R (CAG→CGG)	PP_2646 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,033,557	0	A	G	10.9%	142.0 / 6.8	55	Q33R (CAG→CGG)	PP_2646	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/22);  new base G (3/3);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACGCTTGCGCAACC  >  NC_002947/3033409‑3033703                                                                                                                                                      |                                                                                                                                                   cAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATCTGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgcc‑a                                                                                                                                                    <  3:56256/149‑1 (MQ=255)    gCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgcc‑agcg                                                                                                                                                 >  8:112887/1‑149 (MQ=255)    gCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgcc‑agcg                                                                                                                                                 >  7:75425/1‑149 (MQ=255)              cccTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCgg                                                                                                                                       <  3:314690/149‑1 (MQ=255)                ctctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGgt                                                                                                                                     >  5:227894/1‑149 (MQ=255)                  ctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCAGGACGAAGTGCAGCCCC‑AGCGGCCCATCCgcgggt                                                                                                                                   >  5:66307/1‑144 (MQ=255)                     tatCCGCCGTTTCGGCCAGT‑GTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGa                                                                                                                                <  4:196546/148‑1 (MQ=255)                        ccgccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATc                                                                                                                              >  1:345437/1‑148 (MQ=255)                         cgccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATcc                                                                                                                             <  8:75425/148‑1 (MQ=255)                          gccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGTGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCAc                                                                                                                           >  8:231160/1‑149 (MQ=255)                           ccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACg                                                                                                                          >  2:143782/1‑149 (MQ=255)                              tttCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTg                                                                                                                       >  1:334410/1‑149 (MQ=255)                                     cAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGc                                                                                                                  >  5:178599/1‑147 (MQ=255)                                        tGGTTCTGAGCGATGTAAGGAGAGCAGACCATGCGAGTACTGATATGTGAAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGtt                                                                                                             >  5:190530/1‑149 (MQ=255)                                                aGCGATGCAAGGAGAGCAGACCATGCGAGTACTTATATGTGCAGGGCGCAACTACGCCTACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                                        <  2:208873/146‑1 (MQ=255)                                                aGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtggt                                                                                                      <  1:339994/148‑1 (MQ=255)                                                aGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                                        <  1:143782/146‑1 (MQ=255)                                                 gCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtggtgg                                                                                                    <  5:299176/149‑1 (MQ=255)                                                 gCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtggtgg                                                                                                    >  6:299176/1‑149 (MQ=255)                                                                 aGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTgg                                                                                    <  7:112887/149‑1 (MQ=255)                                                                   aCCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGaaa                                                                                         <  5:102314/142‑1 (MQ=255)                                                                   aCCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGaaa                                                                                         >  6:102314/1‑142 (MQ=255)                                                                   aCCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGgc                                                                                  >  2:330328/1‑149 (MQ=255)                                                                       tGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGgcgcgc                                                                              <  6:257226/149‑1 (MQ=255)                                                                       tGCGACTACTGATATGTGCACGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGTAAATCGAAAGCTGGgcgcgc                                                                              <  6:190530/149‑1 (MQ=255)                                                                          gAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAg                                                                           >  4:103270/1‑149 (MQ=255)                                                                           aGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAg                                                                           >  4:348246/1‑148 (MQ=255)                                                                               cTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATgg                                                                      >  4:26296/1‑149 (MQ=255)                                                                                tGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCCGCGCC‑AGCGGCCGATCCGGGTGGTGAGCCACGGTGGCCGCCAGTTTGTTGGTGGGGAAATCGAGCGCTGGGCGCGCGCGCAGGgc                                                                     >  2:336288/1‑149 (MQ=255)                                                                                tGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCt                                                                                                          >  7:295652/1‑112 (MQ=255)                                                                                tGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCt                                                                                                          <  8:295652/112‑1 (MQ=255)                                                                                  atatGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGgcgc                                                                   >  3:159920/1‑149 (MQ=255)                                                                                  atatGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGgcgc                                                                   >  7:60926/1‑149 (MQ=255)                                                                                  atatGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCCATCCGGGTGTTGATCCACGGTGGCAGCCACTTTATTTGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGgcgc                                                                   >  6:30832/1‑149 (MQ=255)                                                                                  atatGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCCGCCGATCCGGGTGTTGATACACGGTGGCAGCCAGTTTCTTGGTGTGGAAATCGAAAGCTGGGCGCGCGAGCATGgcgc                                                                   >  8:131515/1‑149 (MQ=255)                                                                                   tatGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGcc                                                                  <  2:345437/149‑1 (MQ=255)                                                                                     tgtgCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                                        >  3:120898/1‑109 (MQ=255)                                                                                     tgtgCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                                        <  4:120898/109‑1 (MQ=255)                                                                                       tgCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATc                                                              >  7:174792/1‑149 (MQ=255)                                                                                         cAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTg                                                            <  6:227894/149‑1 (MQ=255)                                                                                            ggCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCAACGGTGGCAGCCAGTGTATTAGTGAGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTg                                                         >  4:329202/1‑149 (MQ=255)                                                                                                 aaCTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTa                                                    >  5:24001/1‑149 (MQ=255)                                                                                                   cTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTAcc                                                  <  1:336288/149‑1 (MQ=255)                                                                                                      cGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCg                                                <  3:26296/148‑1 (MQ=255)                                                                                                          gACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCAACGGTGGCAGCCAATTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCaaa                                           >  6:108874/1‑149 (MQ=255)                                                                                                                    gcTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTCCCCGCCAAACTGGCAGCTg                                 >  5:29753/1‑149 (MQ=255)                                                                                                                       gccgccAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCAGGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCAc                              >  5:329662/1‑149 (MQ=255)                                                                                                                        ccgccAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCTGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGTGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCa                               <  2:334410/147‑1 (MQ=255)                                                                                                                        ccgccAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACg                             <  5:30832/149‑1 (MQ=255)                                                                                                                        ccgccAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGGGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACg                             <  7:231160/149‑1 (MQ=255)                                                                                                                          cgctggCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACgg                            >  7:265498/5‑149 (MQ=37)                                                                                                                           gctggCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACg                             <  8:102512/144‑1 (MQ=37)                                                                                                                           gctggCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACg                             >  7:102512/4‑147 (MQ=37)                                                                                                                            ctggCGCTGGACGTCGTCCAGGGCCTGGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGa                                                                >  7:346140/3‑111 (MQ=37)                                                                                                                            ctggCGCTGGACGTCGTCCAGGGCCTGGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGa                                                                <  8:346140/109‑1 (MQ=37)                                                                                                                                     gacgtcGTCCAGGGCCTGGCGGCAGCGCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGc                   <  8:265498/147‑1 (MQ=38)                                                                                                                                      acgacgTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGaa               <  3:348246/149‑1 (MQ=255)                                                                                                                                           gTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGAGCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACGCtt          <  5:89336/149‑1 (MQ=255)                                                                                                                                             ccagcgcc‑agcgGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACGCTTgc        <  3:103270/149‑1 (MQ=255)                                                                                                                                               agcgcc‑agcgGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACGCTTgcgc      <  6:66307/149‑1 (MQ=255)                                                                                                                                                    cgagcgGCCGTTCCGGGTGTTGATCCACGTTGGCAGCAAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCAGATCTGGTGCGTTACCCGCCAAACTGTCAGCTGCACGGCAAGCAGGCCGAACGCTTGCGCAAcc  <  8:174792/147‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACGCTTGCGCAACC  >  NC_002947/3033409‑3033703

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