Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,938,664 T→G 6.6% S409A (TCG→GCG)  rbbA → ribosome‑associated ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,938,6640TG 6.6% 275.0 / 7.8 91S409A (TCG→GCG) rbbAribosome‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base T (44/41);  new base G (4/2);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GCAGAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTCCAC  >  NC_002947/5938516‑5938809
                                                                                                                                                    |                                                                                                                                                                     
gCAGAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGt                                                                                                                                                                       >  7:167917/1‑149 (MQ=255)
   gAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGc                                                                                                                                                                    <  3:193213/149‑1 (MQ=255)
    aaCCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCt                                                                                                                                                                   >  6:367366/1‑149 (MQ=255)
    aaCCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACAGGCGCt                                                                                                                                                                   >  7:225650/1‑149 (MQ=255)
      ccTGGCCCTGCCTGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGCCCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTgg                                                                                                                                                                 >  7:174566/1‑149 (MQ=255)
      ccTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTgg                                                                                                                                                                 >  2:297901/1‑149 (MQ=255)
          gCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACAGGCGCTGGcg                                                                                                                                                               <  8:225650/147‑2 (MQ=255)
           cccTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTg                                                                                                                                                       <  5:117416/149‑1 (MQ=255)
             cTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑Gc                                                                                                                                                <  8:167917/149‑1 (MQ=255)
                cATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑g‑‑‑‑‑g                                                                                                                                    <  7:163861/148‑1 (MQ=255)
                cATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑g‑‑‑‑‑g                                                                                                                                    <  3:234861/148‑1 (MQ=255)
                 aTGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑Gtt                                                                                                                                  <  5:158102/149‑1 (MQ=255)
                  tGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTc                                                                                                                                 <  5:367366/149‑1 (MQ=255)
                      cGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCa                                                                                                                             <  2:149685/149‑1 (MQ=255)
                           gTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGcc                                                                                                                        >  8:338784/1‑149 (MQ=255)
                              cGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCg                                                                                                                      <  6:319376/148‑1 (MQ=255)
                               gACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCgg                                                                                                                     >  6:404749/1‑148 (MQ=255)
                               gACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCgg                                                                                                                     <  7:15072/148‑1 (MQ=255)
                                aCTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGaa                                                                                                                   <  2:12363/149‑1 (MQ=255)
                                  ttGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGt                                                                                                                 >  1:106403/1‑149 (MQ=255)
                                    gCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCACCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCt                                                                                                               >  2:185642/1‑149 (MQ=255)
                                     cccAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCtt                                                                                                              >  2:143752/1‑149 (MQ=255)
                                     cccAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCt                                                                                                               <  5:254208/148‑1 (MQ=255)
                                      ccAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTg                                                                                                             <  2:168432/149‑1 (MQ=255)
                                      ccAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCACCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCtt                                                                                                              >  7:215568/1‑148 (MQ=255)
                                       cAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGa                                                                                                            <  6:206642/149‑1 (MQ=255)
                                       cAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGa                                                                                                            >  5:140999/1‑149 (MQ=255)
                                        aaaGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGAt                                                                                                           >  6:278267/1‑149 (MQ=255)
                                            ccGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCGGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTc                                                                                                       <  2:107620/149‑1 (MQ=255)
                                             cGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTc                                                                                                       >  5:15824/1‑148 (MQ=255)
                                               agagCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTcgac                                                                                                    <  3:100049/149‑1 (MQ=255)
                                                 agCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTcgacga                                                                                                  <  5:278267/149‑1 (MQ=255)
                                                      ccAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGa                                                                                             >  1:59743/1‑149 (MQ=255)
                                                       cAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCACCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGAc                                                                                            <  8:215568/149‑1 (MQ=255)
                                                            cATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCg                                                                                       <  1:297901/149‑1 (MQ=255)
                                                                gATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGt                                                                                   <  3:17043/149‑1 (MQ=255)
                                                                 aggAGCTGAGCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTc                                                                                  <  6:219634/147‑1 (MQ=255)
                                                                 aTGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTc                                                                                  <  6:15824/149‑1 (MQ=255)
                                                                     gCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGAcc                                                                              >  4:366874/1‑149 (MQ=255)
                                                                         aTCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCccgc                                                                          >  6:225773/1‑149 (MQ=255)
                                                                              gcgcTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGccgccc                                                                     >  1:138688/1‑149 (MQ=255)
                                                                                gcTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGccgcccgc                                                                   >  8:305830/1‑149 (MQ=255)
                                                                                 cTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCcgcg                                                                  <  5:225773/149‑1 (MQ=255)
                                                                                  ttCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAG‑‑‑‑GNGGGGGTGCGGGAGG‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGAcccc                                                                            >  1:141978/1‑149 (MQ=255)
                                                                                  ttCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGAATGAGGGGGT‑‑‑‑GC‑‑‑‑‑GATAGAGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGAcccc                                                                            >  1:136748/1‑149 (MQ=255)
                                                                                      aCCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCgacgac                                                             >  7:132563/1‑149 (MQ=255)
                                                                                       ccTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACt                                                            >  5:136958/1‑149 (MQ=255)
                                                                                        cTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACtt                                                           >  4:186370/1‑149 (MQ=255)
                                                                                           gccgcTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCt                                                         <  5:404749/148‑1 (MQ=255)
                                                                                            ccgcTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTgg                                                       <  7:338784/149‑1 (MQ=255)
                                                                                               cTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGgcgc                                                    <  1:143752/149‑1 (MQ=255)
                                                                                                    gCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCctgct                                               >  7:114700/1‑149 (MQ=255)
                                                                                                    gCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCctgct                                               >  6:192227/1‑149 (MQ=255)
                                                                                                     ccGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCctgctg                                              >  5:315179/1‑149 (MQ=255)
                                                                                                       gACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTgg                                             >  3:353165/1‑148 (MQ=255)
                                                                                                        aCCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTa                                           <  2:138688/149‑1 (MQ=255)
                                                                                                          cAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGCCCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAg                                          >  8:157519/1‑148 (MQ=255)
                                                                                                          cAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAg                                          >  6:297878/1‑148 (MQ=255)
                                                                                                           aaCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGaa                                        >  1:15258/1‑149 (MQ=255)
                                                                                                           aaCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGaa                                        >  1:34692/1‑149 (MQ=255)
                                                                                                             ccGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACt                                      <  5:325308/149‑1 (MQ=255)
                                                                                                             ccGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACt                                      <  2:106403/149‑1 (MQ=255)
                                                                                                                   ggCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCcg                                >  5:367935/1‑149 (MQ=255)
                                                                                                                     cgcCCTGCCCCTTGGTTTGCGCGAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑TTTCTGCATCGCGCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCTACGACTTCTGGCGCCTGCTGGTAGAACTGTCCcgc                               <  8:217594/148‑1 (MQ=255)
                                                                                                                            ccccTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAgg                        <  6:136958/148‑1 (MQ=255)
                                                                                                                            ccccTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAgg                        <  2:15258/148‑1 (MQ=255)
                                                                                                                              ccTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCg                     >  3:257007/1‑149 (MQ=255)
                                                                                                                               cTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAAc                                       >  8:299041/1‑130 (MQ=255)
                                                                                                                               cTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAAc                                       <  7:299041/130‑1 (MQ=255)
                                                                                                                               cTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGt                    <  6:140999/149‑1 (MQ=255)
                                                                                                                               cTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGt                    <  2:239398/149‑1 (MQ=255)
                                                                                                                               cTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGt                    >  5:81356/1‑149 (MQ=255)
                                                                                                                                  ggTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGa                  <  8:132563/148‑1 (MQ=255)
                                                                                                                                   gTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGAc                 <  1:118589/148‑1 (MQ=255)
                                                                                                                                    tttGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCa               >  3:6998/1‑149 (MQ=255)
                                                                                                                                      actgcCAGCGACAGGCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCAtt             >  7:218244/4‑149 (MQ=255)
                                                                                                                                       ctgcCAGCGACAGGCGCTGGCG‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCAttt            >  5:432487/3‑149 (MQ=255)
                                                                                                                                        tgcCAGCGACAGGCGCTGGCG‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGa                  <  8:302019/141‑1 (MQ=255)
                                                                                                                                        tgcCAGCGACAGGCGCTGGCG‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGa                  >  7:302019/2‑142 (MQ=255)
                                                                                                                                        cgcCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCAtttt           >  4:107023/1‑149 (MQ=255)
                                                                                                                                         gcCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTccc                                 >  7:178624/1‑126 (MQ=255)
                                                                                                                                         gcCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTccc                                 <  8:178624/126‑1 (MQ=255)
                                                                                                                                         gcCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTc          <  6:81356/149‑1 (MQ=255)
                                                                                                                                           cAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGCCCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCct        >  4:127464/1‑149 (MQ=255)
                                                                                                                                           cAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCct        <  7:305830/149‑1 (MQ=255)
                                                                                                                                           cAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCct        >  7:139607/1‑149 (MQ=255)
                                                                                                                                           cAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCct        >  5:47099/1‑149 (MQ=255)
                                                                                                                                             gcgcCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCctct      >  5:405870/1‑149 (MQ=255)
                                                                                                                                                ccTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCCTTTTCCTCTcca   >  4:315450/1‑149 (MQ=255)
                                                                                                                                                 cTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTccac  >  2:207053/1‑149 (MQ=255)
                                                                                                                                                   gTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGgc                                                      >  6:55795/1‑95 (MQ=255)
                                                                                                                                                   gTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGgc                                                      <  5:55795/95‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                                     
GCAGAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCA‑‑‑‑‑GTG‑‑‑‑‑GC‑‑‑‑‑G‑‑‑‑‑GTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTCCAC  >  NC_002947/5938516‑5938809

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: