Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,658,116 T→G 8.6% H57Q (CAT→CAG sahR → methionine metabolism transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,658,1160TG 8.6% 246.8 / 13.2 81H57Q (CAT→CAGsahRmethionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base G (2/5);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTAC  >  NC_002947/5657970‑5658264
                                                                                                                                                   |                                                                                                                                                     
cccGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATc                                                                                                                                                      <  3:307324/148‑1 (MQ=255)
 ccGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCAc                                                                                                                                                    <  4:236379/149‑1 (MQ=255)
  cGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCAcc                                                                                                                                                   >  1:74925/1‑149 (MQ=255)
   gAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCt                                                                                                                                                  <  6:108779/149‑1 (MQ=255)
       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                                                                                                                       >  1:224185/1‑128 (MQ=255)
                                                  gctgcgcCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCg                                                                                                                       <  2:224185/128‑1 (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                                                                              >  2:311082/1‑149 (MQ=255)
                                                                      gcgcCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGt                                                                              >  8:190722/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                                                       <  1:364241/149‑1 (MQ=255)
                                                                                             ggCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCt                                                       <  7:196068/149‑1 (MQ=255)
                                                                                             ggCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCt                                                       <  7:190722/149‑1 (MQ=255)
                                                                                               cGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                     <  7:267700/149‑1 (MQ=255)
                                                                                               cGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                     <  6:356934/149‑1 (MQ=255)
                                                                                               cGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                     <  6:362121/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                                                     <  1:304762/134‑1 (MQ=255)
                                                                                                              gcAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                     >  2:304762/1‑134 (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)
                                                                                                                                                   |                                                                                                                                                     
CCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTAC  >  NC_002947/5657970‑5658264

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: