Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,306,923 A→C 6.9% I232S (ATC→AGC)  livK ← branched‑chain amino acids ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,306,9230AC 6.9% 269.9 / 9.7 87I232S (ATC→AGC) livKbranched‑chain amino acids ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (41/40);  new base C (2/4);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

GTTCGCCGGGTCCGCGTCGAACGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTGGTGCCTTTGCTTTCCAGGGT  >  NC_002947/1306775‑1307068
                                                                                                                                                     |                                                                                                                                                   
gTTCGCCGGGTCCGCGTCGAACGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑a                                                                                                                                                     >  6:398930/1‑149 (MQ=255)
   cGCCGGGTCCGCGTCGAACGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑At                                                                                                                                                    >  1:287640/1‑147 (MQ=255)
    gCCGGGTCCGCGTCGAACGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAg                                                                                                                                                 >  8:51068/1‑149 (MQ=255)
     ccGGGTCCGCGTCGAACGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGATCCTGCGCCAGG‑CTCAgg                                                                                                                                                <  3:285757/149‑1 (MQ=255)
     ccGGGTCCGCGTCGAACGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGATCCTGCGCCAGG‑CTCAgg                                                                                                                                                <  5:99490/149‑1 (MQ=255)
        ggTCCGCGTCGAACGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGc‑c                                                                                                                                             >  4:111313/1‑148 (MQ=255)
          tCCGCGTCGAACGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAg                                                                                                                                          <  2:332522/149‑1 (MQ=255)
           ccGCGTCGAACGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGc                                                                                                                                         >  3:84686/1‑149 (MQ=255)
               gTCGAACGACTTCGGCAGGGTGACCCGCAGGCCTTCTGCGGCGTTCGGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAgg                                                                                                                                                <  6:118357/139‑1 (MQ=255)
               gTCGAACGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTGTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAgg                                                                                                                                                >  5:118357/1‑139 (MQ=255)
                   aaCGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGggtg                                                                                                                                 >  3:68922/1‑149 (MQ=255)
                   aaCGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGggtg                                                                                                                                 <  8:18432/149‑1 (MQ=255)
                   aaCGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGggt                                                                                                                                  <  3:336078/148‑1 (MQ=255)
                   aaCGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGggt                                                                                                                                  >  4:336078/1‑148 (MQ=255)
                     cGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGggtggt                                                                                                                               <  5:27061/149‑1 (MQ=255)
                      gACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGggtgg                                                                                                                                >  2:432883/1‑147 (MQ=255)
                       aCTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAg                                                                                                                             >  3:139320/1‑149 (MQ=255)
                           cGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAgccg                                                                                                                          <  6:420792/148‑1 (MQ=255)
                                 ggTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCgtagt                                                                                                                   <  5:341905/149‑1 (MQ=255)
                                 ggTGACCAGCAGGCCTTCGGAGGCGTTCAGGGTGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCgtag                                                                                                                    >  1:248317/1‑148 (MQ=255)
                                     aCCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGAc                                                                                                               >  4:111434/1‑149 (MQ=255)
                                       cagcagGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACg                                                                                                              >  4:185210/1‑148 (MQ=255)
                                       cagcagGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACg                                                                                                              >  1:49493/1‑148 (MQ=255)
                                        agcagGCCTTCGGAGGCGTTCTGGGTGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGaa                                                                                                            <  2:248317/149‑1 (MQ=255)
                                        agcagGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACg                                                                                                              >  3:294089/1‑147 (MQ=255)
                                        agcagGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGaa                                                                                                            <  2:287640/149‑1 (MQ=255)
                                         tcagGCCTTCGGTGGAGGTCTGGTGGATCTGTGAAATGTAGTGGTTAGCCAGGGGTTCCGGACCCTTGAACTTGGCTTTCAGTCCTTTTTCCTGATCCTGGCGCAGT‑ATTAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAg                                                                                                           <  4:139320/148‑1 (MQ=255)
                                         gcagGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAg                                                                                                           <  1:351155/149‑1 (MQ=255)
                                           agGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTc                                                                                                         <  1:172041/149‑1 (MQ=255)
                                           agGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTc                                                                                                         <  2:49493/149‑1 (MQ=255)
                                               cTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACg                                                                                                     <  3:185210/149‑1 (MQ=255)
                                               cTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CAAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACg                                                                                                     >  6:49945/1‑149 (MQ=255)
                                                  cGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACgttg                                                                                                  >  8:278380/1‑149 (MQ=255)
                                                  cGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACgt                                                                                                    <  1:361320/147‑1 (MQ=255)
                                                  cGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACgt                                                                                                    <  4:68922/147‑1 (MQ=255)
                                                      ggCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑Ct                                                                                             <  8:356452/149‑1 (MQ=255)
                                                            cTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGCTTCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCtt                                                                                         <  1:107780/148‑1 (MQ=255)
                                                                      tGCGAAATGGAGTTGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGgat                                                                             >  2:181916/1‑149 (MQ=255)
                                                                        cGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCt                                                                                 <  1:432883/143‑1 (MQ=255)
                                                                        cGAAATCGAGTCGGTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGga                                                                              <  6:12711/146‑1 (MQ=255)
                                                                            aTGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCga                                                                       <  8:222060/149‑1 (MQ=255)
                                                                            aTGGAGTCGTTACCCACGCCTGCCGGACCGATGAACTTGGCTTTCATGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCga                                                                       <  7:350901/149‑1 (MQ=255)
                                                                             tGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCgag                                                                      <  4:84686/149‑1 (MQ=255)
                                                                               gcgTCGTTACCCACGCCTTCCTGACACATCAACATGGCTTTCAGTCCTTTTTCCTCAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTCGCCGCCGTAGTAGACCAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCgagga                                                                    <  8:132924/147‑1 (MQ=255)
                                                                               gAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCgag                                                                      >  3:387900/1‑147 (MQ=255)
                                                                                aGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCgaggag                                                                   >  2:141056/1‑149 (MQ=255)
                                                                                aGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCgag                                                                      >  7:319012/1‑146 (MQ=255)
                                                                                    gTTACCCACGGCTTCCGGACCCATGAACTTTGCTTTCACGCCTTTTTCCTGATGCTTGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTATTAGACGAAGTCCACGTTGTTCCTGCTTGAGCTTCTGGATGATCGAGGAGaa                                                                 <  1:260655/148‑1 (MQ=255)
                                                                                     ttACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAggc                                                              <  7:278380/149‑3 (MQ=255)
                                                                                      tACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTc                                                              <  1:181916/148‑1 (MQ=255)
                                                                                        cccACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCtt                                                           <  4:294089/149‑1 (MQ=255)
                                                                                        cccACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCtt                                                           <  7:51068/149‑1 (MQ=255)
                                                                                        cccACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCtt                                                           <  3:111313/149‑1 (MQ=255)
                                                                                           aCGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGt                                                         >  3:241565/1‑148 (MQ=255)
                                                                                              ccTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTGGACGTTGTT‑GTGCTTGGGCTGCGGGAGGATCGAGGAGATGTCCTTGTCAcg                                                     >  2:335886/1‑148 (MQ=255)
                                                                                               cTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCAccc                                                    >  6:249747/1‑149 (MQ=255)
                                                                                                 tCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCcgc                                                  >  4:11219/1‑149 (MQ=255)
                                                                                                  ccGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCcgcg                                                 >  5:121718/1‑149 (MQ=255)
                                                                                                          aTGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGcc                                         >  2:73902/1‑149 (MQ=255)
                                                                                                          aTGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGcc                                         >  7:190086/1‑149 (MQ=255)
                                                                                                          aTGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGcc                                         >  4:183285/1‑149 (MQ=255)
                                                                                                          aTGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGACGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCATTGTCACAAGCGTTCAGGca                                         >  3:237097/1‑148 (MQ=255)
                                                                                                           tGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCt                                        <  4:412211/149‑1 (MQ=255)
                                                                                                               cTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCt                                        >  4:162871/1‑145 (MQ=255)
                                                                                                               cTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCt                                        <  3:162871/145‑1 (MQ=255)
                                                                                                                ttGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGaa                                   >  5:334366/1‑149 (MQ=255)
                                                                                                                       tCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGc                             <  4:214866/148‑1 (MQ=255)
                                                                                                                       tCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGc                             >  4:255337/1‑148 (MQ=255)
                                                                                                                         aGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAAc                          >  1:417136/1‑149 (MQ=255)
                                                                                                                           gCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCt                        >  7:270419/1‑149 (MQ=255)
                                                                                                                            gcctgatCCTGCGCCAGGCTCAGGAAAAAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACg                               >  6:61514/7‑142 (MQ=255)
                                                                                                                            gcctgatCCTGCGCCAGGCTCAGGAAAAAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACg                               <  5:61514/136‑1 (MQ=255)
                                                                                                                              tttttCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTgg                     >  5:256811/1‑149 (MQ=255)
                                                                                                                              tttttCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTgg                     >  3:368072/1‑149 (MQ=255)
                                                                                                                               ttttCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTc                                              <  5:316481/123‑1 (MQ=255)
                                                                                                                               ttttCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTc                                              >  6:316481/1‑123 (MQ=255)
                                                                                                                                tttCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTGGt                    <  6:334366/148‑1 (MQ=255)
                                                                                                                                   ccTGATCCTGCGCCAGG‑CTCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCt                                        <  6:209395/125‑1 (MQ=255)
                                                                                                                                   ccTGATCCTGCGCCAGG‑CTCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCt                                        <  6:215231/125‑1 (MQ=255)
                                                                                                                                   ccTGATCCTGCGCCAGG‑CTCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCt                                        >  5:209395/1‑125 (MQ=255)
                                                                                                                                   ccTGATCCTGCGCCAGG‑CTCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCt                                        >  5:215231/1‑125 (MQ=255)
                                                                                                                                   ccTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTGGTGCCt                >  1:426184/1‑149 (MQ=255)
                                                                                                                                     tGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTCGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTACTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTGGTGCCttt              >  8:193925/1‑149 (MQ=255)
                                                                                                                                      gAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTGGTGCCTTTg             >  4:259318/1‑149 (MQ=255)
                                                                                                                                       aGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTGGTGCCTTTGc            <  6:256811/149‑1 (MQ=255)
                                                                                                                                               gcAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTGGTGCCTTTGCTTTCCAgg    <  8:190086/149‑1 (MQ=255)
                                                                                                                                                cAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTGGTGCCTTTGCTTTCCAggg   <  2:417136/149‑1 (MQ=255)
                                                                                                                                                 aGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTGGTGCCTTTGCTTTCCAGGGt  <  7:217742/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTTCGCCGGGTCCGCGTCGAACGACTTCGGCAGGGTGACCAGCAGGCCTTCGGAGGCGTTCTGGGCGATCTGCGAAATGGAGTCGTTACCCACGCCTTCCGGACCCATGAACTTGGCTTTCAGGCCTTTTTCCTGAGCCTGGCGCAGG‑ATCAGGC‑CCAGCTCTGGGTGGTAGCCGCCGTAGTAGACGAAGTCGACGTTGTT‑CTGCTTGAGCTTCTGGATGATCGAGGAGAAGTCCTTGTCACCCGCGTTCAGGCCTTCGAACACGGCAACCTTGGTGCCTTTGCTTTCCAGGGT  >  NC_002947/1306775‑1307068

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑