Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,855,539	T→G	8.6%	I106M (ATT→ATG)	PP_4270 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,855,539	0	T	G	8.6%	171.7 / 7.8	58	I106M (ATT→ATG)	PP_4270	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/24);  new base G (2/3);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

AAGTGGGCGCTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGTGA  >  NC_002947/4855391‑4855685                                                                                                                                                     |                                                                                                                                                      aaGTGGGCGCTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGAtt                                                                                                                                                       <  8:210824/149‑1 (MQ=255)      ggCGCTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATtcatca                                                                                                                                                  <  2:254083/149‑1 (MQ=255)        cgcTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAgc                                                                                                                                                <  5:213393/149‑1 (MQ=255)          cTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATGAATCAgcgg                                                                                                                                              <  8:160091/149‑2 (MQ=255)            cGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATGAATCAGCGGGa                                                                                                                                            <  6:189111/149‑1 (MQ=255)             gCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATGAATCAgcggg                                                                                                                                             <  7:203755/147‑3 (MQ=255)               cccGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACg                                                                                                                                          >  7:308514/1‑148 (MQ=255)                                    cgcctgcgcctgGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCcgacga                                                                                                                    <  1:229387/149‑1 (MQ=255)                                            cctgGAGGATGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCTAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCAccc                                                                                                            <  1:14898/149‑1 (MQ=255)                                                gaggagGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGAt                                                                                                        >  1:170603/1‑149 (MQ=255)                                                         aTGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGa                                                                                               <  8:308514/149‑1 (MQ=255)                                                                gcggcgATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTg                                                                                        >  6:176463/1‑149 (MQ=255)                                                                    cgATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCAc                                                                                    >  6:52138/1‑149 (MQ=255)                                                                    cgATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCAc                                                                                    >  6:52069/1‑149 (MQ=255)                                                                        atTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCAcccc                                                                                 <  5:52069/148‑1 (MQ=255)                                                                             tgtgCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGa                                                                        <  5:52138/149‑1 (MQ=255)                                                                               tgCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATc                                                                      <  2:198346/149‑1 (MQ=255)                                                                               tgCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATc                                                                      <  5:176463/149‑1 (MQ=255)                                                                                gCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCg                                                                     >  7:204741/1‑149 (MQ=255)                                                                                     gCTAAGTCAGTACCCAGTGGCTCAACGCACGACCAAGGGATTTTTTGTGGAGTATACCGCTGATTCATAAGCGCGACGTACTGGTTGGCCCCCCCCGCGACAGCCACCCGCTGCCCCCCGAGCAACTGCCACACCAA‑‑‑GCGATCGGCa                                                                  >  4:294712/1‑147 (MQ=255)                                                                                          gTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGc                                                                   <  6:4527/141‑1 (MQ=255)                                                                                          gTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGc                                                                   >  5:4527/1‑141 (MQ=255)                                                                                            cAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTCCTgg                                                         >  8:266530/1‑149 (MQ=255)                                                                                             aGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTggg                                                        <  2:170603/149‑1 (MQ=255)                                                                                                  ccAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTaca                                                   >  1:138183/1‑149 (MQ=255)                                                                                                    aGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACAcc                                                 >  4:160094/1‑149 (MQ=255)                                                                                                        gCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGCGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTAtt                                             >  2:320530/1‑149 (MQ=255)                                                                                                         cTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTc                                            <  4:157333/149‑1 (MQ=255)                                                                                                          tCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCg                                                                                                                        >  5:255429/1‑75 (MQ=255)                                                                                                          tCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCg                                                                                                                        <  6:255429/75‑1 (MQ=255)                                                                                                          tCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTc                                            >  2:107327/1‑148 (MQ=255)                                                                                                          tCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTc                                            >  5:59719/1‑148 (MQ=255)                                                                                                               gACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGc                                                       >  2:353561/1‑132 (MQ=255)                                                                                                               gACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGc                                                       <  1:353561/132‑1 (MQ=255)                                                                                                                 ccGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTg                                     >  2:225666/1‑148 (MQ=255)                                                                                                                  cGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGc                                    >  6:324501/1‑148 (MQ=255)                                                                                                                      cAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCa                               <  1:225666/149‑1 (MQ=255)                                                                                                                      cAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCa                               <  2:138183/149‑1 (MQ=255)                                                                                                                       aaGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGACGCCCCGGCGACAGCACCCCGATGAACCCGGAGCACCTGCCACCCCAG‑‑‑GAGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCAc                              >  5:158447/1‑149 (MQ=255)                                                                                                                        aGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAg                                                                                                                                                 >  4:130670/1‑36 (MQ=255)                                                                                                                        aGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAg                                                                                                                                                 <  3:130670/36‑1 (MQ=255)                                                                                                                        aGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGAt                                                                                                        <  8:87783/77‑1 (MQ=255)                                                                                                                        aGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGAt                                                                                                        >  7:87783/1‑77 (MQ=255)                                                                                                                        aGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCAct                             <  8:204741/149‑1 (MQ=255)                                                                                                                        aGGGAATTTGTGTGGATTATCCCGCTGATTCATCAGATCGACGTACTGGTTGTCCGACCCGGCGACAGCACCCCGATGCCCCCGGACCAACTGCCACCCCAT‑‑‑GCCAACGGCACGGTGCAGGGCTACACCTATTCAACCCTGCAACCAct                             >  8:94588/1‑149 (MQ=255)                                                                                                                         gggATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActg                            >  7:183332/1‑149 (MQ=255)                                                                                                                         gggATTTTGTGTGGAGAATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActg                            >  4:180967/1‑149 (MQ=255)                                                                                                                          ggATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActgc                           <  6:59719/149‑1 (MQ=255)                                                                                                                            aTTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActgctg                         <  8:183332/149‑1 (MQ=255)                                                                                                                                    tgGAGTATCGCGCTGATTCATCATCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCTATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTTCAACCACTGCTGGACCACg                  <  6:158447/148‑1 (MQ=255)                                                                                                                                      gAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGc                >  7:95384/1‑148 (MQ=255)                                                                                                                                      gAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCACCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGc                >  6:28073/1‑148 (MQ=255)                                                                                                                                      gAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGAAGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCAATTCAACCCTGCAACCACTGCTCGACCACGGc                >  6:22833/1‑148 (MQ=255)                                                                                                                                         tATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGt            <  7:266530/149‑1 (MQ=255)                                                                                                                                          gtcgcGCTGATGAATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGtt           >  5:328770/5‑149 (MQ=255)                                                                                                                                             gcGCTGATGAATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCa        >  4:226823/2‑149 (MQ=255)                                                                                                                                             ccGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCa        >  8:157835/1‑149 (MQ=255)                                                                                                                                               gCTGATTCATCAGCGCGACGTACAGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGCTACCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGc         <  1:107327/149‑1 (MQ=255)                                                                                                                                                 tGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGt    >  5:350231/1‑149 (MQ=255)                                                                                                                                                   aTTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGTGa  <  3:160094/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      AAGTGGGCGCTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGTGA  >  NC_002947/4855391‑4855685

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