Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 964,481 T→A 6.5% L96Q (CTG→CAG) ‡ phnE → phosphonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947964,4810TA 6.5% 223.8 / 4.3 77L96Q (CTG→CAG) ‡phnEphosphonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (34/37);  new base A (3/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

CAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGCGCT  >  NC_002947/964333‑964629
                                                                                                                                                    |                                                                                                                                                    
cagATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCt                                                                                                                                                      <  8:94932/149‑1 (MQ=255)
cagATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCt                                                                                                                                                      >  6:304915/1‑149 (MQ=255)
  gATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTgg                                                                                                                                                    <  5:103871/149‑1 (MQ=255)
       ccAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCCAGCAAGCCAGc                                                                                                                                               <  8:30918/149‑1 (MQ=255)
         aGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCCAGCAAGCCAGCCg                                                                                                                                             <  5:91017/149‑1 (MQ=255)
           cTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTg                                                                                                                                           >  8:180376/1‑149 (MQ=255)
           cTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCAGGCCCTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCAGGCTTGCTGGCCAGCCGTg                                                                                                                                           >  7:263808/1‑149 (MQ=255)
              cACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCt                                                                                                                                        >  3:121403/1‑149 (MQ=255)
                ccAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGcc                                                                                                                                              <  6:320785/141‑1 (MQ=255)
                ccAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGcc                                                                                                                                              >  5:320785/1‑141 (MQ=255)
                 cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTg                                                                                                                                                     >  8:205228/1‑133 (MQ=255)
                 cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTg                                                                                                                                                     <  7:205228/133‑1 (MQ=255)
                 cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGt                                                                                                                                     >  7:3553/1‑149 (MQ=255)
                  aGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                    <  4:316310/149‑1 (MQ=255)
                        ttCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACACTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCg                                                                                                                              >  5:295235/1‑149 (MQ=255)
                             gccgccGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGtg                                                                                                                            <  6:345097/146‑1 (MQ=255)
                                  cgGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCa                                                                                                                    >  2:82404/1‑149 (MQ=255)
                                        atgatgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcggc                                                                                                              >  2:249864/1‑149 (MQ=255)
                                           atgCCGATTTTCTTGCACTGCTGTATGACGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTACTTCACGcggcgg                                                                                                            <  7:100235/148‑1 (MQ=255)
                                           atgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGcggcgg                                                                                                            <  7:160631/148‑1 (MQ=255)
                                              ccGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcgggcg                                                                                                         <  5:304915/148‑1 (MQ=255)
                                               cGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGcgggcggg                                                                                                       <  1:82404/149‑1 (MQ=255)
                                                gATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCgggcgggc                                                                                                      <  8:3553/149‑1 (MQ=255)
                                                 aTTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGcgggcggg                                                                                                       >  8:135078/1‑147 (MQ=255)
                                                 aTTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGcgggcggg                                                                                                       <  7:135078/147‑1 (MQ=255)
                                                 aTTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggcc                                                                                                     <  1:249864/149‑1 (MQ=255)
                                                  ttttCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggccg                                                                                                    <  4:121403/149‑1 (MQ=255)
                                                             tgctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCAGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTcg                                                                                         >  6:242649/1‑149 (MQ=255)
                                                                    tGAAGCGACCTTGCAGACCCAGGCCGTGCCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTg                                                                                   <  7:74771/148‑1 (MQ=255)
                                                                     gAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGtt                                                                                 >  1:232491/1‑149 (MQ=255)
                                                                      aaGCTAGCTTGCTGACGCTGGCCTTGGCTTCTGCCGTCATGGCGCTGTCGTTGTTTTTGGCGATCCCGGCTTGCTTGCTGGCCAGCCGTGGCTTGTCCTTGCGTGCTGCTTCACGCGGCGGTCGGGCCGGGCTCTGGTCGCGACTGTTg                                                                                <  2:357786/149‑1 (MQ=255)
                                                                       aGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTgc                                                                               >  1:337518/1‑149 (MQ=255)
                                                                           aCCTTGCAGACTCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGa                                                                                      >  3:144064/1‑138 (MQ=255)
                                                                           aCCTTGCAGACACTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGa                                                                                      <  4:144064/138‑1 (MQ=255)
                                                                             cTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGc                                                                         <  2:177241/149‑1 (MQ=255)
                                                                             cTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGc                                                                         >  7:64541/1‑149 (MQ=255)
                                                                             cTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGc                                                                         <  7:180376/149‑1 (MQ=255)
                                                                              ttGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGcc                                                                        >  4:89621/1‑149 (MQ=255)
                                                                                gCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCg                                                                       <  8:64541/148‑1 (MQ=255)
                                                                                  aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcggc                                                                                                              <  8:26023/107‑1 (MQ=255)
                                                                                  aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcggc                                                                                                              >  7:26023/1‑107 (MQ=255)
                                                                                      ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTg                                                                >  4:67444/1‑149 (MQ=255)
                                                                                      ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTg                                                                >  7:16133/1‑149 (MQ=255)
                                                                                      ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTg                                                                >  7:21721/1‑149 (MQ=255)
                                                                                      ccTGGCCGTGGCGACTGCCGGCAAGGCGCTGGCGTTGTTGTTGGCGATCCCGTCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGGTGCGCCTGCCGGTGCGTg                                                                >  7:272400/1‑149 (MQ=255)
                                                                                        tGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTggg                                                              >  5:272888/1‑149 (MQ=255)
                                                                                         ggCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGc                                                             <  5:240993/149‑1 (MQ=255)
                                                                                           ccGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTGGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGTCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTa                                                           <  6:295235/149‑1 (MQ=255)
                                                                                             gTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACt                                                         <  3:347516/149‑1 (MQ=255)
                                                                                               ggCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGa                                                       <  6:272888/149‑1 (MQ=255)
                                                                                                  gACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGAttt                                                    >  2:118034/1‑149 (MQ=255)
                                                                                                    cTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTc                                                  >  2:192827/1‑149 (MQ=255)
                                                                                                             tGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGtgctg                                                                                                                         <  4:89945/69‑1 (MQ=255)
                                                                                                             tGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGtgctg                                                                                                                         >  3:89945/1‑69 (MQ=255)
                                                                                                              ggCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGAGTGGGCTACTGATTTTCCTGCGTAGc                                         <  1:5792/148‑1 (MQ=255)
                                                                                                                 ggtGGGGTTGTGGTTGGCGATCCCGGCTGGCTTGCGGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGGGGGCGGGGCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGc                                     <  8:263808/147‑1 (MQ=255)
                                                                                                                  cTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGcc                                    <  3:89621/149‑1 (MQ=255)
                                                                                                                    ggCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCgg                                  <  5:242649/149‑1 (MQ=255)
                                                                                                                    ggCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCg                                   >  1:266941/1‑148 (MQ=255)
                                                                                                                    ggCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCg                                   <  2:266941/148‑1 (MQ=255)
                                                                                                                     gCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCgg                                  >  3:294146/1‑148 (MQ=255)
                                                                                                                           ttgttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGt                           <  1:192827/149‑1 (MQ=255)
                                                                                                                              ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCCCGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTg                        <  2:232491/149‑1 (MQ=255)
                                                                                                                                   gATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGTCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGGAGGAAATCGTCTGGGCGt                   >  2:31191/1‑149 (MQ=255)
                                                                                                                                    aTCCAGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCTGGCGGGCCTGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTGCAGCGTAGCTTGCCGGAAATCGTCTGGGCGtt                  <  1:31191/149‑1 (MQ=255)
                                                                                                                                     tCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTg                 >  3:202473/1‑149 (MQ=255)
                                                                                                                                     gcacGGCTGGCCAGCAAGCCAGCCGGGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTg                 >  3:368414/4‑149 (MQ=255)
                                                                                                                                      cccGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTg                 <  2:337518/148‑1 (MQ=255)
                                                                                                                                      cacGGCTGGCCAGCAAGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGc                >  7:211630/3‑149 (MQ=255)
                                                                                                                                       ccGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGc                                     <  3:229469/127‑1 (MQ=255)
                                                                                                                                       ccGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGc                                     >  4:229469/1‑127 (MQ=255)
                                                                                                                                       ccGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCt               >  3:291812/1‑149 (MQ=255)
                                                                                                                                       acGGCTGGCCAGCAAGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCAAGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTACTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCt               >  6:39789/2‑149 (MQ=255)
                                                                                                                                        cGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGtt                                                                                 >  6:33499/1‑82 (MQ=255)
                                                                                                                                        cGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGtt                                                                                 <  5:33499/82‑1 (MQ=255)
                                                                                                                                        cGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTg              >  6:147501/1‑149 (MQ=255)
                                                                                                                                           ctggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGNTGCTGTTc           <  8:21721/149‑1 (MQ=255)
                                                                                                                                           ctggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGGGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTc           <  8:16133/149‑1 (MQ=255)
                                                                                                                                                 tGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTgcgc     >  5:109849/1‑149 (MQ=255)
                                                                                                                                                    tGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGaa                                >  6:75507/1‑119 (MQ=255)
                                                                                                                                                    tGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGCGCt  >  8:330009/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGCGCT  >  NC_002947/964333‑964629

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: