Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,658,123 A→G 8.8% K60E (AAG→GAG)  sahR → methionine metabolism transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,658,1230AG 8.8% 221.6 / 10.2 79K60E (AAG→GAG) sahRmethionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base G (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 = 8.86e-01

AACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGCCC  >  NC_002947/5657977‑5658269
                                                                                                                                                   |                                                                                                                                                   
aaCGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAg                                                                                                                                                   <  2:74925/149‑1 (MQ=255)
aaCGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAg                                                                                                                                                   <  6:66040/149‑1 (MQ=255)
     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                                                                                                                             <  4:66446/147‑1 (MQ=255)
                       tgCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTc                                                                                                                             <  6:380588/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:190722/1‑149 (MQ=255)
                                                               gcgcCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGt                                                                                   >  2:311082/1‑149 (MQ=255)
                                                               gcgcCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGt                                                                                   >  8:196068/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                                                          <  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                                                                      >  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)
                                                                                                                                                cTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGcc   >  6:211727/1‑149 (MQ=255)
                                                                                                                                                   aaGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGccc  <  3:38347/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGCCC  >  NC_002947/5657977‑5658269

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: