Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,938,985 A→T 5.5% L474Q (CTG→CAG)  PP_3467 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,938,9850AT 5.5% 905.8 / 26.9 264L474Q (CTG→CAG) PP_3467sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (135/114);  new base T (4/11);  total (139/125)
Fisher's exact test for biased strand distribution p-value = 5.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

CAGCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCAGCCGC  >  NC_002947/3938837‑3939134
                                                                                                                                                    |                                                                                                                                                     
cagcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAg                                                                                                                                                      <  2:3980926/150‑1 (MQ=255)
 agcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgc                                                                                                                                                     <  2:2487835/150‑1 (MQ=255)
 agcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgc                                                                                                                                                     <  1:4326732/150‑1 (MQ=255)
  gcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgcg                                                                                                                                                    >  1:2964931/1‑150 (MQ=255)
   cagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgcgc                                                                                                                                                   <  1:4968974/150‑1 (MQ=255)
   cagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgcgc                                                                                                                                                   <  1:4759508/150‑1 (MQ=255)
   cagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgcgc                                                                                                                                                   <  1:2285337/150‑1 (MQ=255)
     ggCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCTGATCTGCGCCa                                                                                                                                                 <  1:243609/150‑1 (MQ=255)
     ggCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCTGATCTGCGCCa                                                                                                                                                 <  1:3462542/150‑1 (MQ=255)
        ccTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGt                                                                                                                                              >  1:1953897/1‑150 (MQ=255)
        ccTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGt                                                                                                                                              <  1:4648648/150‑1 (MQ=255)
        ccTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGt                                                                                                                                              <  1:1815331/150‑1 (MQ=255)
         cTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCTGATCTGCGCCAccca                                                                                                                                             <  1:422238/150‑4 (MQ=255)
         cTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGTCACGTc                                                                                                                                             <  2:5499138/150‑1 (MQ=255)
         cTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGTCACGTc                                                                                                                                             <  2:2637049/150‑1 (MQ=255)
         cTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTc                                                                                                                                             >  2:4012345/1‑150 (MQ=255)
         cTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTc                                                                                                                                             <  1:1452068/150‑1 (MQ=255)
               gatCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTc                                                                                                                                       <  2:3267825/150‑1 (MQ=255)
                 tCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCt                                                                                                                                     >  1:1995831/1‑150 (MQ=255)
                   gccgccAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGc                                                                                                                                   >  1:956400/1‑150 (MQ=255)
                     cgccAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTc                                                                                                                                 <  2:2167349/150‑1 (MQ=255)
                       ccAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGt                                                                                                                               >  1:5363898/1‑150 (MQ=255)
                           gtcgtcGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCAcc                                                                                                                           <  2:2232759/150‑1 (MQ=255)
                           gtcgtcGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCAcc                                                                                                                           <  2:2791742/150‑1 (MQ=255)
                           gtcgtcGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCAcc                                                                                                                           <  1:4405548/150‑1 (MQ=255)
                                cgcgTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGccc                                                                                                                      <  2:2089818/150‑1 (MQ=255)
                                   gTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCa                                                                                                                   >  1:699850/1‑150 (MQ=255)
                                   gTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCa                                                                                                                   >  1:1156345/1‑150 (MQ=255)
                                       aCGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAggg                                                                                                               <  1:3649035/150‑1 (MQ=255)
                                         gCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTg                                                                                                             >  1:1284797/1‑150 (MQ=255)
                                          cAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGc                                                                                                            <  1:41023/150‑1 (MQ=255)
                                            gggCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTa                                                                                                          >  1:5405623/1‑150 (MQ=255)
                                              gCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcg                                                                                                        <  1:382427/150‑1 (MQ=255)
                                                        gcAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCt                                                                                              <  2:1727538/150‑1 (MQ=255)
                                                        gcAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCt                                                                                              >  2:1693435/1‑150 (MQ=255)
                                                           gCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGACGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCa                                                                                           <  2:1313594/150‑1 (MQ=255)
                                                               cacacCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTAc                                                                                       >  1:3365191/1‑150 (MQ=255)
                                                                acacCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTAcc                                                                                      >  2:4093068/1‑150 (MQ=255)
                                                                  acCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGccc                                                                                                                      <  2:2874104/116‑1 (MQ=255)
                                                                  acCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCg                                                                                    >  1:480520/1‑150 (MQ=255)
                                                                  acCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCg                                                                                    >  1:1790019/1‑150 (MQ=255)
                                                                    cGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCa                                                                                  >  1:4732654/1‑150 (MQ=255)
                                                                           ccGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCCAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCg                                                                           >  2:5241272/1‑150 (MQ=255)
                                                                           ccGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCg                                                                           >  1:3148503/1‑150 (MQ=255)
                                                                                 gAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTg                                                                     <  2:1022151/150‑1 (MQ=255)
                                                                                   aCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCa                                                                   <  2:2351709/150‑1 (MQ=255)
                                                                                   aCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCa                                                                   <  2:4835512/150‑1 (MQ=255)
                                                                                     aGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAgg                                                                 <  1:2816401/150‑1 (MQ=255)
                                                                                       gTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGc                                                               >  2:5531919/1‑150 (MQ=255)
                                                                                         cTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGa                                                             >  1:1762873/1‑150 (MQ=255)
                                                                                         cTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGa                                                             <  2:2319100/150‑1 (MQ=255)
                                                                                          ttCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAg                                                            >  2:5054988/1‑150 (MQ=255)
                                                                                          ttCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAg                                                            >  1:5261081/1‑150 (MQ=255)
                                                                                           tCACGGAAGTTCTGCTCCTGCACCTGGCGCTTGAGGTCGCGGTGGGTGGCGGAGATCATCGCCTCGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGATCCTGCAGTACCCTCTAATGGCGTGCTTGCAGTGCGAGc                                                           <  2:2257901/150‑1 (MQ=255)
                                                                                           tCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGc                                                           <  1:638827/150‑1 (MQ=255)
                                                                                           tCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGc                                                           >  2:2860079/1‑150 (MQ=255)
                                                                                             aCGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCgg                                                         >  2:4903318/1‑150 (MQ=255)
                                                                                             aCGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCgg                                                         >  2:5078651/1‑150 (MQ=255)
                                                                                             aCGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCgg                                                         >  1:2655665/1‑150 (MQ=255)
                                                                                              cGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGc                                                        <  2:1442863/150‑1 (MQ=255)
                                                                                                gAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCAt                                                      <  1:2337073/150‑1 (MQ=255)
                                                                                                gAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCAt                                                      <  1:3247572/150‑1 (MQ=255)
                                                                                                   gTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTc                                                   <  2:4304631/150‑1 (MQ=255)
                                                                                                    tgctgcTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCg                                                  >  2:119086/1‑150 (MQ=255)
                                                                                                       tgctCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCg                                               >  1:866458/1‑150 (MQ=255)
                                                                                                       tgctCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCg                                               >  2:2091229/1‑150 (MQ=255)
                                                                                                        gctCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGa                                              >  2:5583780/1‑150 (MQ=255)
                                                                                                        gctCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGa                                              >  1:2737098/1‑150 (MQ=255)
                                                                                                         ctCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGAt                                             >  1:5024534/1‑150 (MQ=255)
                                                                                                             tGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTtcg                                         >  1:5254338/1‑150 (MQ=255)
                                                                                                              gCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATAAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTtcgt                                        >  1:2640595/1‑150 (MQ=255)
                                                                                                               cACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTtcgtc                                       <  1:4527778/150‑1 (MQ=255)
                                                                                                                aCCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTtcgtcg                                      <  2:4204982/150‑1 (MQ=255)
                                                                                                                 ccAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGAACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGa                                     <  1:3175572/150‑1 (MQ=255)
                                                                                                                  cAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCACGGGTGCTACGCGGCGCTCCTGCAGTACCCGCACCAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAg                                    >  2:159345/1‑150 (MQ=255)
                                                                                                                  cAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAg                                    >  1:2725257/1‑150 (MQ=255)
                                                                                                                  cAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAg                                    >  1:2887908/1‑150 (MQ=255)
                                                                                                                  cAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAg                                    >  1:451654/1‑150 (MQ=255)
                                                                                                                  cAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAg                                    >  1:4427235/1‑150 (MQ=255)
                                                                                                                  cAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAg                                    <  2:1406140/150‑1 (MQ=255)
                                                                                                                   aGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAgg                                   <  2:3647346/150‑1 (MQ=255)
                                                                                                                        cTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGt                                                                          >  2:1792135/1‑106 (MQ=255)
                                                                                                                         ttGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGc                             >  1:433966/1‑150 (MQ=255)
                                                                                                                         ttGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATGTCGTCGAGGAACAGc                             >  2:867539/1‑150 (MQ=255)
                                                                                                                           gAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGt                           <  1:5595754/150‑1 (MQ=255)
                                                                                                                             ggTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGc                             <  2:370247/146‑1 (MQ=255)
                                                                                                                             ggTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGc                             >  1:370247/1‑146 (MQ=255)
                                                                                                                                   ggtgggtgGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCAt                   >  2:3895068/1‑150 (MQ=255)
                                                                                                                                   ggtgggtgGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCAt                   >  2:2525010/1‑150 (MQ=255)
                                                                                                                                    gtgggtggCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATg                  <  2:2965079/150‑1 (MQ=255)
                                                                                                                                     tgggtggCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATgg                 <  2:4609185/150‑1 (MQ=255)
                                                                                                                                      acgtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATggg                >  1:5111346/3‑150 (MQ=255)
                                                                                                                                      acgtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATggg                >  1:5001897/3‑150 (MQ=255)
                                                                                                                                        gtggCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGcc              >  1:4194153/1‑150 (MQ=255)
                                                                                                                                         tggCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCt             <  1:2193416/150‑1 (MQ=255)
                                                                                                                                           gcgcAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGc           >  1:3070653/1‑150 (MQ=255)
                                                                                                                                           gcgcAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGc           >  1:2675533/1‑150 (MQ=255)
                                                                                                                                             gcAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCAt         >  1:632607/1‑150 (MQ=255)
                                                                                                                                                 aTCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCAgc     >  2:1271561/1‑150 (MQ=255)
                                                                                                                                                    aGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCAgccgc  <  2:4855148/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCAGCCGC  >  NC_002947/3938837‑3939134

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