Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,090,051	C→T	9.2%	G296S (GGC→AGC) ‡	astA‑II ←	arginine N‑succinyltransferase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,090,051	0	C	T	9.2%	180.1 / 6.8	65	G296S (GGC→AGC) ‡	astA‑II	arginine N‑succinyltransferase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base T (3/3);  total (32/33)
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

TTGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTGCTCGGTC  >  NC_002947/5089907‑5090199                                                                                                                                                  |                                                                                                                                                     ttGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCg                                                                                                                                                     <  2:238846/146‑1 (MQ=255) ttGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTg                                                                                                                                                  >  1:23778/1‑149 (MQ=255)    ccAGACCTCAAACCGCGACCCGGCGCACGCTGGCGCCCTCCCCCACGCCCAGCGAGTCGGCTGCTGCCACGCTCAGGCTGCCCGGCTTGCCAGGCGCCCAGTCCACGTCCAGCA‑ACCCAGCGCGGTAGTCCTGCAACAGGCCGGTGca                                                                                                                                                >  4:323386/1‑148 (MQ=255)    ccAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGcac                                                                                                                                               >  4:85641/1‑149 (MQ=255)     cAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACGGCCTGTTGcagg                                                                                                                                              <  1:322904/149‑3 (MQ=255)      aGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACGGCCTGTTGcag                                                                                                                                               >  6:109001/1‑146 (MQ=255)      aGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACGGCCTGTTGcag                                                                                                                                               <  5:109001/147‑2 (MQ=255)      aGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCacac                                                                                                                                             >  2:230233/1‑149 (MQ=255)       gACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGTTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACAcc                                                                                                                                            <  4:259762/149‑1 (MQ=255)       gACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACAcc                                                                                                                                            <  1:230233/149‑1 (MQ=255)           tCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACGGCCTGTTGcag                                                                                                                                               >  8:9574/1‑141 (MQ=255)           tCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACGGCCTGTTGcag                                                                                                                                               <  7:9574/142‑2 (MQ=255)           tCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAgg                                                                                                                                         >  1:82643/1‑148 (MQ=255)            cAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTa                                                                                                                                       >  8:99519/1‑149 (MQ=255)              aaCCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACg                                                                                                                                     >  4:128491/1‑149 (MQ=255)               aCCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACgg                                                                                                                                    <  5:197944/149‑1 (MQ=255)                ccGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACggg                                                                                                                                   >  8:6311/1‑149 (MQ=255)                                gCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCg                                                                                                                   >  3:303742/1‑149 (MQ=255)                                         cTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCTTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCtt                                                                                                          <  1:241094/149‑1 (MQ=255)                                             cccACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGAt                                                                                                      <  5:2309/149‑1 (MQ=255)                                              ccACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATc                                                                                                     >  8:166901/1‑149 (MQ=255)                                              ccACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATc                                                                                                     >  1:317254/1‑149 (MQ=255)                                                   cccAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCAc                                                                                                <  3:134823/149‑1 (MQ=255)                                                   cccAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCAGCCGTTGACCGGGGCCTCTTCGATCTTCAc                                                                                                >  7:252139/1‑149 (MQ=255)                                                     cAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCg                                                                                              <  3:273093/149‑1 (MQ=255)                                                         gcgTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGcac                                                                                          <  3:85641/149‑1 (MQ=255)                                                           gTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGcacca                                                                                        <  7:99519/149‑1 (MQ=255)                                                                    tgcCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGctct                                                                               >  4:161063/1‑149 (MQ=255)                                                                      ccACGCTCAGTCTGCCCGGCTTGCCAGTCGCGCATTCGAGGTCCAGCA‑ACAGCGCGGGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACGCGGCTCTg                                                                              <  6:233019/148‑1 (MQ=255)                                                                       cACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTggg                                                                            <  3:128491/149‑1 (MQ=255)                                                                         cGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCg                                                                          <  4:303742/149‑1 (MQ=255)                                                                           cTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGAt                                                                        >  5:43332/1‑149 (MQ=255)                                                                              aGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATc                                                                       >  1:44050/1‑147 (MQ=255)                                                                              aGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATc                                                                       <  6:351345/147‑1 (MQ=255)                                                                                gCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGc                                                                   <  2:82643/149‑1 (MQ=255)                                                                                      cGGCTTGCCAGGCACCCAGTCCACGTCCAGCA‑ACACCGCGCGGTAGTCCTCCAACACGCCGTTGGACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCGCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGc                                                             <  3:323386/149‑1 (MQ=255)                                                                                        gCTTGCCAGGCACCCAGTCCAGGTCCAGCAGACACCCCGCGGTATTCATGCAACAGGCCGTTGCACACCAGTTACGGCCGCCCGCCTTTGCCCGGGGCCTCTTCGATATTCCCGGGCCCCACCCGGATCTGGGCGAGGGACCGCATGc                                                             >  1:205741/1‑148 (MQ=255)                                                                                           tGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCaga                                                         <  3:179051/148‑1 (MQ=255)                                                                                           tGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCa                                                           >  4:179051/1‑146 (MQ=255)                                                                                            gCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCagag                                                        <  5:180432/148‑1 (MQ=255)                                                                                              cAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCttt                                                                                                                        >  4:72917/1‑82 (MQ=255)                                                                                              cAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCttt                                                                                                                        <  3:72917/82‑1 (MQ=255)                                                                                                 gCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTgcgc                                                  >  3:72020/1‑149 (MQ=255)                                                                                                 gCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTgcgc                                                  <  7:166901/149‑1 (MQ=255)                                                                                                        gtccaggtccagCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATg                                             >  2:317917/1‑147 (MQ=255)                                                                                                           caggtccagCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAgg                                         <  2:23778/148‑1 (MQ=255)                                                                                                                 cagCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCagag                                                        >  4:44861/1‑127 (MQ=255)                                                                                                                 cagCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCagag                                                        <  3:44861/127‑1 (MQ=255)                                                                                                                 cagCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGgc                                  <  2:317254/149‑1 (MQ=255)                                                                                                                 cagCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGgc                                  >  4:209827/1‑149 (MQ=255)                                                                                                                    cA‑ACACCGCGCGGTAGACCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTACGATCTTCACCGCCACCACCCGGCTCTGCGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTc                                      >  6:186137/1‑142 (MQ=255)                                                                                                                     cctccgcgcgcgGTAGTCCTGCAACCGGCCGTTGCACACCAGGTACGGGCGCCGGCGTGTGAGCGGGGCCTCTTCGATCTTCACCGGCACCACCGGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGGATGCAGGGTc                                      <  5:186137/136‑1 (MQ=255)                                                                                                                           cgcgcgGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCg                          >  5:206091/1‑148 (MQ=255)                                                                                                                                ggTAGTCCTGCACCTGGCCGTTGCACCCCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGTCCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATg                    <  8:252139/149‑1 (MQ=255)                                                                                                                                      tgtGCAACGGCCTGTTACACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATg              >  2:163611/3‑149 (MQ=255)                                                                                                                                       cTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGa                                                                     >  4:127151/1‑93 (MQ=255)                                                                                                                                       cTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGa                                                                     <  3:127151/93‑1 (MQ=255)                                                                                                                                       cTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATg              >  5:74345/1‑148 (MQ=255)                                                                                                                                        tGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATg              >  1:60980/1‑147 (MQ=255)                                                                                                                                        tGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATg              >  5:131700/1‑147 (MQ=255)                                                                                                                                         gCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAg           >  4:114266/1‑149 (MQ=255)                                                                                                                                          cAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTa            <  1:317917/147‑1 (MQ=255)                                                                                                                                          cAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGACGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAg           <  1:276642/148‑1 (MQ=255)                                                                                                                                              aGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTGCTc      <  7:6311/149‑1 (MQ=255)                                                                                                                                               ggCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTGCTCg     <  4:72020/149‑1 (MQ=255)                                                                                                                                                  cGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTGCTCGGtc  >  8:76183/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TTGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCA‑ACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTGCTCGGTC  >  NC_002947/5089907‑5090199

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