Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,003,644 C→A 5.5% G98V (GGC→GTC) ‡ amaB ← L‑piperidine‑6‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,003,6440CA 5.5% 295.6 / 5.4 90G98V (GGC→GTC) ‡amaBL‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (46/39);  new base A (2/3);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

CTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGC  >  NC_002947/6003497‑6003789
                                                                                                                                                    |                                                                                                                                                 
cTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGTCCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGc                                                                                                                                                   <  1:343431/149‑1 (MQ=255)
       cGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                           >  6:302883/1‑149 (MQ=255)
       cGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                           <  4:212882/149‑1 (MQ=255)
        gCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTggg                                                                                                                                          >  2:420486/1‑149 (MQ=255)
        gCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTggg                                                                                                                                          <  4:411868/149‑1 (MQ=255)
            aCCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTGGGTGaa                                                                                                                                      <  8:360443/149‑2 (MQ=255)
             ccACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTGGGTGaag                                                                                                                                     <  8:257157/149‑3 (MQ=255)
              cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGt                                                                                                                                         <  1:215818/144‑1 (MQ=255)
              cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGt                                                                                                                                         >  2:215818/1‑144 (MQ=255)
              cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCt                                                                                                                                    >  7:378764/1‑149 (MQ=255)
              cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCt                                                                                                                                    >  7:97620/1‑149 (MQ=255)
              cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACATCACCCAGGAAGGCCTGGGTGaagt                                                                                                                                    <  6:425896/149‑4 (MQ=255)
                   ccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCg                                                                                                                               <  8:81703/149‑1 (MQ=255)
                   ccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCg                                                                                                                               >  8:252010/1‑149 (MQ=255)
                   ccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCg                                                                                                                               >  1:310291/1‑149 (MQ=255)
                   ccAGCGGGTGCCAGGTTGAACGCATGTGGTGGCCCGGGCGCTCGGCAGCGAGGGTCAGGCCGTACCGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCg                                                                                                                               <  4:17951/149‑1 (MQ=255)
                       cGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCtt                                                                                                                           <  2:424933/149‑1 (MQ=255)
                            gCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatc                                                                                                                      >  2:344671/1‑149 (MQ=255)
                              cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCATACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                                                                                                                    <  7:371505/149‑1 (MQ=255)
                                 gTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGAc                                                                                                                 <  4:124867/149‑1 (MQ=255)
                                 gTTTCACGCATGTGGTGGCCCGGGCACTCGGAAGCGATGGTCAGGCCGTAAAGCAGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACACAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgaga                                                                                                                  >  3:175675/1‑148 (MQ=255)
                                      aCGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCt                                                                                                            >  1:431361/1‑149 (MQ=255)
                                      aCGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCt                                                                                                            >  6:50806/1‑149 (MQ=255)
                                      aCGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCt                                                                                                            >  3:146036/1‑149 (MQ=255)
                                      aCGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCt                                                                                                            >  1:281081/1‑149 (MQ=255)
                                        gCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCg                                                                                                          <  1:344671/149‑1 (MQ=255)
                                           tgtgGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTg                                                                                                                                        >  4:409201/1‑116 (MQ=255)
                                           tgtgGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTg                                                                                                                                        <  3:409201/116‑1 (MQ=255)
                                           tgtgGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCg                                                                                                       <  6:129702/149‑1 (MQ=255)
                                             tggtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAg                                                                                                     >  7:173287/1‑149 (MQ=255)
                                             tggtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAACTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAg                                                                                                     <  4:369980/149‑1 (MQ=255)
                                                 ggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCg                                                                                                 <  2:431361/149‑1 (MQ=255)
                                                  gCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCg                                                                                                 <  6:108929/148‑1 (MQ=255)
                                                  gCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCg                                                                                                 <  7:22826/148‑1 (MQ=255)
                                                   cccGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGc                                                                                               <  4:11036/149‑1 (MQ=255)
                                                   cccGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGc                                                                                               <  2:399096/149‑1 (MQ=255)
                                                    ccGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGcc                                                                                              <  1:420486/149‑1 (MQ=255)
                                                      gggCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCtt                                                                                            <  5:302883/149‑1 (MQ=255)
                                                      gggCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCtt                                                                                            >  7:386189/1‑149 (MQ=255)
                                                      gggCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCtt                                                                                            <  7:252010/149‑1 (MQ=255)
                                                         cgcTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtgtg                                                                                         <  3:40172/149‑1 (MQ=255)
                                                          gcTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGc                                                                                        >  1:111813/1‑149 (MQ=255)
                                                          gcTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGc                                                                                        <  3:152619/149‑1 (MQ=255)
                                                           cTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCt                                                                                       <  2:310291/149‑1 (MQ=255)
                                                            tCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCt                                                                                       >  2:3100/1‑148 (MQ=255)
                                                            tCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCt                                                                                       >  8:275874/1‑148 (MQ=255)
                                                              ggAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCTACGGCGAAGTCTCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCAc                                                                                    <  7:148981/149‑1 (MQ=255)
                                                               gAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTc                                                                                                  <  8:131789/134‑1 (MQ=255)
                                                               gAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTc                                                                                                  >  7:131789/1‑134 (MQ=255)
                                                                        gTCAGGCTGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACtt                                                                           <  2:20164/148‑1 (MQ=255)
                                                                          cAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGa                                                                                                      <  4:241960/119‑1 (MQ=255)
                                                                          cAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGa                                                                                                      >  3:241960/1‑119 (MQ=255)
                                                                           aGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGcc                                                                       <  8:386189/149‑1 (MQ=255)
                                                                           aGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGcc                                                                       >  6:275182/1‑149 (MQ=255)
                                                                               cGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACCTCGCCGAAc                                                                   <  2:281081/149‑1 (MQ=255)
                                                                                  aCAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAgg                                                                >  1:144922/1‑149 (MQ=255)
                                                                                    aGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGgcg                                                              >  1:371507/1‑149 (MQ=255)
                                                                                          gcgcGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCa                                                         <  4:175675/149‑1 (MQ=255)
                                                                                             cgACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGtt                                                      <  8:378764/149‑1 (MQ=255)
                                                                                              gACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACa                                                                  >  1:175199/1‑136 (MQ=255)
                                                                                                cAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGc                                                   >  5:349835/1‑149 (MQ=255)
                                                                                                   gCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACg                                                                                   <  5:278429/114‑1 (MQ=255)
                                                                                                   gCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACg                                                                                   >  6:278429/1‑114 (MQ=255)
                                                                                                   gCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCAc                                                >  2:106132/1‑149 (MQ=255)
                                                                                                   gCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCATGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCAc                                                >  1:21706/1‑149 (MQ=255)
                                                                                                     cAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACgg                                              >  5:266829/1‑149 (MQ=255)
                                                                                                      aaCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgcagc                                                                               <  3:51700/115‑1 (MQ=255)
                                                                                                      aaCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgcagc                                                                               >  4:51700/1‑115 (MQ=255)
                                                                                                       aCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTc                                                     >  2:318356/1‑140 (MQ=255)
                                                                                                       aCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTc                                                     <  1:318356/140‑1 (MQ=255)
                                                                                                           cGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGcgcg                                          >  4:211015/1‑147 (MQ=255)
                                                                                                                   cAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACAAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGca                                 >  5:258664/1‑148 (MQ=255)
                                                                                                                     gATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCacac                              >  5:320411/1‑149 (MQ=255)
                                                                                                                            aTCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTgcgc                         >  5:396147/1‑147 (MQ=255)
                                                                                                                             tCATTTCCTGTACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAg                      >  1:409728/1‑149 (MQ=255)
                                                                                                                               aTTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAgg                     <  7:275874/148‑1 (MQ=255)
                                                                                                                                 ttCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCtt                  >  5:85019/1‑149 (MQ=255)
                                                                                                                                  tCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCtt                  <  2:144922/148‑1 (MQ=255)
                                                                                                                                  tCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCtt                  <  6:349835/148‑1 (MQ=255)
                                                                                                                                       cACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAgg            >  8:102982/1‑149 (MQ=255)
                                                                                                                                       cACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAg             <  4:146036/148‑1 (MQ=255)
                                                                                                                                         gatCACCCAGGAAGGCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGgcg          >  8:100422/3‑149 (MQ=255)
                                                                                                                                         gatCACCCAGGAAGGCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGgcg          >  5:418599/3‑149 (MQ=255)
                                                                                                                                         cTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGgcg          >  8:300869/1‑149 (MQ=255)
                                                                                                                                             aCCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGctc       <  3:211015/148‑1 (MQ=255)
                                                                                                                                              cccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACAATGCGCCAGGCTTCGAAGGCGCTCTg     >  2:415444/1‑149 (MQ=255)
                                                                                                                                               ccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTgg    >  8:361093/1‑149 (MQ=255)
                                                                                                                                               ccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTgg    >  7:168583/1‑149 (MQ=255)
                                                                                                                                               ccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTgg    >  3:394552/1‑149 (MQ=255)
                                                                                                                                                 aGGCCTTCCTGGGTGATCTTGTCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGc  >  1:224311/1‑149 (MQ=255)
                                                                                                                                                 aGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGc  >  6:156081/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGC  >  NC_002947/6003497‑6003789

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: