Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,216,203	A→G	6.5%	F168S (TTC→TCC)	PP_4596 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,216,203	0	A	G	6.5%	175.4 / 2.4	61	F168S (TTC→TCC)	PP_4596	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base G (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

ATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGACATAGTGGTCGGTGGCGGC  >  NC_002947/5216056‑5216351                                                                                                                                                    |                                                                                                                                                      aTGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGaa                                                                                                                                                      <  7:139430/149‑1 (MQ=255)    aCGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACa                                                                                                                                                    <  6:147681/148‑1 (MQ=255)         tatCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGt                                                                                                                                              <  1:363205/149‑1 (MQ=255)             gcgTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTaaa                                                                                                                                           <  5:156740/148‑1 (MQ=255)             gcgTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAg                                                                                                                                          <  4:241884/149‑1 (MQ=255)                     cAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAg                                                                                                                                          <  8:44382/141‑1 (MQ=255)                     cAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAg                                                                                                                                          >  7:44382/1‑141 (MQ=255)                            aGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGtt                                                                                                                          <  1:394210/149‑1 (MQ=255)                               tCGTGGAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCt                                                                                                                       <  4:353249/149‑1 (MQ=255)                                  tGGAACACCACGTTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTGGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTggc                                                                                                                    <  8:374128/149‑1 (MQ=255)                                    gAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACTCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTggccc                                                                                                                  >  2:170061/1‑149 (MQ=255)                                      acacCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTggcccgg                                                                                                                <  5:271670/149‑1 (MQ=255)                                         ccacGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTggcccggccc                                                                                                             >  5:271727/1‑149 (MQ=255)                                          cacGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCAGCT‑TCGTTGCTggcccggccc                                                                                                             >  8:173539/1‑148 (MQ=255)                                            cGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAg                                                                                                          >  1:244938/1‑149 (MQ=255)                                              gTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGgt                                                                                                        >  4:282687/1‑149 (MQ=255)                                                gCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTCGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGgtg                                                                                                       >  5:332603/1‑148 (MQ=255)                                                  ccTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGgtg                                                                                                       >  6:307268/1‑146 (MQ=255)                                                  ccTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAg                                                                                                    <  1:190944/149‑1 (MQ=255)                                                   cTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAg                                                                                                    >  8:299350/1‑148 (MQ=255)                                                        ccACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCg                                                                                              <  5:326177/149‑1 (MQ=255)                                                         cACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCgg                                                                                             <  1:170061/149‑1 (MQ=255)                                                         cACCCCGAGGGTTACCGACATGCCCGCGCCCCTGGCTTCGGGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGCAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCgg                                                                                             <  3:294928/149‑1 (MQ=255)                                                           ccccGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCt                                                                                           >  8:131519/1‑149 (MQ=255)                                                           ccccGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCt                                                                                           <  6:332603/149‑1 (MQ=255)                                                                        cGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACt                                                                              >  6:353412/1‑149 (MQ=255)                                                                          aCATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTggcc                                                                                                                   <  6:17298/110‑1 (MQ=255)                                                                          aCATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTggcc                                                                                                                   >  5:17298/1‑110 (MQ=255)                                                                             tGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAgc                                                                         >  7:187782/1‑149 (MQ=255)                                                                              gCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGGACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAgcg                                                                        >  3:169565/1‑149 (MQ=255)                                                                              gCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGGACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAgcg                                                                        >  3:175503/1‑149 (MQ=255)                                                                              gCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCcgcg                                                                        >  2:123142/1‑149 (MQ=255)                                                                               cccGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAgcgc                                                                       <  7:299350/149‑1 (MQ=255)                                                                                ccGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTGGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGcgca                                                                      >  6:95221/1‑149 (MQ=255)                                                                                ccGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGcgca                                                                      >  4:182810/1‑149 (MQ=255)                                                                                   cgccgcTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGcgcacgc                                                                   <  5:307268/149‑1 (MQ=255)                                                                                       gcTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGt                                                               >  4:235318/1‑149 (MQ=255)                                                                                           gCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTGGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCt                                                           <  5:95221/149‑1 (MQ=255)                                                                                           gCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGc                                                             <  3:182810/147‑1 (MQ=255)                                                                                            cTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTg                                                          <  3:282687/149‑1 (MQ=255)                                                                                             ttCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAGGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTgg                                                         >  8:46750/1‑149 (MQ=255)                                                                                                 tGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGaa                                                                                                                                                      >  8:140974/1‑53 (MQ=255)                                                                                                 tGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGaa                                                                                                                                                      <  7:140974/53‑1 (MQ=255)                                                                                                   gCACATGTTCGAGCAATCCGGGAAGTTGTTGCTGCCGTAGGCACGCTCGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGCTGAAGAATTCGGCTTCATCCGGCGACTGAAGCGCACGCAGGTGCCTGGCGAc                                                     >  5:290571/1‑147 (MQ=255)                                                                                                     acaTGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGc                                                 >  5:202036/1‑149 (MQ=255)                                                                                                      caTGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTTGCGACCATCt                                                >  2:270071/1‑149 (MQ=255)                                                                                                            cGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAg                                           <  3:235318/148‑1 (MQ=255)                                                                                                             gAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGc                                         >  4:275294/1‑149 (MQ=255)                                                                                                                     cGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGc                                  >  6:225009/1‑148 (MQ=255)                                                                                                                          aGTTGTTGGTGCCGTAGGCACGCACTAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAggtg                            <  7:173539/149‑1 (MQ=255)                                                                                                                          aGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAggtg                            <  5:136130/149‑1 (MQ=255)                                                                                                                           gttgttGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAggtgg                           >  1:302868/1‑149 (MQ=255)                                                                                                                            ttgttgGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCTGTCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAggtgg                           >  1:92273/1‑149 (MQ=255)                                                                                                                               tttgTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGGCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATGCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTc                       <  6:290571/146‑1 (MQ=255)                                                                                                                                 ggTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGa                     >  6:248426/1‑149 (MQ=255)                                                                                                                                  gTGCCGTAGGCACGCACGGACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGAc                    <  4:169565/149‑1 (MQ=255)                                                                                                                                  gTGCCGTAGGCACGCACGGACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGAc                    <  4:175503/149‑1 (MQ=255)                                                                                                                                    gCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGACAt                  <  5:248426/149‑1 (MQ=255)                                                                                                                                     cctacggCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGACATa                 <  1:293186/144‑1 (MQ=255)                                                                                                                                     ccGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGACa                   >  2:211895/1‑147 (MQ=255)                                                                                                                                        tAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGACATAGTg              <  2:328702/149‑1 (MQ=255)                                                                                                                                          ggCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGACATAGTGGt            <  2:92273/149‑1 (MQ=255)                                                                                                                                               gcacgAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGACATAGTGGTCGGt        <  5:323247/148‑1 (MQ=255)                                                                                                                                               gcacgAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGACATAGTGGTCGGt        >  6:323247/1‑148 (MQ=255)                                                                                                                                                    aaCAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGACATAGTGGTCGGTggcggc  <  5:225009/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      ATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGAGGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCT‑TCGTTGCTGGCCCGGCCCGAGGTGTAGAATTCGGCTTCATCCGGCGACTGCAGCGCACGCAGGTGCCTGGCGACCAGCTCGAAGGCCTCCTGCCAGGTGGTTTCGACATAGTGGTCGGTGGCGGC  >  NC_002947/5216056‑5216351

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