Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 285,287 G→C 5.5% Q55E (CAG→GAG)  atsK ← alpha‑ketoglutarate‑dependent sulfonate dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947285,2870GC 5.5% 233.4 / 2.6 72Q55E (CAG→GAG) atsKalpha‑ketoglutarate‑dependent sulfonate dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base C (2/2);  total (34/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.01e-01

AAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAGTGCTGGGCTTAGGGGA  >  NC_002947/285143‑285432
                                                                                                                                                |                                                                                                                                                 
aaaaGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGt                                                                                                                                               >  2:211876/1‑149 (MQ=255)
      tCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCgg                                                                                                                                          <  1:398522/148‑1 (MQ=255)
         cGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATGGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTg                                                                                                                                      <  8:29365/149‑1 (MQ=255)
         cGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTg                                                                                                                                      <  4:16459/149‑1 (MQ=255)
             ggTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCg                                                                                                                                  <  7:174841/149‑1 (MQ=255)
                   ccACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGt                                                                                                                            >  3:222681/1‑149 (MQ=255)
                    cacaCGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTc                                                                                                                           >  6:230230/1‑149 (MQ=255)
                     acacGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTCCTGTTCCGGGTTGATCGGCTGGTcg                                                                                                                          >  1:267964/1‑149 (MQ=255)
                     acacGGCGTTGTCGCGCACATCGGTAACCGCCGTGACGAGCACAAGCACCTGTGGGGTGTCTGGCACGGTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTCCTGTGCCGGGTTGCTCGGCTGGTcg                                                                                                                          >  1:273858/1‑149 (MQ=255)
                       acGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTcgcg                                                                                                                        >  8:140986/1‑149 (MQ=255)
                         ggCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTcgcg                                                                                                                        >  5:79695/1‑147 (MQ=255)
                             ttGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGaa                                                                                                                  <  3:257272/149‑1 (MQ=255)
                               gTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAAc                                                                                                                 <  1:316833/148‑1 (MQ=255)
                               gTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAAc                                                                                                                 <  6:385760/148‑1 (MQ=255)
                                tCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACa                                                                                                                <  3:303959/148‑1 (MQ=255)
                                            gTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTg                                                                                                   >  4:211233/1‑149 (MQ=255)
                                                 cgccgTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCa                                                                                              <  1:211876/149‑1 (MQ=255)
                                                 cgccgGGTCGAGCACAAGCCCCTGTGGGGTGTCTGGCACGTGGGGGTAACGCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGGCGCGGAGGAACAGCACCTGGGGCTGGAGCa                                                                                              <  6:209561/149‑1 (MQ=255)
                                                   ccgTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTCCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGt                                                                                            <  2:267964/149‑1 (MQ=255)
                                                   ccgTGTCGAGCACAAGCACCTGTGGGGTGGCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCTGCAAAGCGCGCCTCCTGGTCCGGGTGGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGt                                                                                            <  2:273858/149‑1 (MQ=255)
                                                     gtgtCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGc                                                                                          <  7:216446/149‑1 (MQ=255)
                                                      tgtCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGcc                                                                                         <  1:50708/149‑1 (MQ=255)
                                                           aGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTc                                                                                    >  5:193050/1‑149 (MQ=255)
                                                                aaGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCggag                                                                                                                     <  2:135485/111‑1 (MQ=255)
                                                                aaGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCAcc                                                                                                           <  2:228112/121‑1 (MQ=255)
                                                                aaGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCAcc                                                                                                           >  1:228112/1‑121 (MQ=255)
                                                                aaGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGc                                                                                          <  4:225259/138‑1 (MQ=255)
                                                                aaGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGc                                                                                          >  3:225259/1‑138 (MQ=255)
                                                                aaGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATgg                                                                               <  6:193050/149‑1 (MQ=255)
                                                                aaGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATgg                                                                               <  7:140986/149‑1 (MQ=255)
                                                                      cTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCg                                                                                   <  4:244135/139‑1 (MQ=255)
                                                                      cTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCg                                                                                   >  3:244135/1‑139 (MQ=255)
                                                                           gggTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGctct                                                                    <  2:90216/149‑1 (MQ=255)
                                                                             gtgtCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCtcttc                                                                  >  1:398038/1‑149 (MQ=255)
                                                                             gtgtCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCtcttc                                                                  <  5:230230/149‑1 (MQ=255)
                                                                                tCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTg                                                                                                      >  1:319255/1‑110 (MQ=255)
                                                                                tCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTg                                                                                                      <  2:319255/110‑1 (MQ=255)
                                                                                  tGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTg                                                                                                   >  7:278269/1‑111 (MQ=255)
                                                                                  tGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTg                                                                                                   <  8:278269/111‑1 (MQ=255)
                                                                                  tGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTgcgc                                                             >  3:107264/1‑149 (MQ=255)
                                                                                  tGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTgcgc                                                             >  1:197601/1‑149 (MQ=255)
                                                                                  tGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTgcgc                                                             >  7:390650/1‑149 (MQ=255)
                                                                                  tGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGACTGTTCGATGGCATCGCGCTCTTCTgcgc                                                             >  5:275912/1‑149 (MQ=255)
                                                                                   ggCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCt                                                            >  8:43774/1‑149 (MQ=255)
                                                                                    gCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGtt                                                                                                      >  8:260319/1‑105 (MQ=255)
                                                                                    gCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTgcgc                                                             <  6:79695/147‑1 (MQ=255)
                                                                                    gCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTg                                                           <  7:43774/149‑1 (MQ=255)
                                                                                    gCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTg                                                           <  3:211233/149‑1 (MQ=255)
                                                                                         ttGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATc                                                      >  2:284037/1‑149 (MQ=255)
                                                                                           ggggTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGaa                                                                                                                  <  8:388478/87‑1 (MQ=255)
                                                                                           ggggTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGaa                                                                                                                  >  7:388478/1‑87 (MQ=255)
                                                                                                             gtgCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCAcgccgc                                  >  2:412821/1‑149 (MQ=255)
                                                                                                              ggCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCAcgccgc                                  <  6:374921/147‑1 (MQ=255)
                                                                                                                      gCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTggg                         >  6:164547/1‑149 (MQ=255)
                                                                                                                      gCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTggg                         >  2:212833/1‑149 (MQ=255)
                                                                                                                      gCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTgg                          >  4:44845/1‑148 (MQ=255)
                                                                                                                      gCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTgg                          <  3:44845/148‑1 (MQ=255)
                                                                                                                       ccGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGgc                        >  4:127909/1‑149 (MQ=255)
                                                                                                                            gCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGcc                     >  4:61966/1‑147 (MQ=255)
                                                                                                                            gCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGa                   <  8:390650/149‑1 (MQ=255)
                                                                                                                             cGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAg                  >  3:233881/1‑149 (MQ=255)
                                                                                                                             cGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAg                  <  7:266217/149‑1 (MQ=255)
                                                                                                                              gTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAGt                 <  3:127909/149‑1 (MQ=255)
                                                                                                                              gTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAGt                 >  6:157336/1‑149 (MQ=255)
                                                                                                                                 cGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAGTGCt              >  7:239095/1‑149 (MQ=255)
                                                                                                                                  ggCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAGTGCTg             <  8:239095/149‑1 (MQ=255)
                                                                                                                                  ggCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAGTGCTg             >  1:404758/1‑149 (MQ=255)
                                                                                                                                  ggCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCCGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAGTGCTg             <  4:107264/149‑1 (MQ=255)
                                                                                                                                         cgcgcCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTgg                          >  2:45535/1‑129 (MQ=255)
                                                                                                                                         cgcgcCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTgg                          <  1:45535/129‑1 (MQ=255)
                                                                                                                                          gcgcCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAGTGCTGGGCTTAgg     <  1:284037/149‑1 (MQ=255)
                                                                                                                                             ccTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAGTGCTGGGCTTAGGGGa  <  6:275912/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGGTGGATGTGCAGGTCGCCGAAGCGTGCGGCAAAGCGCGCCTGCTGTTCCGGGTTGATCGGCTGGTCGCGGAGGAACAGCACCTGGTGCTGGAGCAGTGCCTGTTCGATGGCATCGCGCTCTTCTGCGCTGATATCGCGGCTGATGTCCACGCCGCTGATCTGGGCGCCGAGTGCTGGGCTTAGGGGA  >  NC_002947/285143‑285432

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