Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 964,480 T→C 6.4% L95L (TTG→CTG) ‡ phnE → phosphonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947964,4800TC 6.4% 1004.3 / 36.4 313L95L (TTG→CTG) ‡phnEphosphonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (179/113);  new base C (13/7);  total (193/120)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

Only 100 of 339 total aligned reads displayed.

GTCAGCAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGC  >  NC_002947/964331‑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)
            ccAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGcc                                                                                                                                           >  1:1295071/1‑150 (MQ=255)
            ccAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCCAGCAagccagcc                                                                                                                                           <  1:7239493/150‑1 (MQ=255)
             cAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCg                                                                                                                                          >  2:6077549/1‑150 (MQ=255)
               gCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCCAGCAAGCCAGCCggg                                                                                                                                        <  1:1455625/150‑3 (MQ=255)
                cTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCACGGCTGGCCAGCAAGCCAGCCggga                                                                                                                                       <  1:1650350/150‑4 (MQ=255)
                     ccAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGt                                                                                                                                  <  2:2508340/150‑1 (MQ=255)
                      cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                 >  1:998310/1‑150 (MQ=255)
                      cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                 >  2:2111231/1‑150 (MQ=255)
                      cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                 >  2:2156882/1‑150 (MQ=255)
                      cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                 >  2:4785352/1‑150 (MQ=255)
                      cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                 >  2:5782220/1‑150 (MQ=255)
                         cgcCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCtt                                                                                                                              <  1:4987605/150‑1 (MQ=255)
                           ccTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGc                                                                                                                            <  2:5605856/150‑1 (MQ=255)
                           ccTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGc                                                                                                                            >  1:5941243/1‑150 (MQ=255)
                             ttCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGt                                                                                                                          <  2:4040765/150‑1 (MQ=255)
                             ttCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGt                                                                                                                          >  2:244383/1‑150 (MQ=255)
                              acTGGCCGCCGGCACAAGATGCCGATTTTCTTGCACTGCAGAATGAAGCTACCTTGCAGACCCTGCCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCCGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCgag                                                                                                                         <  2:6425807/149‑3 (MQ=255)
                                tGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGtgct                                                                                                                       >  1:1509740/1‑150 (MQ=255)
                                   ccgccgGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTgctgct                                                                                                                    >  1:5547619/1‑150 (MQ=255)
                                     gccgGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTc                                                                                                                  >  1:1064637/1‑150 (MQ=255)
                                        tctgacggatgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCAcg                                                                                                               >  1:91074/8‑150 (MQ=255)
                                        ggCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCAcg                                                                                                               >  2:6674508/1‑150 (MQ=255)
                                           acaTGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcgg                                                                                                            >  2:1727455/1‑150 (MQ=255)
                                                atgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcgggc                                                                                                       <  1:827449/150‑1 (MQ=255)
                                                atgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcgggc                                                                                                       >  2:4596732/1‑150 (MQ=255)
                                                     gATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggcc                                                                                                  >  1:6773012/1‑150 (MQ=255)
                                                     gATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggcc                                                                                                  >  2:6203616/1‑150 (MQ=255)
                                                     gATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggcc                                                                                                  >  2:1360942/1‑150 (MQ=255)
                                                       ttttCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggccgg                                                                                                <  2:5018427/150‑1 (MQ=255)
                                                        tttCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggccggg                                                                                               <  2:4909728/150‑1 (MQ=255)
                                                            ttGCACTGCTGTATGTAGCTAGCTTGCAGACCGTGGCCGTGGCGAGTGCCGGCATGGCGGTGGGGTTGTTGTTGGCGATCCCGGGTGGCTTGCTGGCGAGCCGTGCGTTGTGCTTGCGTGCTGCTTCTGGCGGCGGGCGGGGCGGGgtct                                                                                           >  2:3999669/1‑150 (MQ=255)
                                                             tGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTg                                                                                          <  2:4715931/150‑1 (MQ=255)
                                                               cACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGt                                                                                        >  2:3017365/1‑150 (MQ=255)
                                                               cACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGGTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGt                                                                                        >  2:639259/1‑150 (MQ=255)
                                                                  tgctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcgggcg                                                                                                      >  1:3587917/1‑133 (MQ=255)
                                                                  tgctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTcgc                                                                                     >  1:5828522/1‑150 (MQ=255)
                                                                  tgctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTcgc                                                                                     <  2:5347500/150‑1 (MQ=255)
                                                                    ctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCACGGCTGGCCAGCAAGCCAGCCgggatc                                                                                                                                     <  1:4233301/100‑6 (MQ=37)
                                                                    ctgAATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGa                                                                                   >  2:2966773/1‑150 (MQ=255)
                                                                       tATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTg                                                                                >  2:2760212/1‑150 (MQ=255)
                                                                         tGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGtt                                                                              <  2:3699750/150‑1 (MQ=255)
                                                                          gAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTg                                                                             >  1:4335395/1‑150 (MQ=255)
                                                                           aaGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTgc                                                                            <  1:5396315/150‑1 (MQ=255)
                                                                           aaGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTgc                                                                            <  2:4835444/150‑1 (MQ=255)
                                                                               tACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCt                                                                        <  2:3619119/150‑1 (MQ=255)
                                                                                 ccTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGGTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGc                                                                      >  2:324987/1‑150 (MQ=255)
                                                                                   ttGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCg                                                                    >  2:6564228/1‑150 (MQ=255)
                                                                                   ttGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCg                                                                    >  2:5866312/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)
                                                                                       aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGc                                                                <  2:4292306/150‑1 (MQ=255)
                                                                                       aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGc                                                                >  2:5928126/1‑150 (MQ=255)
                                                                                       aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGc                                                                >  1:424889/1‑150 (MQ=255)
                                                                                           ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTgg                                                            >  1:1828601/1‑150 (MQ=255)
                                                                                           ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTgg                                                            >  2:3532648/1‑150 (MQ=255)
                                                                                            cTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTggg                                                           <  1:7000143/150‑1 (MQ=255)
                                                                                             tGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGAGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGc                                                          <  2:1168510/150‑1 (MQ=255)
                                                                                                ccGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTAc                                                       >  1:5884576/1‑150 (MQ=255)
                                                                                                  gTGGCGACTGCCGGCATGGCGGTGGCGTTGGTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTg                                                     >  2:5236/1‑150 (MQ=255)
                                                                                                       gACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGAtttt                                                >  2:187098/1‑150 (MQ=255)
                                                                                                       gACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGAtttt                                                >  2:615071/1‑150 (MQ=255)
                                                                                                         cTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTcc                                              <  2:7437132/150‑1 (MQ=255)
                                                                                                          tGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGTCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGTCTACTGATTTTCCt                                             <  1:4565206/150‑1 (MQ=255)
                                                                                                          tGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCt                                             <  1:2714124/150‑1 (MQ=255)
                                                                                                               gCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTa                                        >  1:6180237/1‑150 (MQ=255)
                                                                                                               gCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGACCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTACTGCGTa                                        >  1:239812/1‑150 (MQ=255)
                                                                                                                  tGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTTCGTAGCg                                     >  1:5720914/1‑150 (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)
                                                                                                                         ggCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGa                              >  2:1819747/1‑150 (MQ=255)
                                                                                                                          gCGTTGTGGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGGGCCTTGTCCGTGCGTGGTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGGTGCGCCAGCCGGTGCGTGGGCTCCTGATTTTCCTGCGTAGCGTGCCGGac                             >  2:420351/1‑149 (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)
                                                                                                                                   ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                    >  1:1408223/1‑150 (MQ=255)
                                                                                                                                   ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                    <  1:717274/150‑1 (MQ=255)
                                                                                                                                   ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                    <  1:4890958/150‑1 (MQ=255)
                                                                                                                                   ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                    >  1:6415433/1‑150 (MQ=255)
                                                                                                                                   ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                    >  1:5512600/1‑150 (MQ=255)
                                                                                                                                   ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCGCTGGGCGCGGCTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                    >  1:1001430/1‑150 (MQ=255)
                                                                                                                                      gCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCg                 >  2:2658081/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)
                                                                                                                                            gcGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCAGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGGTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTg           >  2:2060737/2‑150 (MQ=255)
                                                                                                                                            acGGCTGGCCAGCAAGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTg           >  2:696085/2‑150 (MQ=255)
                                                                                                                                               gctggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCAGCCGGTGCGTGGGCTACTGATTTTACTGCGTAGCGTGCCGGAAATCGTCTGGGCGATGCTGTTc        >  1:4741270/1‑150 (MQ=255)
                                                                                                                                                 tggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGt      <  1:5399035/150‑1 (MQ=255)
                                                                                                                                                 tggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGt      >  2:5459508/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCAGCAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGC  >  NC_002947/964331‑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: