Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,658,122 C→T 8.8% L59L (CTC→CTT sahR → methionine metabolism transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,658,1220CT 8.8% 225.9 / 9.0 79L59L (CTC→CTTsahRmethionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (37/35);  new base T (2/5);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

AACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGCC  >  NC_002947/5657977‑5658268
                                                                                                                                                  |                                                                                                                                                   
aaCGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAg                                                                                                                                                  <  2:74925/149‑1 (MQ=255)
aaCGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAg                                                                                                                                                  <  6:66040/149‑1 (MQ=255)
     gAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCt                                                                                                                                            >  3:109850/1‑149 (MQ=255)
     gAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCt                                                                                                                                            >  6:123396/1‑149 (MQ=255)
      aaaCATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTg                                                                                                                                           <  8:30074/149‑1 (MQ=255)
        aCATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCAGCACCTTGAGGTGATGGc                                                                                                                                         <  1:287173/149‑1 (MQ=255)
             ggCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGCCTGAGCCAGCACCTTGAGGTGATGGCTCatg                                                                                                                                    <  1:23237/148‑3 (MQ=38)
              gCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCCTGAGCCAGCACCTTGAGGTGATGGCTCATGc                                                                                                                                   <  7:260018/149‑1 (MQ=39)
               cAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGc                                                                                                                                   <  1:91816/149‑1 (MQ=255)
                   ccTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGAcc                                                                                                                              >  4:293821/1‑149 (MQ=255)
                       tgCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTc                                                                                                                            <  6:380588/147‑1 (MQ=255)
                       tgCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTc                                                                                                                            <  4:66446/147‑1 (MQ=255)
                        gCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTg                                                                                                                         >  2:281354/1‑149 (MQ=255)
                        gCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCAGCACCTTGAGGTGATGGc                                                                                                                                         >  6:4306/1‑133 (MQ=255)
                        gCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCAGCACCTTGAGGTGATGGc                                                                                                                                         <  5:4306/133‑1 (MQ=255)
                         cAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCAGGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTg                                                                                                                         <  3:331723/148‑1 (MQ=255)
                          aaaGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTgg                                                                                                                        >  3:103327/1‑148 (MQ=255)
                           aaGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGc                                                                                                                       >  7:189651/1‑148 (MQ=255)
                                aGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACccg                                                                                                                 <  5:133919/149‑1 (MQ=255)
                                 gTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACccgc                                                                                                                >  5:106649/1‑149 (MQ=255)
                                          cgctgcgcCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAgg                                                                                                        <  5:123396/148‑1 (MQ=255)
                                           gctgcgcCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCg                                                                                                                           <  2:224185/128‑1 (MQ=255)
                                           gctgcgcCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCg                                                                                                                           >  1:224185/1‑128 (MQ=255)
                                               cgcCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAAc                                                                                                    <  6:347842/147‑1 (MQ=255)
                                                  cTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCa                                                                                                >  2:364241/1‑148 (MQ=255)
                                                  cTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCAt                                                                                               >  4:65330/1‑149 (MQ=255)
                                                     aaCGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCAtctt                                                                                            <  6:348926/149‑1 (MQ=255)
                                                     aaCGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCAtct                                                                                             >  2:196180/1‑148 (MQ=255)
                                                     aaCGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCAtc                                                                                              <  6:106649/147‑1 (MQ=255)
                                                       cGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCAtcttct                                                                                          <  5:88569/149‑1 (MQ=255)
                                                         tATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTCGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTa                                                                                         <  1:281354/149‑1 (MQ=255)
                                                             gcgcgcCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTAccgcc                                                                                    >  4:349822/1‑149 (MQ=255)
                                                               gcgcCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGt                                                                                  >  8:196068/1‑149 (MQ=255)
                                                               gcgcCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGt                                                                                  >  2:311082/1‑149 (MQ=255)
                                                               gcgcCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGt                                                                                  >  8:190722/1‑149 (MQ=255)
                                                                 gcCTTGTCCAATGACTCGTTCGGAGTGCTGTAACTTACGCAGATCTACGACATCGGCCAGTAGGGCATGAGCCATCACCTCAACTTGCTGTATCAGTCCGACATCGTGGCAACCCTCCGCGAAGGCAAAGCAATCTAATACAGCCGtgc                                                                                >  8:158338/1‑149 (MQ=255)
                                                                    ttGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGtgccct                                                                             >  5:337866/1‑149 (MQ=255)
                                                                       gCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGtgccctgcc                                                                          <  3:293821/149‑1 (MQ=255)
                                                                       gCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGtgccctgcc                                                                          <  1:196180/149‑1 (MQ=255)
                                                                             gACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATgg                                                                    >  4:337280/1‑149 (MQ=255)
                                                                              aCTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGc                                                                   >  8:267700/1‑149 (MQ=255)
                                                                               cTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGcc                                                                  >  2:350424/1‑149 (MQ=255)
                                                                                  gTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGCCCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTgc                                                               >  1:187497/1‑149 (MQ=255)
                                                                                      ggCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCt                                                           <  1:364241/149‑1 (MQ=255)
                                                                                      ggCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCt                                                           <  7:190722/149‑1 (MQ=255)
                                                                                      ggCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCt                                                           <  7:196068/149‑1 (MQ=255)
                                                                                        cGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                         <  7:267700/149‑1 (MQ=255)
                                                                                        cGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                         <  6:362121/149‑1 (MQ=255)
                                                                                        cGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                         <  6:356934/149‑1 (MQ=255)
                                                                                          tGCTGGAACTGGCGCAGATCTTCGACATCGTCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGt                                                       <  8:189651/149‑1 (MQ=255)
                                                                                          tGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGt                                                       <  4:103327/149‑1 (MQ=255)
                                                                                            cTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTgg                                                     >  5:32289/1‑149 (MQ=255)
                                                                                               gAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCg                                                  <  1:350424/149‑1 (MQ=255)
                                                                                                       gcAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                         >  2:304762/1‑134 (MQ=255)
                                                                                                       gcAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                         <  1:304762/134‑1 (MQ=255)
                                                                                                       gcAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCAcg                                          >  8:158507/1‑149 (MQ=255)
                                                                                                        cAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCAcgc                                         <  6:32289/149‑1 (MQ=255)
                                                                                                           aTCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACgcggc                                      >  4:143052/1‑149 (MQ=255)
                                                                                                            tCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATg                                                                     <  3:59483/117‑1 (MQ=255)
                                                                                                            tCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATg                                                                     >  4:59483/1‑117 (MQ=255)
                                                                                                                       cGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCg                              >  8:152333/1‑146 (MQ=255)
                                                                                                                        ggCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaag                           >  6:330417/1‑148 (MQ=255)
                                                                                                                        ggCCAGGCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaag                           <  7:158338/148‑1 (MQ=255)
                                                                                                                         gCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaagaa                         <  5:13301/149‑1 (MQ=255)
                                                                                                                         gCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaag                           >  7:175968/1‑147 (MQ=255)
                                                                                                                            aGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTc                      <  4:109850/149‑1 (MQ=255)
                                                                                                                                  gCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCgacgac                >  4:38347/1‑149 (MQ=255)
                                                                                                                                     tGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaa                            >  1:185308/1‑134 (MQ=255)
                                                                                                                                     tGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaa                            <  2:185308/134‑1 (MQ=255)
                                                                                                                                      gAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTgg            >  4:315281/1‑149 (MQ=255)
                                                                                                                                        gCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTTCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTggac          >  7:180532/1‑147 (MQ=255)
                                                                                                                                        gCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGcc          <  3:315281/149‑1 (MQ=255)
                                                                                                                                         ccATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGTAGGCAACGCCCTCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGccc         >  1:90964/1‑149 (MQ=255)
                                                                                                                                         ccATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaa                            >  7:269013/1‑130 (MQ=255)
                                                                                                                                         ccATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaa                            <  8:269013/130‑1 (MQ=255)
                                                                                                                                         ccAGCACCTTGAGGTGATGGCTCAGGCCGACCTCGTGGCAACGCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTc                      <  8:346760/136‑1 (MQ=255)
                                                                                                                                         ccAGCACCTTGAGGTGATGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTc                      >  7:346760/1‑136 (MQ=255)
                                                                                                                                            tCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTAc      >  1:150133/1‑149 (MQ=255)
                                                                                                                                                cTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGcc  >  6:211727/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGCC  >  NC_002947/5657977‑5658268

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: