Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,658,129 C→A 10.8% L62M (CTG→ATG)  sahR → methionine metabolism transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,658,1290CA10.8% 250.9 / 16.4 83L62M (CTG→ATG) sahRmethionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base A (3/6);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

GAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGCCCGATGTG  >  NC_002947/5657982‑5658275
                                                                                                                                                     |                                                                                                                                                  
gAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCt                                                                                                                                                   >  6:123396/1‑149 (MQ=255)
gAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCt                                                                                                                                                   >  3:109850/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                                                                                 <  1:196180/149‑1 (MQ=255)
                                                                  gCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGtgccctgcc                                                                                 <  3:293821/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                                                                  <  7:196068/149‑1 (MQ=255)
                                                                                 ggCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCt                                                                  <  7:190722/149‑1 (MQ=255)
                                                                                 ggCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCt                                                                  <  1:364241/149‑1 (MQ=255)
                                                                                   cGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                                <  6:362121/149‑1 (MQ=255)
                                                                                   cGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                                <  6:356934/149‑1 (MQ=255)
                                                                                   cGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                                <  7:267700/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                                                                            >  4:59483/1‑117 (MQ=255)
                                                                                                       tCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATg                                                                            <  3:59483/117‑1 (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)
                                                                                                                                          ccttgaGGTGATGGCTCATGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCgg                                                        <  2:338300/98‑1 (MQ=255)
                                                                                                                                          ccttgaGGTGATGGCTCATGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCgg                                                        >  1:338300/6‑103 (MQ=255)
                                                                                                                                           cTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGcc         >  6:211727/1‑149 (MQ=255)
                                                                                                                                              aaGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGccc        <  3:38347/147‑1 (MQ=255)
                                                                                                                                                ggTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACgcgg                                              >  2:232427/1‑107 (MQ=255)
                                                                                                                                                ggTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACgcgg                                              <  1:232427/107‑1 (MQ=255)
                                                                                                                                                  tGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGCCCGAtgtg  <  3:143052/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGCCCGATGTG  >  NC_002947/5657982‑5658275

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: