Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,750,399 A→G 7.6% H388H (CAT→CAC)  thiI ← tRNA sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,750,3990AG 7.6% 228.0 / 7.3 79H388H (CAT→CAC) thiItRNA sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base G (2/4);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

ACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTG  >  NC_002947/5750254‑5750546
                                                                                                                                                 |                                                                                                                                                   
aCGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGa                                                                                                                                                  <  4:18082/149‑1 (MQ=255)
    cATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGa                                                                                                                                              >  5:268400/1‑149 (MQ=255)
          cTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGaaa                                                                                                                                        >  4:248765/1‑149 (MQ=255)
             gACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGATgagctg                                                                                                                                       <  7:75315/147‑6 (MQ=255)
               cTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAAc                                                                                                                                   <  5:104063/149‑1 (MQ=255)
                ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACt                                                                                                                                  >  7:169507/1‑149 (MQ=255)
                   gATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTgg                                                                                                                               >  5:303756/1‑149 (MQ=255)
                   gATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTgg                                                                                                                               >  4:184403/1‑149 (MQ=255)
                        cTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGgcgcgc                                                                                                                          >  6:24416/1‑149 (MQ=255)
                         ttCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCt                                                                                                                         >  5:302178/1‑149 (MQ=255)
                            cAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCaa                                                                                                                      <  8:169507/149‑1 (MQ=255)
                            cAGCACCAGAGGTTCGGCCTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCaa                                                                                                                      <  4:252393/149‑1 (MQ=255)
                             aGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGAt                                                                                                                                             >  2:49195/1‑125 (MQ=255)
                             aGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGAt                                                                                                                                             <  1:49195/125‑1 (MQ=255)
                             aGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTc                                                                                                                        >  2:71781/1‑146 (MQ=255)
                              gCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCaa                                                                                                                      >  5:365662/1‑147 (MQ=255)
                                    gagGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGAt                                                                                                                                             <  1:92385/118‑1 (MQ=255)
                                    gagGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGAATTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGAt                                                                                                                                             >  2:92385/1‑118 (MQ=255)
                                     agGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgcc                                                                                                             <  3:184403/149‑2 (MQ=255)
                                          cGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAg                                                                                                         <  3:253580/148‑1 (MQ=255)
                                              ttCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACgg                                                                                                    >  5:78463/1‑149 (MQ=255)
                                                 cTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCa                                                                                                 >  3:26585/1‑149 (MQ=255)
                                                  tGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCATGGCCAt                                                                                                >  8:147455/1‑149 (MQ=255)
                                                   gggCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATg                                                                                               >  7:26484/1‑149 (MQ=255)
                                                    ggCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCAt                                                                                                <  7:160557/147‑1 (MQ=255)
                                                     gCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAaga                                                                                                                    <  3:204005/126‑1 (MQ=255)
                                                     gCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAaga                                                                                                                    >  4:204005/1‑126 (MQ=255)
                                                     gCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTc                                                                                             <  6:302178/149‑1 (MQ=255)
                                                        tCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATgg                                                                                                                                           >  5:28614/1‑100 (MQ=255)
                                                        tCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATgg                                                                                                                                           <  6:28614/100‑1 (MQ=255)
                                                        tCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGaa                                                                                          <  1:71781/149‑1 (MQ=255)
                                                           gggTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTg                                                                                       >  7:202637/1‑149 (MQ=255)
                                                                  gAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttct                                                                                <  4:26585/149‑1 (MQ=255)
                                                                  gAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttct                                                                                >  3:357608/1‑149 (MQ=255)
                                                                  gAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttct                                                                                <  6:78463/149‑1 (MQ=255)
                                                                    aTGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCATGGCCATGTCGAACTGCTTttcttc                                                                              <  7:147455/149‑1 (MQ=255)
                                                                    aTGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttcttc                                                                              <  3:66280/149‑1 (MQ=255)
                                                                     tGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCg                                                                             >  6:34368/1‑149 (MQ=255)
                                                                       tCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGt                                                                            <  6:365662/148‑1 (MQ=255)
                                                                           tGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTc                                                                        >  5:91513/1‑148 (MQ=255)
                                                                                 tCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCgg                                                                 >  8:264891/1‑149 (MQ=255)
                                                                                    ggcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGa                                                              >  7:48152/1‑149 (MQ=255)
                                                                                     gcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGc                                                                   <  5:297142/143‑1 (MQ=255)
                                                                                     gcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCgg                                                                 <  8:202637/145‑1 (MQ=255)
                                                                                     gcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCg                                                                  <  6:303756/144‑1 (MQ=255)
                                                                                          gcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTg                                                        >  6:275839/1‑149 (MQ=255)
                                                                                          gcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTg                                                        >  3:185109/1‑149 (MQ=255)
                                                                                           cAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTgg                                                       <  2:160088/149‑1 (MQ=255)
                                                                                            aGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGc                                                      <  7:264891/149‑1 (MQ=255)
                                                                                                ccTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATgg                                                  <  6:91513/149‑1 (MQ=255)
                                                                                                ccTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATgg                                                  >  5:169667/1‑149 (MQ=255)
                                                                                                ccTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATgg                                                  <  3:248765/149‑1 (MQ=255)
                                                                                                 cTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATggg                                                 <  6:268400/149‑1 (MQ=255)
                                                                                                  tCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTGTTCTTCGTGCTGCATGCGGTGACGCTTGGCATGggt                                                >  4:224269/1‑149 (MQ=255)
                                                                                                      ccACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTg                                                        <  5:329577/137‑1 (MQ=255)
                                                                                                      ccACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTg                                                        >  6:329577/1‑137 (MQ=255)
                                                                                                        aCTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCttt                                                                                   >  1:360716/1‑108 (MQ=255)
                                                                                                        aCTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCttt                                                                                   <  2:360716/108‑1 (MQ=255)
                                                                                                          ttCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGt                                        >  2:294256/1‑149 (MQ=255)
                                                                                                              tCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGggtg                                               >  8:285761/1‑138 (MQ=255)
                                                                                                              tCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGggtg                                               <  7:285761/138‑1 (MQ=255)
                                                                                                              tCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTAc                                    >  7:46639/1‑149 (MQ=255)
                                                                                                               cGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATggg                                                 <  8:184844/135‑1 (MQ=255)
                                                                                                               cGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATggg                                                 >  7:184844/1‑135 (MQ=255)
                                                                                                                gATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTAcc                                   <  5:275839/148‑1 (MQ=255)
                                                                                                                  tttCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCgag                                <  5:24416/149‑1 (MQ=255)
                                                                                                                    tCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCggg                                         >  3:348194/1‑138 (MQ=255)
                                                                                                                    tCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCggg                                         <  4:348194/138‑1 (MQ=255)
                                                                                                                      gATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCa                            >  6:380578/1‑149 (MQ=255)
                                                                                                                          tCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCAcccc                        >  7:161697/1‑149 (MQ=255)
                                                                                                                            cTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTg                                                        <  3:27337/115‑1 (MQ=255)
                                                                                                                            cTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTg                                                        >  4:27337/1‑115 (MQ=255)
                                                                                                                               gCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAg                    <  4:357608/148‑1 (MQ=255)
                                                                                                                               gCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAg                    <  2:131089/148‑1 (MQ=255)
                                                                                                                               gCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAg                    <  1:294256/148‑1 (MQ=255)
                                                                                                                                     tccATCGATCATGTGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCg             <  2:237438/147‑1 (MQ=255)
                                                                                                                                     cTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACCCCGCAGTATTCCg             <  8:161697/149‑1 (MQ=255)
                                                                                                                                           gATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAtgt       <  8:46639/149‑1 (MQ=255)
                                                                                                                                            aTCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAtgtg      >  8:18274/1‑149 (MQ=255)
                                                                                                                                               acaTGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCtt   <  4:185109/149‑1 (MQ=255)
                                                                                                                                                caTGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGATCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTg  <  3:224269/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTG  >  NC_002947/5750254‑5750546

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