Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,356,515	T→G	6.0%	E749A (GAG→GCG)	infB ←	translation initiation factor IF‑2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,356,515	0	T	G	6.0%	208.6 / 4.2	67	E749A (GAG→GCG)	infB	translation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base T (26/36);  new base G (2/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CCTTCGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGAAGCCGAA  >  NC_002947/5356371‑5356660                                                                                                                                                  |                                                                                                                                                  ccTTCGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTcgcg                                                                                                                                               <  7:36654/149‑1 (MQ=255)  cTTCGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTcgcg                                                                                                                                               >  1:216355/1‑148 (MQ=255)      gAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACAt                                                                                                                                          <  7:218195/149‑1 (MQ=255)       aaGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATc                                                                                                                                         >  2:253538/1‑149 (MQ=255)         gATTACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCg                                                                                                                                        >  8:241383/1‑148 (MQ=255)           tCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTg                                                                                                                                     >  8:292011/1‑149 (MQ=255)             acaACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGcc                                                                                                                                   >  4:196457/1‑149 (MQ=255)             acaACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGc                                                                                                                                    <  7:247771/148‑1 (MQ=255)                  gTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACTCCGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGc                                                                                                                               <  1:96394/149‑1 (MQ=255)                      tcaCGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGc                                                                                                                           <  7:154332/148‑1 (MQ=255)                      tcaCGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGc                                                                                                                           >  8:154332/1‑148 (MQ=255)                       caCGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCg                                                                                                                         >  6:11342/1‑149 (MQ=255)                             gTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAgg                                                                                                                   <  7:241383/149‑1 (MQ=255)                                cgcgAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGc                                                                                                                 <  7:226759/148‑1 (MQ=255)                                   gAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCcttc                                                                                                             >  1:112017/1‑149 (MQ=255)                                     aTCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCAAGCCGGTCAGGGCcttctt                                                                                                           <  6:346269/149‑1 (MQ=255)                                        ggACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGCATGTACACGCTAACAACCCTCCCGAGCATGCCGGAATGGGCCTTCTTGAc                                                                                                        >  3:154738/1‑149 (MQ=255)                                        ggAAGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCATCGATGGAGCCGAACTTCGGCGAAATGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGAc                                                                                                        <  7:130020/149‑1 (MQ=255)                                         gACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGa                                                                                                         <  6:66020/147‑1 (MQ=255)                                         gACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTACTTGa                                                                                                         >  5:66020/1‑147 (MQ=255)                                             gTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCtt                                                                                                   <  3:81464/149‑1 (MQ=255)                                              ttGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTc                                                                                                  <  2:216355/149‑1 (MQ=255)                                                gCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCg                                                                                                 <  5:11342/148‑1 (MQ=255)                                                 cGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGTATGCCGGTCAGGGCCTTCTTGACGTCTTCgat                                                                                               <  7:292011/149‑1 (MQ=255)                                                  ggTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCgat                                                                                               >  1:308714/1‑148 (MQ=255)                                                       acGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTc                                                                                         <  4:256560/149‑1 (MQ=255)                                                         ggTATCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGACGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCg                                                                                        >  4:271539/1‑148 (MQ=255)                                                         ggTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCg                                                                                        >  4:318460/1‑148 (MQ=255)                                                          gTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTa                                                                                      <  2:112017/149‑1 (MQ=255)                                                                tCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCgat                                                                                               <  4:122403/134‑1 (MQ=255)                                                                tCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCgat                                                                                               >  3:122403/1‑134 (MQ=255)                                                                   aTGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGtt                                                                             <  4:127378/149‑1 (MQ=255)                                                                   aTGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGGCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGtt                                                                             <  1:253538/149‑1 (MQ=255)                                                                        cATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCATGTTGtaat                                                                        <  8:330687/149‑1 (MQ=255)                                                                            cagccagcGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAAc                                                                     <  3:196457/148‑1 (MQ=255)                                                                              gccagcGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCa                                                                  <  4:64924/149‑1 (MQ=255)                                                                               ccagcGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCa                                                                  >  4:147358/1‑148 (MQ=255)                                                                               ccagcGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCa                                                                  >  2:190803/1‑148 (MQ=255)                                                                                cagcGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCata                                                                <  3:318460/149‑1 (MQ=255)                                                                                cagcGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCata                                                                >  3:232425/1‑149 (MQ=255)                                                                                             gAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTTGATCACGTTGTAATAACGCATATCCAGACCTTCCt                                                   <  1:240624/149‑1 (MQ=255)                                                                                              aaCTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTACTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTg                                                  >  1:91356/1‑149 (MQ=255)                                                                                                      cGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTcgacga                                          >  6:374187/1‑149 (MQ=255)                                                                                                       gAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCCACCAAGCCGGTCAGCGCCTTCTAGACGTCTACGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTcgacga                                          >  1:70132/1‑148 (MQ=255)                                                                                                            gAAC‑ACGTCACGTACTTCGGCGTCACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTgc                                    <  2:351707/149‑1 (MQ=255)                                                                                                             aaC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCtt                                      <  1:83590/146‑1 (MQ=255)                                                                                                             aaC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCtt                                      >  2:83590/1‑146 (MQ=255)                                                                                                                    tcacgtACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCAcc                             >  2:54692/1‑149 (MQ=255)                                                                                                                      acgtACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCg                            <  1:190803/148‑1 (MQ=255)                                                                                                                      acgtACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCg                            <  4:98112/148‑1 (MQ=255)                                                                                                                           cTTCGGCGACTCCAAGGATGTTCTCGCGAACATGGCTGCCGAGCATGGCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTTTAATAACGCATATCCAGACCTTGCTGCTCGACGATCTTGCGCGCACCGGCATc                       <  2:70132/148‑1 (MQ=255)                                                                                                                            ttCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCgg                     >  2:184347/1‑149 (MQ=255)                                                                                                                              cGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCgtc                    <  3:271539/148‑3 (MQ=255)                                                                                                                                  gACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGgcacgca               >  5:227444/1‑149 (MQ=255)                                                                                                                                   acacCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGgcacgcac              <  6:227444/149‑1 (MQ=255)                                                                                                                                    cacCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGcacgcacg             >  1:326974/1‑149 (MQ=255)                                                                                                                                     acCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGt            <  1:54692/149‑1 (MQ=255)                                                                                                                                       ccAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTc                                                                                         >  6:272937/1‑70 (MQ=255)                                                                                                                                       ccAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTc                                                                                         <  5:272937/70‑1 (MQ=255)                                                                                                                                       ccAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGtt           >  8:287661/1‑148 (MQ=255)                                                                                                                                        agcgATGTTCGCGAGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTg          >  2:300752/4‑148 (MQ=255)                                                                                                                                         gcgATGTTCGCGAGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTg          >  4:167506/3‑147 (MQ=255)                                                                                                                                         gcgATGTTCGCGAGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGaa        <  3:167506/147‑1 (MQ=255)                                                                                                                                         aGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGaa        <  7:287661/149‑1 (MQ=255)                                                                                                                                         aGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGaa        <  3:64846/149‑1 (MQ=255)                                                                                                                                          cgATGTTCGCGAGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGACTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGAAg       <  1:300752/148‑1 (MQ=255)                                                                                                                                               ttCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGAAGCCGaa  <  2:308714/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CCTTCGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAAC‑ACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGAAGCCGAA  >  NC_002947/5356371‑5356660

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