Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,882,092 A→G 6.7% Q218R (CAG→CGG) ‡ serA → bifunctional D‑3‑phosphoglycerate dehydrogenase/alpha‑ketoglutarate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,882,0920AG 6.7% 218.2 / 4.6 74Q218R (CAG→CGG) ‡serAbifunctional D‑3‑phosphoglycerate dehydrogenase/alpha‑ketoglutarate reductase
Reads supporting (aligned to +/- strand):  ref base A (37/32);  new base G (3/2);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

AGAACATGGGCATGCAGGTGTACTTCTTCGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGCGCCGCCATCGACGTGTTCC  >  NC_002947/5881951‑5882238
                                                                                                                                             |                                                                                                                                                  
agaACATGGGCATGCAGGTGTACTTCTTCGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCATCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGgat                                                                                                                                              <  1:306581/148‑1 (MQ=255)
   aCATGGGCATGCAGGTGTACTTCTTCGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACTGGGTGgatggca                                                                                                                                          <  7:298187/149‑4 (MQ=255)
     aTGGGCATGCAGGTGTACTTCTTCGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCgg                                                                                                                                        <  7:73746/149‑1 (MQ=255)
      tGGGCATGCAGGTGTACTTCTTCGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGc                                                                                                                                       <  3:348980/149‑1 (MQ=255)
       gggCATGCAGGTGTACTTCTTCGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGc                                                                                                                                       >  5:293231/1‑148 (MQ=255)
             gCAGGTGTACTTCTTCGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAg                                                                                                                                 >  1:14189/1‑148 (MQ=255)
               aGGTGTACTTCTTCGACCCGCTGACCAAGCTACCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGaa                                                                                                                              <  4:327353/149‑1 (MQ=255)
                      cttcttCGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCt                                                                                                                       >  3:93347/1‑149 (MQ=255)
                       ttcttcGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTc                                                                                                                      >  5:378341/1‑149 (MQ=255)
                       ttcttcGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTc                                                                                                                      >  3:194881/1‑149 (MQ=255)
                         cttcGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCg                                                                                                                                      <  5:188356/131‑1 (MQ=255)
                         cttcGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCg                                                                                                                                      >  6:188356/1‑131 (MQ=255)
                          ttcGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACTGGGTGgatggca                                                                                                                                          >  7:38411/1‑123 (MQ=255)
                          ttcGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACTGGGTGgatggca                                                                                                                                          <  8:38411/126‑4 (MQ=255)
                          ttcGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGt                                                                                                                                                  >  2:339329/1‑118 (MQ=255)
                          ttcGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGt                                                                                                                                                  <  1:339329/118‑1 (MQ=255)
                            cGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATg                                                                                                                  >  1:102849/1‑148 (MQ=255)
                              aCCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATg                                                                                                                  >  5:79040/1‑146 (MQ=255)
                                 cGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGaagaag                                                                                                            <  2:22515/149‑1 (MQ=255)
                                 cGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCCTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGaagaag                                                                                                            >  4:35073/1‑149 (MQ=255)
                                  gCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAgg                                                                                                           <  6:293231/149‑1 (MQ=255)
                                      aCCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGgcgc                                                                                                       >  1:214350/1‑149 (MQ=255)
                                      aCCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGgc                                                                                                         <  6:79040/147‑1 (MQ=255)
                                        cAAGCTGCAACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGAGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGa                                                                                                     >  5:226986/1‑149 (MQ=255)
                                               ccACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGa                                                                                               <  3:35073/148‑1 (MQ=255)
                                               ccACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGa                                                                                               <  2:14189/148‑1 (MQ=255)
                                               ccACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGAt                                                                                              <  6:99693/149‑1 (MQ=255)
                                                     ggCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGt                                                                                                                                                  >  2:376852/1‑91 (MQ=255)
                                                     ggCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGt                                                                                                                                                  <  1:376852/91‑1 (MQ=255)
                                                     ggCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAAcgc                                                                                        >  1:117504/1‑149 (MQ=255)
                                                      gCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAAcgcc                                                                                       >  6:75270/1‑149 (MQ=255)
                                                            ccGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGt                                                                                 >  8:153617/1‑149 (MQ=255)
                                                             cGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTg                                                                                >  3:83673/1‑149 (MQ=255)
                                                                       aCCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGaa                                                                                                                                  >  6:60109/1‑89 (MQ=255)
                                                                       aCCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGaa                                                                                                                                  <  5:60109/89‑1 (MQ=255)
                                                                         cAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTc                                                                     >  6:221272/1‑148 (MQ=255)
                                                                           gCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAg                                                                  <  6:226986/149‑1 (MQ=255)
                                                                             cTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCt                                                                >  6:170253/1‑149 (MQ=255)
                                                                              tGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTc                                                                     <  4:239576/143‑1 (MQ=255)
                                                                              tGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTc                                                                     >  3:239576/1‑143 (MQ=255)
                                                                               gAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATTCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTg                                                               >  1:5379/1‑148 (MQ=255)
                                                                                   gAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTgg                                                                               <  4:222345/128‑1 (MQ=255)
                                                                                   gAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTgg                                                                               >  3:222345/1‑128 (MQ=255)
                                                                                   gAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCGGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTg                                                                                <  7:245882/127‑1 (MQ=255)
                                                                                      ctgctgGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCt                                                       <  7:153617/149‑1 (MQ=255)
                                                                                        gctgGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCCAGCTCGaccacggtg                                                     <  8:152767/149‑5 (MQ=255)
                                                                                                     gACATCGTCTCGCTTCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCaa                                        <  7:206679/149‑1 (MQ=255)
                                                                                                        aTCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGa                                     <  3:82098/149‑1 (MQ=255)
                                                                                                        aTCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGACGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGa                                     >  1:176834/1‑149 (MQ=255)
                                                                                                        aTCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACACAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGa                                     >  6:359194/1‑149 (MQ=255)
                                                                                                               cGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCAc                              >  7:316278/1‑149 (MQ=255)
                                                                                                               cGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGCGCGCGATCCTGATCAACGACGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCAGCAAGGACAAGCAc                              >  3:37987/1‑149 (MQ=255)
                                                                                                                gCTGCACGTGCCTGAGCTGCCATGCCCCCAGTGGATGATCGGCGAGAATGAAATCCGCTCGATTTAGAAGGGCGCTATGCTTATCAACGCCTCCCGTGGCTCCGTGGTCGAGCTGTACCACCTGGCCGCTGCCATCAAGGACAAGCAc                              <  6:378341/148‑1 (MQ=255)
                                                                                                                gCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCAc                              <  4:194881/148‑1 (MQ=255)
                                                                                                                gCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCAc                              <  5:221272/148‑1 (MQ=255)
                                                                                                                gCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCAc                              <  2:214350/148‑1 (MQ=255)
                                                                                                                 cTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCt                            >  8:143154/1‑149 (MQ=255)
                                                                                                                     aCGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATc                        >  2:40062/1‑149 (MQ=255)
                                                                                                                     aCGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATc                        >  2:370401/1‑149 (MQ=255)
                                                                                                                          ccTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGgcg                    <  2:117504/148‑1 (MQ=255)
                                                                                                                          ccTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGgccc                   >  5:33242/1‑147 (MQ=255)
                                                                                                                           cTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGcgcc                  <  1:370401/149‑1 (MQ=255)
                                                                                                                            tGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGcgccg                 >  1:110154/1‑149 (MQ=255)
                                                                                                                             gAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAAACCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGACCGTGGCACCGTGTTCGAGCTGGACCACCTGGCCGCTGTCATCAAGGACAATCAACTCAACGGcgccgc                >  4:210545/1‑149 (MQ=255)
                                                                                                                              aGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGcgccgc                >  2:71761/1‑148 (MQ=255)
                                                                                                                               gCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACaa                                  <  6:207501/129‑1 (MQ=255)
                                                                                                                               gCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACaa                                  >  5:207501/1‑129 (MQ=255)
                                                                                                                               gCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGCGCCGCCa              >  3:68990/1‑149 (MQ=255)
                                                                                                                                cgatcATCCACTGGGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGCGCCGCCAt             >  8:155338/5‑149 (MQ=255)
                                                                                                                                 gatcATCCACTGGGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGCGCCGCCATc            >  3:333643/4‑149 (MQ=255)
                                                                                                                                   ccATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGCGCCGCCATc            <  5:170253/147‑1 (MQ=255)
                                                                                                                                       ccACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGCGCCGCCATCGACgtg      >  2:297011/1‑149 (MQ=255)
                                                                                                                                          cccAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGCGCCGCCATCGACGTGTTc   <  4:83673/149‑1 (MQ=255)
                                                                                                                                           ccAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGCGCCGCCATCGACGTGTTcc  <  7:143154/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
AGAACATGGGCATGCAGGTGTACTTCTTCGACCCGCTGACCAAGCTGCCACTGGGCAACGCCGTGCAGGTCACCAGCCTGAACGAGCTGCTGGGCCTGGCTGACATCGTCTCGCTGCACGTGCCTGAGCTGCCATCCACCCAGTGGATGATCGGCGAGAAGGAAATCCGCTCGATGAAGAAGGGCGCGATCCTGATCAACGCCGCCCGTGGCACCGTGGTCGAGCTGGACCACCTGGCCGCTGCCATCAAGGACAAGCACCTGATCGGCGCCGCCATCGACGTGTTCC  >  NC_002947/5881951‑5882238

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: