Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 297,352 C→A 5.1% E218* (GAA→TAA)  gshA ← glutamate‑‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947297,3520CA 5.1% 254.2 / 3.2 78E218* (GAA→TAA) gshAglutamate‑‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (33/41);  new base A (2/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

CCCACGGCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTCGCGGTCATTTTCCGTACT  >  NC_002947/297204‑297499
                                                                                                                                                    |                                                                                                                                                    
cccACGGCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGTGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTc                                                                                                                                                      >  1:393329/1‑149 (MQ=255)
cccACGGCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTc                                                                                                                                                      <  5:259259/149‑1 (MQ=255)
cccACGGCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTc                                                                                                                                                      <  6:89916/149‑1 (MQ=255)
 ccACGGCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCg                                                                                                                                                     >  7:348147/1‑149 (MQ=255)
  cACGGCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTc                                                                                                                                                      >  1:386454/1‑147 (MQ=255)
  cACGGCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTc                                                                                                                                                      >  4:261922/1‑147 (MQ=255)
  cACGGCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTc                                                                                                                                                      >  5:171403/1‑147 (MQ=255)
   aCGGCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAg                                                                                                                                                   >  4:190800/1‑149 (MQ=255)
    cGGCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGc                                                                                                                                                  <  7:284571/149‑1 (MQ=255)
      gCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTg                                                                                                                                                <  7:50567/149‑1 (MQ=255)
       ccTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTgg                                                                                                                                               >  7:292964/1‑149 (MQ=255)
        cTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGt                                                                                                                                              <  3:108719/149‑1 (MQ=255)
        cTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCGAAGAGCTgga                                                                                                                                              <  4:160940/149‑2 (MQ=255)
           gcgcAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTgcg                                                                                                                                           >  2:67135/1‑149 (MQ=255)
              cAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTgcgggc                                                                                                                                        <  2:386454/149‑1 (MQ=255)
              cAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTgcgggc                                                                                                                                        <  2:414369/149‑1 (MQ=255)
              cAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGCTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTgcgggc                                                                                                                                        <  3:292354/149‑1 (MQ=255)
               aGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTgcgggcg                                                                                                                                       <  2:218032/149‑1 (MQ=255)
                        gtgtAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAg                                                                                                                              >  8:33033/1‑149 (MQ=255)
                         tgtAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAgg                                                                                                                             <  1:389084/149‑1 (MQ=255)
                              cTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGcc                                                                                                                        >  5:36016/1‑149 (MQ=255)
                               tGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCt                                                                                                                       >  7:29147/1‑149 (MQ=255)
                                ggCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGTTGCGGGCGGCCACGCAGGAAGCCtt                                                                                                                      >  2:195474/1‑149 (MQ=255)
                                 gCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCt                                                                                                                       <  8:277060/147‑1 (MQ=255)
                                   cAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGt                                                                                                                   >  2:54120/1‑149 (MQ=255)
                                      gttgttGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTcc                                                                                                                 <  8:348147/148‑1 (MQ=255)
                                      gttgttGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCa                                                                                                                <  2:393329/149‑1 (MQ=255)
                                        tgttgtAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAgg                                                                                                             >  6:124468/1‑149 (MQ=255)
                                         gttgtAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAgg                                                                                                             <  5:340346/148‑1 (MQ=255)
                                          ttgtAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCa                                                                                                                               <  5:107957/130‑1 (MQ=255)
                                          ttgtAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCa                                                                                                                               >  6:107957/1‑130 (MQ=255)
                                          ttgtAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGc                                                                                                           >  3:377158/1‑149 (MQ=255)
                                              aGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGc                                                                                                       >  1:67587/1‑149 (MQ=255)
                                              aGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGc                                                                                                       >  5:42387/1‑149 (MQ=255)
                                                cAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGa                                                                                                     <  6:171403/149‑1 (MQ=255)
                                                cAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGa                                                                                                     <  5:124468/149‑1 (MQ=255)
                                                   ggCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGgc                                                                                                  <  6:352868/149‑1 (MQ=255)
                                                     cGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTgcgggc                                                                                                                                        >  8:346788/1‑110 (MQ=255)
                                                     cGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTgcgggc                                                                                                                                        <  7:346788/110‑1 (MQ=255)
                                                       tAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCg                                                                                              <  7:202497/149‑1 (MQ=255)
                                                            ccggccTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGTTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAg                                                                                         <  1:195474/149‑1 (MQ=255)
                                                            ccggccTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAg                                                                                         <  8:292964/149‑1 (MQ=255)
                                                              ggccTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGt                                                                                       <  6:42387/149‑1 (MQ=255)
                                                              ggccTGGGCGTTGCGCTGGTACCCCAGGTCGCTCCTGCGCAGGCTGGGGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGt                                                                                       <  5:117131/149‑1 (MQ=255)
                                                                 cTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGa                                                                                                     >  4:163430/1‑132 (MQ=255)
                                                                 cTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGa                                                                                                     <  3:163430/132‑1 (MQ=255)
                                                                   gggCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACa                                                                                          >  4:266649/1‑141 (MQ=255)
                                                                   gggCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACa                                                                                          <  3:266649/141‑1 (MQ=255)
                                                                   gggCGTTGCTCTGGTAACCCAGGGCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTATTCGAGCTTGTGCGGGCGGCCACGCATGAAGCCAATGTCCACGGCCCGCGAAGCGCCGAACATCTACAt                                                                                   >  6:71458/1‑148 (MQ=255)
                                                                       tgggCTCTGGTAAGCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATcagca                                                                              <  8:373288/146‑1 (MQ=255)
                                                                       gTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATcagca                                                                              <  6:36016/149‑1 (MQ=255)
                                                                       gTTGCTCTGGTAACCAAGGTCGGTCATGCGCAGTCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCGCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTTCATcagca                                                                              <  8:302336/149‑1 (MQ=255)
                                                                         tGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCagcagc                                                                            >  7:368131/1‑149 (MQ=255)
                                                                          gCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCagcagc                                                                            <  3:129800/148‑1 (MQ=255)
                                                                           ctctGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAACGCAGGAACAGGCTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTACAGGGCCGGCCAGGCGCCGAACAGGTACATCAGcagcca                                                                          >  5:49466/1‑149 (MQ=255)
                                                                            tctGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGTCGGCGAGGCGCCGAACAGGTACATAAGcagccag                                                                         >  3:6383/1‑149 (MQ=255)
                                                                            tctGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGcagccag                                                                         >  7:261402/1‑149 (MQ=255)
                                                                              tGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCt                                                                       >  8:60690/1‑149 (MQ=255)
                                                                               ggTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTg                                                                      >  5:286472/1‑149 (MQ=255)
                                                                                gTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGc                                                                                                           >  3:17509/1‑111 (MQ=255)
                                                                                gTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGc                                                                                                           <  4:17509/111‑1 (MQ=255)
                                                                                    cccAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAgcg                                                                 >  6:271827/1‑149 (MQ=255)
                                                                                          tCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGaa                                                           <  7:248082/149‑1 (MQ=255)
                                                                                                aTGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTAc                                                      <  3:190800/148‑1 (MQ=255)
                                                                                                aTGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTAc                                                      <  5:326132/148‑1 (MQ=255)
                                                                                                      aGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAg                                               >  6:167079/1‑149 (MQ=255)
                                                                                                        gCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGTCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAgcg                                             <  5:167079/149‑1 (MQ=255)
                                                                                                         cTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAgcgc                                            <  1:67135/149‑1 (MQ=255)
                                                                                                          tggtggCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGcgcg                                           <  8:29147/149‑1 (MQ=255)
                                                                                                            gtggCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGTGGCGGAAGTTACGAATCAGCGCGa                                          <  7:60690/148‑1 (MQ=255)
                                                                                                                             cAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGt                        >  7:371927/1‑149 (MQ=255)
                                                                                                                                ggTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGt                     >  1:37127/1‑149 (MQ=255)
                                                                                                                                  tttCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTcg                   >  5:137349/1‑149 (MQ=255)
                                                                                                                                    tCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTcgcg                 >  5:274590/1‑149 (MQ=255)
                                                                                                                                     cGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTcgcg                 <  7:33033/148‑1 (MQ=255)
                                                                                                                                     cGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTcgcg                 <  2:349681/148‑1 (MQ=255)
                                                                                                                                      ggCATCCATCTCTTCGAGCTGGTGCGTGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGTAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTCGCgg                <  5:71458/148‑1 (MQ=255)
                                                                                                                                        gcaccAGCTCGAAGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTg                          >  6:341774/4‑136 (MQ=255)
                                                                                                                                        gcaccAGCTCGAAGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTg                          <  5:341774/133‑1 (MQ=255)
                                                                                                                                         caccAGCTCGAAGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTCGCGGTCa             >  1:108093/3‑149 (MQ=255)
                                                                                                                                          tCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTCGCGGTCAtt           <  8:368131/149‑1 (MQ=255)
                                                                                                                                            cAGCTCTTCGAGCTGGTGCGGGCTGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCTGCTGTAGCGGCGGTAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTCGCGGTCAtttt         <  3:275598/149‑1 (MQ=255)
                                                                                                                                            cAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTCGCGGTCAtttt         >  7:189233/1‑149 (MQ=255)
                                                                                                                                                   tCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTCGCGGTCATTTTCCGTACt  >  7:279848/1‑149 (MQ=255)
                                                                                                                                                    cGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTCGCGGTCATTTTCCGTACt  <  3:261922/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCACGGCCTTGCGCAGGCTGTCGGTGTAGCTGGCCAGGTTGTTGTAGCAGGGCGTAAGGCCGGCCTGGGCGTTGCTCTGGTAACCCAGGTCGCTCATGCGCAGGCTGGTGGCATAAGGCAGGAACAGGGTTTCGGCATCCAGCTCTTCGAGCTGGTGCGGGCGGCCACGCAGGAAGCCTTTGTCCAGGGCCGGCGAGGCGCCGAACAGGTACATCAGCAGCCAGCTGTAGCGGCGGAAGTTACGAATCAGCGCGATGTAGGCCGACGACTGGTAGTCGCGGTCATTTTCCGTACT  >  NC_002947/297204‑297499

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: