Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 964,485 G→A 5.4% L96L (CTG→CTA) ‡ phnE → phosphonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947964,4850GA 5.4% 1078.0 / 31.6 334L96L (CTG→CTA) ‡phnEphosphonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (193/122);  new base A (13/5);  total (207/127)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 343 total aligned reads displayed.

CAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGCGCTG  >  NC_002947/964336‑964633
                                                                                                                                                     |                                                                                                                                                    
  gATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGc                                                                                                                                                    <  1:3680901/150‑1 (MQ=255)
       ccAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGcc                                                                                                                                               >  1:1295071/1‑150 (MQ=255)
        cAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCg                                                                                                                                              >  2:4631939/1‑150 (MQ=255)
        cAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCg                                                                                                                                              >  2:2143238/1‑150 (MQ=255)
         aGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGGTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGt                                                                                                                                             >  2:5399035/1‑150 (MQ=255)
          gCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCCAGCAAGCCAGCCggg                                                                                                                                            <  1:1455625/150‑3 (MQ=255)
           cTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCACGGCTGGCCAGCAAGCCAGCCggga                                                                                                                                           <  1:1650350/150‑4 (MQ=255)
                ccAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTgg                                                                                                                                                     <  2:4110921/135‑1 (MQ=255)
                ccAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGt                                                                                                                                      <  2:2508340/150‑1 (MQ=255)
                 cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                     >  2:4785275/1‑150 (MQ=255)
                 cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                     >  2:2156882/1‑150 (MQ=255)
                 cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                     >  1:998310/1‑150 (MQ=255)
                 cAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                                     >  2:5782309/1‑150 (MQ=255)
                        ttCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGt                                                                                                                              >  2:244383/1‑150 (MQ=255)
                            ggCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGtgctg                                                                                                                          <  2:6183418/150‑1 (MQ=255)
                                 ccgGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTg                                                                                                                                                      >  2:354540/1‑117 (MQ=255)
                                 ccgGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTg                                                                                                                                                      <  1:354540/117‑1 (MQ=255)
                                   tctgacggatgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCAcg                                                                                                                   >  1:91074/8‑150 (MQ=255)
                                   ggCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCAcg                                                                                                                   >  1:1544952/1‑150 (MQ=255)
                                      acaTGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcgg                                                                                                                >  2:1727455/1‑150 (MQ=255)
                                      acaTGATGACGATATTCTTGCACTGCTGTAAGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcgg                                                                                                                >  2:566409/1‑150 (MQ=255)
                                         tgatgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGcggcgg                                                                                                             <  1:2564977/150‑1 (MQ=255)
                                          gatgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcggg                                                                                                            <  2:2630700/150‑1 (MQ=255)
                                           atgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcgggc                                                                                                           <  1:827449/150‑1 (MQ=255)
                                                gATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggcc                                                                                                      >  2:1360942/1‑150 (MQ=255)
                                                gATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggcc                                                                                                      >  1:6773012/1‑150 (MQ=255)
                                                       ttGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGctct                                                                                               <  2:3698581/150‑1 (MQ=255)
                                                       ttGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGctct                                                                                               >  1:5938149/1‑150 (MQ=255)
                                                       ttGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGctct                                                                                               <  1:1928914/150‑1 (MQ=255)
                                                        tGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTg                                                                                              <  2:3742570/150‑1 (MQ=255)
                                                          cACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGt                                                                                            >  2:3017365/1‑150 (MQ=255)
                                                          cACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGACCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGt                                                                                            >  2:5591639/1‑150 (MQ=255)
                                                               ctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGa                                                                                       <  1:5261905/150‑1 (MQ=255)
                                                               ctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCACGGCTGGCCAGCAAGCCAGCCgggatc                                                                                                                                         <  1:4233301/100‑6 (MQ=37)
                                                                 gTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACt                                                                                     <  1:6315094/150‑1 (MQ=255)
                                                                  tATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTg                                                                                    >  2:2760212/1‑150 (MQ=255)
                                                                   aTGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGt                                                                                   >  1:5027498/1‑150 (MQ=255)
                                                                    tGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGtt                                                                                  >  1:7227916/1‑150 (MQ=255)
                                                                     gAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTg                                                                                 >  1:4335395/1‑150 (MQ=255)
                                                                     gAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTg                                                                                 >  2:6410240/1‑150 (MQ=255)
                                                                        gCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTgcgc                                                                              <  2:4008600/150‑1 (MQ=255)
                                                                            ccTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGc                                                                          >  2:324882/1‑150 (MQ=255)
                                                                              ttGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCg                                                                        >  2:6564228/1‑150 (MQ=255)
                                                                                gCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGt                                                                      <  2:4444056/150‑1 (MQ=255)
                                                                                gCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGt                                                                      >  1:4648334/1‑150 (MQ=255)
                                                                                gCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGt                                                                      >  2:1463116/1‑150 (MQ=255)
                                                                                  aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCtt                                                                                                                       >  2:2965928/1‑99 (MQ=255)
                                                                                  aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCtt                                                                                                                       <  1:2965928/99‑1 (MQ=255)
                                                                                  aGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGc                                                                    >  1:424889/1‑150 (MQ=255)
                                                                                  aGACCCTGGCCATGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCGTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGc                                                                    >  2:6425834/1‑150 (MQ=255)
                                                                                      ccTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTgg                                                                >  1:5512157/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)
                                                                                                  gacagacGGCATGCCGCTGGCGTTGTTGTTGGCGATCCCTGCTCGATAGCTTGCCAGCCGTGCATTGTCCTTGCGTGCTGCGTCACGCGGAGGGCGGGCCGGGCTAGCGTCGCCACTGTTGCGCCTTCAGGAGCGTGGTCAACTTAtttt                                                    >  2:4973919/7‑150 (MQ=255)
                                                                                                  gACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGAtttt                                                    >  2:187098/1‑150 (MQ=255)
                                                                                                  gACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGAtttt                                                    >  2:5908770/1‑150 (MQ=255)
                                                                                                  gACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGAtttt                                                    >  2:615071/1‑150 (MQ=255)
                                                                                                   aCTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGGGCGTGGGCTACTGATTTTc                                                   >  1:3403382/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)
                                                                                                     tGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCt                                                 >  2:3359796/1‑150 (MQ=255)
                                                                                                         ggCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGt                                             <  2:6740573/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                                        <  2:7124312/150‑1 (MQ=255)
                                                                                                              ggCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGt                                        <  1:4076261/150‑1 (MQ=255)
                                                                                                              ggCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGt                                        <  1:5808996/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)
                                                                                                                           ttgttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTc                           <  1:6009933/150‑1 (MQ=255)
                                                                                                                           ttgttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCGCTGGTCGCGACTGTTGCGCCGGCCGGGGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCggc                           >  1:554031/1‑148 (MQ=255)
                                                                                                                             gttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGGGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTg                         >  2:5228908/1‑150 (MQ=255)
                                                                                                                              ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                        <  1:4890958/150‑1 (MQ=255)
                                                                                                                              ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                        >  1:1408223/1‑150 (MQ=255)
                                                                                                                              ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                        >  1:6415433/1‑150 (MQ=255)
                                                                                                                              ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                        <  1:717274/150‑1 (MQ=255)
                                                                                                                              ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                        <  1:2408685/150‑1 (MQ=255)
                                                                                                                              ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCGCTGGGCGCGGCTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTgg                        >  1:1001430/1‑150 (MQ=255)
                                                                                                                                    aTCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAACTCGTCTGGGCGTTg                  >  2:4588858/1‑150 (MQ=255)
                                                                                                                                     gcacGGCTGGCCAGCAAGCCAGCCGGGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGc                 >  1:2447305/4‑150 (MQ=255)
                                                                                                                                      cccGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCt                >  2:3086044/1‑150 (MQ=255)
                                                                                                                                       gcGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCAGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGGTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTg               >  2:2060737/2‑150 (MQ=255)
                                                                                                                                       acGGCTGGCCAGCAAGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGCCTCGGGGCGCGACTGTTGCGCCAGCCGGTGCGTGGGCTACTGAGTTTCCTGCGTAGCGGGCCGGAAAGCGGCGGGGCGGTGCTg               >  2:695727/2‑150 (MQ=255)
                                                                                                                                          gctggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTc            <  1:1199327/150‑1 (MQ=255)
                                                                                                                                          gctggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCAGCCGGTGCGTGGGCTACTGATTTTACTGCGTAGCGTGCCGGAAATCGTCTGGGCGATGCTGTTc            >  1:4741270/1‑150 (MQ=255)
                                                                                                                                           ctggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCGCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGGTACTGATTTTACTGCGGAGCGTGCCGGAAATCGTCTGGGCGTTGAGGTTCg           >  1:2876636/1‑150 (MQ=255)
                                                                                                                                            tggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGt          <  2:4876605/150‑1 (MQ=255)
                                                                                                                                            tggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGt          <  1:5399029/150‑1 (MQ=255)
                                                                                                                                             ggctTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTg         <  2:423607/150‑1 (MQ=255)
                                                                                                                                                    tGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCGCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGCGCtg  >  2:5328691/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGACTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGCGTTGCTGTTCGTGCGCGCTG  >  NC_002947/964336‑964633

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: