Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 964,478 T→A 6.1% L95* (TTG→TAG) ‡ phnE → phosphonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947964,4780TA 6.1% 242.3 / 5.1 82L95* (TTG→TAG) ‡phnEphosphonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (37/39);  new base A (3/2);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

GCAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGC  >  NC_002947/964332‑964626
                                                                                                                                                  |                                                                                                                                                    
gcagATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGc                                                                                                                                                    >  1:177241/1‑149 (MQ=255)
gcagATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGc                                                                                                                                                    >  6:240993/1‑149 (MQ=255)
 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                                                                                                                                                  <  7:205228/133‑1 (MQ=255)
                  cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTg                                                                                                                                                  >  8:205228/1‑133 (MQ=255)
                  cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGt                                                                                                                                  >  7:3553/1‑149 (MQ=255)
                   aGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                 <  4:316310/149‑1 (MQ=255)
                    gcgcCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGc                                                                                                                                                    <  7:232049/129‑1 (MQ=255)
                    gcgcCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGc                                                                                                                                                    >  8:232049/1‑129 (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                                                                      <  7:180376/149‑1 (MQ=255)
                                                                              cTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGc                                                                      >  7:64541/1‑149 (MQ=255)
                                                                              cTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGc                                                                      <  2:177241/149‑1 (MQ=255)
                                                                               ttGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGcc                                                                     >  4:89621/1‑149 (MQ=255)
                                                                                 gCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCg                                                                    <  8:64541/148‑1 (MQ=255)
                                                                                   aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcggc                                                                                                           >  7:26023/1‑107 (MQ=255)
                                                                                   aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcggc                                                                                                           <  8:26023/107‑1 (MQ=255)
                                                                                       ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTg                                                             >  7:21721/1‑149 (MQ=255)
                                                                                       ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTg                                                             >  4:67444/1‑149 (MQ=255)
                                                                                       ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTg                                                             >  7:16133/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                                                                                                                      >  3:89945/1‑69 (MQ=255)
                                                                                                              tGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGtgctg                                                                                                                      <  4:89945/69‑1 (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                                                                              <  5:33499/82‑1 (MQ=255)
                                                                                                                                         cGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGtt                                                                              >  6:33499/1‑82 (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)
                                                                                                                                                  |                                                                                                                                                    
GCAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGC  >  NC_002947/964332‑964626

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: