Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	790,957	T→A	13.5%	Q46L (CAG→CTG)	PP_0677 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	790,957	0	T	A	13.5%	137.3 / 13.0	59	Q46L (CAG→CTG)	PP_0677	lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (27/24);  new base A (2/6);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

AAGCGCTCGTGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAG  >  NC_002947/790810‑791095                                                                                                                                                    |                                                                                                                                           aaGCGCTCGTGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTg                                                                                                                                           >  6:136852/1‑149 (MQ=255)          tGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccggg                                                                                                                                  <  8:181754/149‑4 (MQ=255)          tGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccggg                                                                                                                                  <  5:303904/149‑4 (MQ=255)          tGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccggg                                                                                                                                  <  6:247560/149‑4 (MQ=255)          tGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccggg                                                                                                                                  <  2:351533/149‑4 (MQ=255)           ggTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCActg                                                                                                                                 >  3:367745/1‑149 (MQ=255)            gTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCTTGGCGGCCActg                                                                                                                                 <  7:109833/148‑1 (MQ=255)            gTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccgggc                                                                                                                                 <  3:7190/148‑5 (MQ=255)                      tGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTc                                                                                                                      >  7:154024/1‑149 (MQ=255)                        ttcttcGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTcc                                                                                                                     <  8:157194/148‑1 (MQ=255)                             cGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGTCCACTGCTGCAGCCGTCCTGGaa                                                                                                                <  5:112118/148‑1 (MQ=255)                             cGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAg                                                                                                               >  4:120656/1‑149 (MQ=255)                                     tGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCg                                                                                                        >  1:195769/1‑148 (MQ=255)                                      gCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCa                                                                                                      >  2:303402/1‑149 (MQ=255)                                       cATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAg                                                                                                               >  6:91154/1‑139 (MQ=255)                                       cATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAg                                                                                                               <  5:91154/139‑1 (MQ=255)                                        aTGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGACGTCCTGGAAGCCATCCGCa                                                                                                      >  6:288177/1‑147 (MQ=255)                                        aTGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCCTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATa                                                                                                    <  1:316990/149‑1 (MQ=255)                                                acacTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCaggag                                                                                            <  8:157445/149‑1 (MQ=255)                                                 cacTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCaggagg                                                                                           <  3:294984/149‑1 (MQ=255)                                                  acTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAggg                                                                                          >  7:223975/1‑149 (MQ=255)                                                  acTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGCAggg                                                                                          >  1:177418/1‑149 (MQ=255)                                                   cTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGCTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAggg                                                                                          >  3:105137/1‑148 (MQ=255)                                                        ggAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCCCTGCTGGAGCCGTCCTGGAAGCCATCCGCTGAGGCAGGAGGGtcgccg                                                                                    >  8:353289/1‑144 (MQ=255)                                                           aGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTgg                                                                                 <  1:333991/149‑1 (MQ=255)                                                            gCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGACGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGt                                                                                <  5:288177/149‑1 (MQ=255)                                                              cATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGtt                                                                               >  1:370537/1‑148 (MQ=255)                                                               aTCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGa                                                                             >  4:171793/1‑149 (MQ=255)                                                                   tCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCa                                                                          >  8:73060/1‑148 (MQ=255)                                                                             tttCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCgtgg                                                               <  4:367745/149‑1 (MQ=255)                                                                                 gTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCgtgggtg                                                            <  6:42676/148‑1 (MQ=255)                                                                                      gACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCg                                                      >  4:182165/1‑149 (MQ=255)                                                                                       aCGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTc                                                       <  5:136852/147‑1 (MQ=255)                                                                                         gCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCaa                                                   >  6:109650/1‑149 (MQ=255)                                                                                         gCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCaa                                                   >  7:366164/1‑149 (MQ=255)                                                                                           tCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAAcc                                                 >  8:113881/1‑149 (MQ=255)                                                                                             gTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAAcc                                                 <  5:109650/147‑1 (MQ=255)                                                                                              tGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGa                                              <  3:120656/149‑1 (MQ=255)                                                                                                    tAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCGTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGAc                                             <  1:171692/144‑1 (MQ=255)                                                                                                    tAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCGTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGAc                                             >  2:171692/1‑144 (MQ=255)                                                                                                         cgGTACGTTCTTGCGGAAGTTGCCCACCATGAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACg                                   <  2:177418/149‑1 (MQ=255)                                                                                                         cgGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACg                                   <  6:258073/149‑1 (MQ=255)                                                                                                         cgGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACg                                   <  2:370537/149‑1 (MQ=255)                                                                                                              cGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAAcc                                                 >  2:263356/1‑130 (MQ=255)                                                                                                                      cGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGcc                      <  8:223975/149‑1 (MQ=255)                                                                                                                      cGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGc                       <  3:182165/148‑1 (MQ=255)                                                                                                                      cGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGc                       <  8:154024/148‑1 (MQ=255)                                                                                                                       ggAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGCTTGACCATCTGCTCGTGGGTGGTTTCGCACCCCGACAGCAACAACGCAATGCACACGcc                      >  5:100525/1‑148 (MQ=255)                                                                                                                         aaGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGAcc                                                                           >  8:95673/1‑93 (MQ=255)                                                                                                                         aaGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGAcc                                                                           <  7:95673/93‑1 (MQ=255)                                                                                                                              gCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAg               <  8:366164/148‑1 (MQ=255)                                                                                                                               cccACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTa             >  5:141128/1‑149 (MQ=255)                                                                                                                                ccaccaTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTAc            <  6:141128/149‑1 (MQ=255)                                                                                                                                ccaccaTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTAc            >  5:176200/1‑149 (MQ=255)                                                                                                                                         gCCCGGACGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGCTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCCCAACCCGACAGCAACAACGCAATCCACACTCCGAACCAGAACCGAGACATa   >  2:66997/1‑149 (MQ=255)                                                                                                                                         cagtggCCGCCAGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACa     >  7:316980/5‑147 (MQ=255)                                                                                                                                          cccGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAg  >  4:13829/1‑149 (MQ=255)                                                                                                                                          agtggCCGCCAGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAAt                               <  7:162036/117‑1 (MQ=255)                                                                                                                                          agtggCCGCCAGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAAt                               >  8:162036/4‑120 (MQ=255)                                                                                                                                                    |                                                                                                                                           AAGCGCTCGTGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAG  >  NC_002947/790810‑791095

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