Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,158,943	T→G	9.5%	L321R (CTC→CGC)	opuA →	glycine betaine ABC transporter ATP‑binding protein/permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,158,943	0	T	G	9.5%	170.9 / 5.4	60	L321R (CTC→CGC)	opuA	glycine betaine ABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base T (31/23);  new base G (2/4);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGGTTTTCTTCGACGAGCGGCCGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGT  >  NC_002947/3158796‑3159091                                                                                                                                                    |                                                                                                                                                      tGGCGGTTTTCTTCGACGAGCGGCCGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCGAAAGCGTGTTGGGCCCGGCAGcgctc                                                                                                                                                      >  8:188983/1‑149 (MQ=255)     ggTTTTCTTCGACGAGCGGCCGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTg                                                                                                                                                  >  2:317386/1‑149 (MQ=255)           cttcGACGAGCGGCCGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATGGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGAGCGCTgccggg                                                                                                                                             <  7:67844/148‑3 (MQ=255)            ttcGACGAGCGGCCGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGccgccg                                                                                                                                            >  7:287481/1‑148 (MQ=255)               gacgaGCGGCCGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑Ac                                                                                                                                        <  5:176246/148‑1 (MQ=255)               gacgaGCGGCCGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑Ac                                                                                                                                        >  1:67956/1‑148 (MQ=255)               gacgaGCGGCCGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑Ac                                                                                                                                        >  6:176246/1‑148 (MQ=255)                  gaGCGGCCGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAgg                                                                                                                                    >  1:31781/1‑149 (MQ=255)                       gCCGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGtggtt                                                                                                                               <  3:59652/149‑1 (MQ=255)                        ccGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGCGCTGCCAGCGCGCTCGACAATCCATTGGAGTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAggtggtgg                                                                                                                              >  1:34427/1‑147 (MQ=255)                        ccGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGtggttg                                                                                                                              <  1:317386/149‑1 (MQ=255)                        ccGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGtggttg                                                                                                                              >  1:90048/1‑149 (MQ=255)                        ccGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGtggttg                                                                                                                              <  2:31781/149‑1 (MQ=255)                            atatCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGtggtt                                                                                                                               <  3:121101/144‑1 (MQ=255)                            atatCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGtggtt                                                                                                                               >  4:121101/1‑144 (MQ=255)                            atatCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTc                                                                                                                          >  7:33744/1‑149 (MQ=255)                               tctgCCCGAGCATTGCGGCGTCTAGCGACGAGGTGCTCTGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGGGGGGGGCCCGGCAGCGCGCGCTGCCGCCGG‑ACGAGGTGGTTGGTGct                                                                                                                         <  2:34427/144‑1 (MQ=255)                               tCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCTGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTctgg                                                                                                                       <  2:67956/149‑1 (MQ=255)                                cgtcccGAGCATTGCGTCGTGTTCCCACGCGGTGCTCAGCGAGCCTGCTGCGGCCGTGCTGCCTGCGCGTTCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCTGCGCTGGCTGCGGCCGG‑ACGAGGTGGTTGGTTctgg                                                                                                                       <  8:186343/145‑1 (MQ=255)                                   cccGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTctggctgg                                                                                                                   <  7:26027/149‑1 (MQ=255)                                    ccGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGtggttg                                                                                                                              <  3:159623/137‑1 (MQ=255)                                    ccGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGtggttg                                                                                                                              >  4:159623/1‑137 (MQ=255)                                     cGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGgcg                                                                                                                 >  2:226833/1‑149 (MQ=255)                                            tGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGctgctg                                                                                                           <  2:90048/148‑1 (MQ=255)                                             gCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGcgct                                                                                                                                                       >  7:227781/1‑104 (MQ=255)                                             gCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGcgct                                                                                                                                                       <  8:227781/104‑1 (MQ=255)                                                     tAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTgg                                                                                                 >  5:90945/1‑149 (MQ=255)                                                       gCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTggc                                                                                                <  6:90945/148‑1 (MQ=255)                                                                     aGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCAtgc                                                                                  <  8:287481/148‑1 (MQ=255)                                                                         aGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCAtgcttgct                                                                             >  4:116258/1‑149 (MQ=255)                                                                                    gCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCtggcgcgtggcg                                                                  >  3:6555/1‑149 (MQ=255)                                                                                         gCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAgg                                                             >  7:315105/1‑149 (MQ=255)                                                                                          cTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAggg                                                            <  3:116258/149‑1 (MQ=255)                                                                                             ccAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAgggg                                                           >  4:131848/1‑147 (MQ=255)                                                                                             ccAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGc                                                         >  3:150709/1‑149 (MQ=255)                                                                                              cAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCg                                                        >  5:225378/1‑149 (MQ=255)                                                                                                    gATCGACATTCTATTGGATTCGCCAATGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTgg                                                  <  4:150709/149‑1 (MQ=255)                                                                                                    gATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGTACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTg                                                   >  1:292015/1‑149 (MQ=255)                                                                                                      tCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCgctg                                                                                                              >  6:168671/1‑87 (MQ=255)                                                                                                      tCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCgctg                                                                                                              <  5:168671/87‑1 (MQ=255)                                                                                                        gACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCtggcgcgtg                                                                     <  7:88452/126‑1 (MQ=255)                                                                                                        gACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCtggcgcgtg                                                                     >  8:88452/1‑126 (MQ=255)                                                                                                         aCAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATc                                             >  7:89999/1‑149 (MQ=255)                                                                                                         aCAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTATGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATc                                             >  5:340666/1‑149 (MQ=255)                                                                                                            aTCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGc                                          <  8:33744/149‑1 (MQ=255)                                                                                                                aTTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCca                                        >  7:206992/1‑147 (MQ=255)                                                                                                                aTTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCca                                        >  1:335418/1‑147 (MQ=255)                                                                                                                aTTGGATTCGCCAAAGCGTGTTGGGCCAGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGAATGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCca                                        >  1:300007/1‑147 (MQ=255)                                                                                                                    gATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCtggcgcgtgg                                                                    >  1:134451/1‑115 (MQ=255)                                                                                                                    gATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCtggcgcgtgg                                                                    <  2:134451/115‑1 (MQ=255)                                                                                                                    gATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCAc                                  >  1:30892/1‑149 (MQ=255)                                                                                                                        cGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTgg                               <  8:315105/148‑1 (MQ=255)                                                                                                                         gCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTgg                               >  4:326327/1‑147 (MQ=255)                                                                                                                                cGTGTTGGGCCCGGCAGCGCGCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGtgctgc                      <  7:188983/149‑1 (MQ=255)                                                                                                                                   gtTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTgg                   >  2:27258/1‑149 (MQ=255)                                                                                                                                    ttGGGCCCGGCAGCGCTCGCTGCTGCCGG‑ACGTGGTGGTTGGTTGTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGt                  <  6:340666/149‑1 (MQ=255)                                                                                                                                    ttGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTgg                   <  1:226833/148‑1 (MQ=255)                                                                                                                                     tGGGCGCGGCTGCGCTCGCTGCGGCCGG‑ACGTGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGGGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGt                  <  2:335418/148‑1 (MQ=255)                                                                                                                                         ggcgtcAGCGAGCGCTGCCGCCGG‑ACGAGGTGGTTGTTTCTGGCTGGCGGTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGtgct                        <  4:116731/133‑1 (MQ=255)                                                                                                                                         ggcGGCAGCGAGCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGtgct                        >  3:116731/3‑138 (MQ=255)                                                                                                                                          gcGGCAGCGAGCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCg            >  2:231686/2‑149 (MQ=255)                                                                                                                                              cAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTg        <  4:6555/149‑1 (MQ=255)                                                                                                                                                gcgcTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTggg      >  1:144821/1‑149 (MQ=255)                                                                                                                                                    tcgctGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTtgt  >  4:47454/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TGGCGGTTTTCTTCGACGAGCGGCCGCATATCGGCCCGAGCATTGCGTCGTCTAGCGACGAGGTGCTCAGCGAGCCAGCGTCGGCCGAGCTGCCAGCGCGATCGACAATCCATTGGATTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGG‑ACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGT  >  NC_002947/3158796‑3159091

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