Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,734,772 G→C 7.1% R298G (CGT→GGT)  hutH ← histidine ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,734,7720GC 7.1% 276.3 / 9.8 84R298G (CGT→GGT) hutHhistidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (40/38);  new base C (4/2);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GCGATGGCCAGGGCCAGGTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATC  >  NC_002947/5734624‑5734914
                                                                                                                                                    |                                                                                                                                              
gCGATGGCCAGGGCCAGGTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACg                                                                                                                                                >  8:8413/1‑149 (MQ=255)
  gATGGCCAGGGCCAGGTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCa                                                                                                                                              <  7:165032/149‑1 (MQ=255)
    tGGCCAGGGCCAGGTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGc                                                                                                                                            <  8:68000/149‑1 (MQ=255)
      gccagggccaggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTg                                                                                                                                          >  4:322759/1‑149 (MQ=255)
      gccagggccaggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTg                                                                                                                                          >  2:246312/1‑149 (MQ=255)
       ccagggccaggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTgg                                                                                                                                         >  1:130018/1‑149 (MQ=255)
            gccaggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTTACGCAGCTGGGTcag                                                                                                                                    <  2:96324/149‑1 (MQ=255)
            gccaggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcag                                                                                                                                    <  2:293077/149‑1 (MQ=255)
               aggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcaggca                                                                                                                                 <  5:152223/149‑1 (MQ=255)
               aggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcaggca                                                                                                                                 <  6:368980/149‑1 (MQ=255)
                 gTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcaggcagg                                                                                                                               <  6:27465/149‑1 (MQ=255)
                  ttGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCaggcaggc                                                                                                                              <  5:128062/149‑1 (MQ=255)
                    accggtggCCATGGCCGCCGGTTCGGCGTGGAAGTTGCCGCCGGGGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCAc                                                                                                                            >  4:153921/3‑149 (MQ=255)
                              aTGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTg                                                                                                                  <  2:130018/149‑1 (MQ=255)
                              aTGGCCACCGGTTCCGCGTGGAAGTCGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTCTCCGCCACAGCGTTGGCTTCCATGACCAGTACCTCGGCGGCCTGACGAAGCTGCGTGATACAGGCACCCaagcccgc                                                                                                                  >  6:285201/1‑142 (MQ=255)
                                 gCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCgg                                                                                                               <  3:307391/149‑1 (MQ=255)
                                 gCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCgg                                                                                                               <  2:133655/149‑1 (MQ=255)
                                  ccACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGc                                                                                                              <  3:322759/149‑1 (MQ=255)
                                  ccACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGc                                                                                                              >  2:395271/1‑149 (MQ=255)
                                   cACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCTAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGGCCAGTACCTCGGCGGCCTGACGCATCTGGGTCAGGCATGCACCCATGACCTGCGGCt                                                                                                             <  4:232633/149‑1 (MQ=255)
                                   cACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCt                                                                                                             >  7:129899/1‑149 (MQ=255)
                                     ccGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTgg                                                                                                           <  7:9096/149‑1 (MQ=255)
                                     ccGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTgg                                                                                                           <  1:395271/149‑1 (MQ=255)
                                           cGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcgc                                                                                                     >  7:31615/1‑149 (MQ=255)
                                            ggCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcgc                                                                                                     <  7:354516/148‑1 (MQ=255)
                                            ggCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcg                                                                                                      >  6:57321/1‑147 (MQ=255)
                                             gCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcgcag                                                                                                   <  7:8413/149‑1 (MQ=255)
                                                tGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCg                                                                                                                >  8:303250/1‑133 (MQ=255)
                                                  gAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAg                                                                                               >  4:185116/1‑148 (MQ=255)
                                                  gAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAg                                                                                               <  2:324657/148‑1 (MQ=255)
                                                      ttGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGTCCTGCGGCTGGCAGCGCAGCGAGTACgg                                                                                          <  5:334465/149‑1 (MQ=255)
                                                          cgccgGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTcc                                                                                      >  5:47896/1‑149 (MQ=255)
                                                              ggAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGa                                                                                  >  5:229262/1‑149 (MQ=255)
                                                               gagaTCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGAc                                                                                 >  6:245360/1‑149 (MQ=255)
                                                                agaTCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGAcc                                                                                >  3:128046/1‑149 (MQ=255)
                                                                    cacatcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGt                                                                            >  5:406856/1‑149 (MQ=255)
                                                                     acatcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTc                                                                           >  2:30253/1‑149 (MQ=255)
                                                                      catcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCt                                                                                     <  5:218489/138‑1 (MQ=255)
                                                                      catcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCt                                                                                     >  6:218489/1‑138 (MQ=255)
                                                                       atcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCg                                                                          >  2:118211/1‑148 (MQ=255)
                                                                        tcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCg                                                                          >  6:128751/1‑147 (MQ=255)
                                                                           ccTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGAGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGAAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGt                                                                      >  2:315464/1‑148 (MQ=255)
                                                                                  cAGCGAACACCGGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgc                                                              <  2:232835/149‑1 (MQ=255)
                                                                                   aGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcg                                                             <  5:245360/149‑1 (MQ=255)
                                                                                       aaCACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcgacag                                                         >  4:239191/1‑149 (MQ=255)
                                                                                            cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcgacagcgaca                                                    >  4:145263/1‑149 (MQ=255)
                                                                                            cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcgacagcgaca                                                    >  3:272371/1‑149 (MQ=255)
                                                                                            cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcgacagcgac                                                     >  7:82726/1‑148 (MQ=255)
                                                                                              gCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACt                                                  <  1:30253/149‑1 (MQ=255)
                                                                                                   ttGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCt                                             <  6:47896/149‑1 (MQ=255)
                                                                                                     gTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                            <  1:93617/148‑1 (MQ=255)
                                                                                                      tCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGa                                          <  5:128751/149‑1 (MQ=255)
                                                                                                        cGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCAGCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGt                                        >  5:117408/1‑149 (MQ=255)
                                                                                                           cacaGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGc                                     <  5:52236/149‑1 (MQ=255)
                                                                                                             caGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGccc                                   <  6:229262/149‑1 (MQ=255)
                                                                                                                cGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                            <  5:172909/137‑1 (MQ=255)
                                                                                                                cGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                            >  6:172909/1‑137 (MQ=255)
                                                                                                                cGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcagc                                <  8:82726/149‑1 (MQ=255)
                                                                                                                cGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcag                                 >  1:183766/1‑148 (MQ=255)
                                                                                                                 gTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                            <  1:37736/136‑1 (MQ=255)
                                                                                                                 gTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                            >  2:37736/1‑136 (MQ=255)
                                                                                                                 gTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcag                                 >  1:320588/1‑147 (MQ=255)
                                                                                                                  ttGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcagcag                              >  1:275689/1‑149 (MQ=255)
                                                                                                                     gCTTCGATGCCCAGTACCTCGGCGGCCTGGCGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTc                           >  8:296498/1‑149 (MQ=255)
                                                                                                                          gATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCgg                       >  4:158681/1‑148 (MQ=255)
                                                                                                                            tGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAg                    <  8:129899/149‑1 (MQ=255)
                                                                                                                                cAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAgg                <  4:128046/149‑1 (MQ=255)
                                                                                                                                cAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAgg                >  3:55700/1‑149 (MQ=255)
                                                                                                                                     ccTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCcgcg           <  5:57321/149‑1 (MQ=255)
                                                                                                                                      cTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAg                                                                       <  8:270752/88‑1 (MQ=255)
                                                                                                                                      cTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAg                                                                       >  7:270752/1‑88 (MQ=255)
                                                                                                                                      cTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCgg          >  6:153488/1‑149 (MQ=255)
                                                                                                                                       tCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGgt         >  6:138687/1‑149 (MQ=255)
                                                                                                                                       tCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGgt         >  6:247103/1‑149 (MQ=255)
                                                                                                                                         ggcggcCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGgt         <  3:239191/147‑1 (MQ=255)
                                                                                                                                           ctgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGACAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTcgc                         >  3:257769/3‑129 (MQ=255)
                                                                                                                                           ctgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGACAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTcgc                         <  4:257769/127‑1 (MQ=255)
                                                                                                                                            tgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGa                      <  6:97340/130‑1 (MQ=255)
                                                                                                                                            tgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGa                      >  5:97340/2‑131 (MQ=255)
                                                                                                                                            tgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGa    >  6:259819/2‑149 (MQ=255)
                                                                                                                                            tgcCTGACCCAGCTGCGTCAGGCAGGAACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCg     >  3:402896/2‑148 (MQ=255)
                                                                                                                                            ggcCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGa    <  3:158681/149‑1 (MQ=255)
                                                                                                                                            ggcCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGa    <  3:153921/149‑1 (MQ=255)
                                                                                                                                              ccTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGGAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATc  >  1:128977/1‑149 (MQ=255)
                                                                                                                                              ccTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATc  >  1:274969/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GCGATGGCCAGGGCCAGGTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATC  >  NC_002947/5734624‑5734914

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: ‑