Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 132,883 A→C 7.0% E102D (GAA→GAC)  dsbA → thiol:disulfide interchange protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947132,8830AC 7.0% 258.7 / 7.6 84E102D (GAA→GAC) dsbAthiol:disulfide interchange protein
Reads supporting (aligned to +/- strand):  ref base A (43/35);  new base C (3/3);  total (46/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.55e-01

GGCTGCCCACACTGCTACCACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCG  >  NC_002947/132743‑133030
                                                                                                                                            |                                                                                                                                                   
ggCTGCGGACGCTGCTAGCACTTGGAGCCGACAATCAACGCTTGGCCTGAATAGCTGCCGAAGGACGTCAACTTCAAGCGGGTACCGGCCATGTTCGGTGGCCCGTGGGATGGAGACGGCCTGATGTTCCTGACCCTCGACGcccatgg                                                                                                                                             >  6:244446/1‑144 (MQ=255)
       ctcacTGCTACAACTTCTAGCCGACCATCAACCCTAGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCTGGGCATGCACACGGCCAGATGCTCCTGACCCTCGACGCCCATGGCTTCGAg                                                                                                                                      <  5:244446/147‑1 (MQ=255)
       cacaCTGCTACCACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTAGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAg                                                                                                                                      >  5:110480/1‑149 (MQ=255)
        acacTGCTACCACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAg                                                                                                                                      >  4:24082/1‑148 (MQ=255)
        acacTGCTACCACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAg                                                                                                                                      >  5:235862/1‑148 (MQ=255)
         cacTGCTACCACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGca                                                                                                                                    >  7:354859/1‑149 (MQ=255)
            tGCTACCACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGcac                                                                                                                                   >  3:352579/1‑147 (MQ=255)
              cTACCACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAg                                                                                                                                <  8:139189/148‑1 (MQ=255)
               tACCACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAg                                                                                                                                      <  6:247362/141‑1 (MQ=255)
               tACCACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAg                                                                                                                                      >  5:247362/1‑141 (MQ=255)
               tACCACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGt                                                                                                                              >  5:388085/1‑149 (MQ=255)
                 ccACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGtt                                                                                                                             >  5:22365/1‑148 (MQ=255)
                  cACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGCCCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCa                                                                                                                           >  6:123383/1‑149 (MQ=255)
                   aCTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCAt                                                                                                                          <  4:331392/149‑1 (MQ=255)
                   aCTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCAt                                                                                                                          <  4:85516/149‑1 (MQ=255)
                   aCTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCAt                                                                                                                          >  1:167859/1‑149 (MQ=255)
                    cTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCAtg                                                                                                                         >  4:322603/1‑149 (MQ=255)
                     ttCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCAtgc                                                                                                                        >  7:132633/1‑149 (MQ=255)
                        gAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCAtgctgc                                                                                                                     <  2:280225/149‑1 (MQ=255)
                        gAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCAtgctgc                                                                                                                     >  5:285589/1‑149 (MQ=255)
                         aGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCg                                                                                                                    >  2:39512/1‑149 (MQ=255)
                          gCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCgg                                                                                                                   >  5:44836/1‑149 (MQ=255)
                           ccGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGgt                                                                                                                  >  2:248051/1‑149 (MQ=255)
                              aCCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCgttgtt                                                                                                               <  8:132633/149‑1 (MQ=255)
                              aCCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGtt                                                                                                               <  2:194107/149‑1 (MQ=255)
                                 aTCAACCCTTGGGCTGAAAAGCTGCGGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGACGCCCATGGCTTCGAg                                                                                                                                      <  1:210613/123‑1 (MQ=255)
                                 aTCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGACGCCCATGGCTTCGAg                                                                                                                                      >  2:210613/1‑123 (MQ=255)
                                  tCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGAt                                                                                                           <  6:22365/149‑1 (MQ=255)
                                  tCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCAAGTTAGGTGGCCCGTGGGATGCACACGGCCAGATGTTACTGACCCTCGAACCCATGGGCGTCCAGCACAAGGTTCATGCTGCGGTGTTCGAt                                                                                                           >  1:305674/1‑149 (MQ=255)
                                    aaCCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGaa                                                                                                                                                     <  7:99973/105‑1 (MQ=255)
                                    aaCCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGaa                                                                                                                                                     >  8:99973/1‑105 (MQ=255)
                                    aaCCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATg                                                                                                          <  1:15233/148‑1 (MQ=255)
                                       ccTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCAt                                                                                                      >  3:106896/1‑149 (MQ=255)
                                         ttGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGACGCCCATGGCTTCGAgg                                                                                                                                     >  8:209921/1‑115 (MQ=255)
                                         ttGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGACGCCCATGGCTTCGAgg                                                                                                                                     <  7:209921/116‑2 (MQ=255)
                                         ttGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCAtg                                                                                                                         <  5:120987/128‑1 (MQ=255)
                                         ttGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCAtg                                                                                                                         >  6:120987/1‑128 (MQ=255)
                                           gggCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATcc                                                                                                    <  8:354859/147‑1 (MQ=255)
                                              cTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAAc                                                                                               <  6:285589/149‑1 (MQ=255)
                                              cTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAAc                                                                                               <  6:44836/149‑1 (MQ=255)
                                                    aGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGcac                                                                                         <  6:235862/149‑1 (MQ=255)
                                                    aGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGca                                                                                          >  5:197580/1‑148 (MQ=255)
                                                      cTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACaa                                                                                       <  4:350645/149‑1 (MQ=255)
                                                           gAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTc                                                                                  >  6:150303/1‑149 (MQ=255)
                                                           gAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTc                                                                                  <  1:176811/149‑1 (MQ=255)
                                                            aaGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCt                                                                                 <  3:24082/149‑1 (MQ=255)
                                                             aGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTGGATGCCATCCAGAACCAGCACAAGCGTCt                                                                                 >  8:336634/1‑148 (MQ=255)
                                                              ggACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGa                                                                               <  2:167859/149‑1 (MQ=255)
                                                                aCGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGAcc                                                                             >  1:149516/1‑149 (MQ=255)
                                                                 cGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCg                                                                            <  3:322603/149‑1 (MQ=255)
                                                                  gTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGa                                                                           >  7:280340/1‑149 (MQ=255)
                                                                    cAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCAt                                                                                                                          >  5:241123/1‑100 (MQ=255)
                                                                    cAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCAt                                                                                                                          <  6:241123/100‑1 (MQ=255)
                                                                          cAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAAc                                                                    <  7:176032/148‑1 (MQ=255)
                                                                                       gCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACtt                                                      >  5:396486/1‑149 (MQ=255)
                                                                                         cATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTcc                                                    >  3:111112/1‑149 (MQ=255)
                                                                                          aTGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCt                                                   >  3:102195/1‑149 (MQ=255)
                                                                                          aTGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCt                                                   >  3:107630/1‑149 (MQ=255)
                                                                                          aTGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCt                                                   >  5:169547/1‑149 (MQ=255)
                                                                                           tGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTc                                                  <  5:150303/149‑1 (MQ=255)
                                                                                             ttCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGGCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGTCTTCGTCGc                                                >  6:128106/1‑149 (MQ=255)
                                                                                             ttCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGc                                                >  5:21930/1‑149 (MQ=255)
                                                                                             ttCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGc                                                <  1:248051/149‑1 (MQ=255)
                                                                                             ttCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGACGACTTCCTCGc                                                <  1:39512/149‑1 (MQ=255)
                                                                                                     cccGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAgg                                                                                                                               >  5:238047/1‑62 (MQ=255)
                                                                                                     cccGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAgg                                                                                                                               <  6:238047/62‑1 (MQ=255)
                                                                                                     cccGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAg                                        <  6:110480/149‑1 (MQ=255)
                                                                                                      ccGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAgg                                       >  2:137933/1‑149 (MQ=255)
                                                                                                      ccGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAgg                                       >  2:204641/1‑149 (MQ=255)
                                                                                                      ccGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAgg                                       >  8:150235/1‑149 (MQ=255)
                                                                                                        gTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGgtg                                     >  2:58804/1‑149 (MQ=255)
                                                                                                         tGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGgtgt                                    >  6:62921/1‑149 (MQ=255)
                                                                                                           ggATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCTTGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGtgtg                                   <  6:388085/148‑1 (MQ=255)
                                                                                                           ggATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGtgtg                                   <  4:106896/148‑1 (MQ=255)
                                                                                                           ggATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCGCGCCACTCAAGGtgtg                                   <  4:111112/148‑1 (MQ=255)
                                                                                                                cacaCGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTggacaa                              <  6:21930/148‑1 (MQ=255)
                                                                                                                 acacGGCCAGATGTTCCTGACCCTCGAAGACATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTggacaag                             <  6:396486/148‑1 (MQ=255)
                                                                                                                   acGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTggacaagga                           >  4:242980/1‑148 (MQ=255)
                                                                                                                      gCCAGAGGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACACGGTGCCTGCTGCGGTGTGCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGGGgacaaggacaag                       <  6:182105/149‑1 (MQ=255)
                                                                                                                            tGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTc                 <  2:149516/149‑1 (MQ=255)
                                                                                                                              ttCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGc               >  6:227433/1‑149 (MQ=255)
                                                                                                                                 cTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCAc            >  4:282640/1‑149 (MQ=255)
                                                                                                                                 cTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGT‑TTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCAcc           >  1:392970/1‑149 (MQ=255)
                                                                                                                                    aCCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCtt         <  6:169547/149‑1 (MQ=255)
                                                                                                                                     cccTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTc        >  1:113989/1‑149 (MQ=255)
                                                                                                                                          gAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGAc     <  3:242980/147‑1 (MQ=255)
                                                                                                                                           aaGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGACGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGACACCTTCGACTCg  >  6:231051/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GGCTGCCCACACTGCTACCACTTCGAGCCGACCATCAACCCTTGGGCTGAAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCG  >  NC_002947/132743‑133030

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