Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 964,481 T→A 6.2% 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,4810TA 6.2% 1117.3 / 42.1 331L95* (TTG→TAG) ‡phnEphosphonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (187/122);  new base A (14/7);  total (202/129)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 338 total aligned reads displayed.

TCAGCAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGC  >  NC_002947/964332‑964629
                                                                                                                                                     |                                                                                                                                                    
tCAGCAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCtt                                                                                                                                                      >  1:290739/1‑150 (MQ=255)
  agcagATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGTTTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGc                                                                                                                                                    >  1:6010158/1‑150 (MQ=255)
   gcagATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCt                                                                                                                                                   <  1:2805090/150‑1 (MQ=255)
    cagATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTg                                                                                                                                                  <  1:2545851/150‑1 (MQ=255)
        ttGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCa                                                                                                                                              <  2:4189241/150‑1 (MQ=255)
           ccAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGcc                                                                                                                                           >  1:1295071/1‑150 (MQ=255)
            cAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCg                                                                                                                                          >  2:4631939/1‑150 (MQ=255)
              gCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCCAGCAAGCCAGCCggg                                                                                                                                        <  1:1455625/150‑3 (MQ=255)
               cTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGc                                                                                                                                       >  1:905719/1‑150 (MQ=255)
               cTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCACGGCTGGCCAGCAAGCCAGCCggga                                                                                                                                       <  1:1650350/150‑4 (MQ=255)
                 tCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCACGGCTGGCCAGCAAGCCAGCCgggatc                                                                                                                                     <  1:2699699/150‑6 (MQ=255)
                     ccgcgcCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGTCCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                 >  2:2489019/3‑150 (MQ=255)
                     cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                 >  2:4785275/1‑150 (MQ=255)
                     cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                 >  1:998310/1‑150 (MQ=255)
                     cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                 <  1:3551739/150‑1 (MQ=255)
                     cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                 >  2:2111231/1‑150 (MQ=255)
                        cgcCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCtt                                                                                                                              <  1:4987605/150‑1 (MQ=255)
                         gcCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTg                                                                                                                             >  2:1710939/1‑150 (MQ=255)
                            ttCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGt                                                                                                                          >  2:244383/1‑150 (MQ=255)
                             acTGGCCGCCGGCACAAGATGCCGATTTTCTTGCACTGCAGAATGAAGCTACCTTGCAGACCCTGCCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCCGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCgag                                                                                                                         <  2:6425807/149‑3 (MQ=255)
                               tGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGtgct                                                                                                                       >  1:1509740/1‑150 (MQ=255)
                                ggCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGtgctg                                                                                                                      <  2:6183418/150‑1 (MQ=255)
                                  ccgccgGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTgctgct                                                                                                                    >  1:5547619/1‑150 (MQ=255)
                                    gccgGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTc                                                                                                                  >  1:1064637/1‑150 (MQ=255)
                                       ggCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCAcg                                                                                                               >  2:6674508/1‑150 (MQ=255)
                                        gCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCAcgc                                                                                                              >  2:567074/1‑150 (MQ=255)
                                            atgatgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcggcg                                                                                                          >  2:1926425/1‑150 (MQ=255)
                                              gatgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcggg                                                                                                        <  2:2630700/150‑1 (MQ=255)
                                               atgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcgggc                                                                                                       <  1:827449/150‑1 (MQ=255)
                                                tgCCGATTTTCTTGCACTGATGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcgggcg                                                                                                      >  1:472637/1‑150 (MQ=255)
                                                    gATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggcc                                                                                                  >  1:6773012/1‑150 (MQ=255)
                                                    gATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggcc                                                                                                  >  1:6401697/1‑150 (MQ=255)
                                                    gATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggcc                                                                                                  >  2:1360942/1‑150 (MQ=255)
                                                       tttCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggccggg                                                                                               <  2:4909728/150‑1 (MQ=255)
                                                           ttGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGctct                                                                                           >  1:5938149/1‑150 (MQ=255)
                                                           ttGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGctct                                                                                           <  2:3698581/150‑1 (MQ=255)
                                                              cACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGACCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGt                                                                                        >  2:5591639/1‑150 (MQ=255)
                                                              cACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGGTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGt                                                                                        >  2:639259/1‑150 (MQ=255)
                                                                 tgctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcgggcg                                                                                                      <  2:3587917/133‑1 (MQ=255)
                                                                 tgctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTcgc                                                                                     <  2:5347500/150‑1 (MQ=255)
                                                                   ctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGa                                                                                   >  1:6992887/1‑150 (MQ=255)
                                                                      tATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTg                                                                                >  2:2993108/1‑150 (MQ=255)
                                                                      tATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTg                                                                                >  2:2760212/1‑150 (MQ=255)
                                                                        tGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGtt                                                                              >  1:7227916/1‑150 (MQ=255)
                                                                         gAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTg                                                                             >  1:4335395/1‑150 (MQ=255)
                                                                          aaGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTgc                                                                            <  1:5396315/150‑1 (MQ=255)
                                                                            gCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTgcgc                                                                          <  2:4008600/150‑1 (MQ=255)
                                                                                ccTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGc                                                                      >  2:324882/1‑150 (MQ=255)
                                                                                  ttGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCg                                                                    >  2:6564228/1‑150 (MQ=255)
                                                                                    gCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGt                                                                  >  1:4648334/1‑150 (MQ=255)
                                                                                    gCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGt                                                                  >  2:1463116/1‑150 (MQ=255)
                                                                                    gCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGt                                                                  <  2:4444161/150‑1 (MQ=255)
                                                                                      aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCtt                                                                                                                   <  1:2965928/99‑1 (MQ=255)
                                                                                      aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCtt                                                                                                                   >  2:2965928/1‑99 (MQ=255)
                                                                                      aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGc                                                                >  1:424889/1‑150 (MQ=255)
                                                                                      aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGc                                                                >  2:5928126/1‑150 (MQ=255)
                                                                                          ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTgg                                                            >  1:1828601/1‑150 (MQ=255)
                                                                                          ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTgg                                                            >  1:5512157/1‑150 (MQ=255)
                                                                                          ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTgg                                                            >  2:3532648/1‑150 (MQ=255)
                                                                                          ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCGCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTgg                                                            >  2:5407849/1‑150 (MQ=255)
                                                                                              gCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTa                                                        <  2:1862107/150‑1 (MQ=255)
                                                                                              gCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTa                                                        >  2:5857559/1‑150 (MQ=255)
                                                                                               ccGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTAc                                                       >  1:5884576/1‑150 (MQ=255)
                                                                                                      gacagacGGCATGCCGCTGGCGTTGTTGTTGGCGATCCCTGCTCGATAGCTTGCCAGCCGTGCATTGTCCTTGCGTGCTGCGTCACGCGGAGGGCGGGCCGGGCTAGCGTCGCCACTGTTGCGCCTTCAGGAGCGTGGTCAACTTAtttt                                                >  2:4973919/7‑150 (MQ=255)
                                                                                                       aCTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTc                                               <  2:6146644/150‑1 (MQ=255)
                                                                                                        cTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTcc                                              <  2:7437132/150‑1 (MQ=255)
                                                                                                         tGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGTCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGTCTACTGATTTTCCt                                             <  1:4565206/150‑1 (MQ=255)
                                                                                                              gCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTa                                        >  1:6180237/1‑150 (MQ=255)
                                                                                                              gCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGACCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTACTGCGTa                                        >  1:239812/1‑150 (MQ=255)
                                                                                                               cATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAg                                       <  1:4196231/150‑1 (MQ=255)
                                                                                                                aTGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGc                                      <  2:6072179/150‑1 (MQ=255)
                                                                                                                aTGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGc                                      <  2:4792625/150‑1 (MQ=255)
                                                                                                                 tGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTTCGTAGCg                                     >  1:5720914/1‑150 (MQ=255)
                                                                                                                  ggCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGt                                    <  1:5808996/150‑1 (MQ=255)
                                                                                                                  ggCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGt                                    <  1:4076261/150‑1 (MQ=255)
                                                                                                                  ggCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGt                                    <  2:7124312/150‑1 (MQ=255)
                                                                                                                     gcTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGcc                                 >  1:2205452/1‑150 (MQ=255)
                                                                                                                     gcTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGcc                                 >  1:4013598/1‑150 (MQ=255)
                                                                                                                       tGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCgg                               >  2:395114/1‑150 (MQ=255)
                                                                                                                        ggCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGa                              >  2:1819747/1‑150 (MQ=255)
                                                                                                                          cGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGaaa                            <  1:6702034/150‑1 (MQ=255)
                                                                                                                          cGTTGTTGTTGGCGAGCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGaaa                            <  2:7033835/150‑1 (MQ=255)
                                                                                                                           gttgttgttGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAAt                           <  1:2231406/150‑1 (MQ=255)
                                                                                                                            ttgttgttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATc                          <  1:753208/150‑1 (MQ=255)
                                                                                                                              gttgttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGt                        >  2:7226889/1‑150 (MQ=255)
                                                                                                                                 gttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGGGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTg                     >  2:5228908/1‑150 (MQ=255)
                                                                                                                                  ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                    >  1:1408223/1‑150 (MQ=255)
                                                                                                                                  ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                    <  1:717274/150‑1 (MQ=255)
                                                                                                                                  ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                    <  1:4890958/150‑1 (MQ=255)
                                                                                                                                  ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCGCTGGGCGCGGCTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                    >  1:1001430/1‑150 (MQ=255)
                                                                                                                                        aTCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAACTCGTCTGGGCGTTg              >  2:4588858/1‑150 (MQ=255)
                                                                                                                                         gcccGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGc             >  2:3748912/2‑150 (MQ=255)
                                                                                                                                         gcacGGCTGGCCAGCAAGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGc             >  1:3510328/4‑150 (MQ=255)
                                                                                                                                         gcacGGCTGGCCAGCAAGCCAGCCGGGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGc             >  1:2447305/4‑150 (MQ=255)
                                                                                                                                          cccGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCt            >  2:3086044/1‑150 (MQ=255)
                                                                                                                                           acGGCTGGCCAGCAAGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTg           >  2:696085/2‑150 (MQ=255)
                                                                                                                                              gctggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTc        <  2:1088419/150‑1 (MQ=255)
                                                                                                                                              gctggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTc        >  2:215162/1‑150 (MQ=255)
                                                                                                                                                tggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGt      <  1:5399035/150‑1 (MQ=255)
                                                                                                                                                 ggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTg     <  2:423607/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCAGCAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGC  >  NC_002947/964332‑964629

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: