Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,141,753 T→A 5.7% intergenic (+48/+22) PP_3640 → / ← PP_3641 AraC family transcriptional regulator/cytosine/purine/uracil/thiamine/allantoin permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,141,7530TA 5.7% 499.8 / 15.4 159intergenic (+48/+22)PP_3640/PP_3641AraC family transcriptional regulator/cytosine/purine/uracil/thiamine/allantoin permease family protein
Reads supporting (aligned to +/- strand):  ref base T (70/79);  new base A (5/4);  total (76/83)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

Only 100 of 171 total aligned reads displayed.

GGGTATCAGTCGGTGGCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTCGACAATGTTC  >  NC_002947/4141604‑4141901
                                                                                                                                                     |                                                                                                                                                    
 ggTATCAGTCGGTGGCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCttt                                                                                                                                                     <  2:4706870/150‑1 (MQ=255)
 ggTATCAGTCGGTGGCCGCATTCACCAGGGAGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCttt                                                                                                                                                     <  1:3611203/150‑1 (MQ=255)
  gTATCAGTCGGTGGCCGCATTCACCAGGGCGTTCGCCTAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCTGTCTGCATGACGGGCGCGGGCATTCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTg                                                                                                                                                    <  2:354281/150‑1 (MQ=255)
    aTCAGTCGGTGGCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCa                                                                                                                                                  <  1:1564772/150‑1 (MQ=255)
       aGTCGGTGGCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGcc                                                                                                                                               <  1:833004/150‑1 (MQ=255)
       aGTCGGTGGCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGcc                                                                                                                                               <  2:2011333/150‑1 (MQ=255)
       aGTCGGTGGCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCGTTTCAGcc                                                                                                                                               >  1:318397/1‑150 (MQ=255)
        gTCGGTGGACGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGccc                                                                                                                                              >  2:4525752/1‑150 (MQ=255)
         tCGGTGGCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCAAAGCAGcccc                                                                                                                                             <  1:2731995/150‑1 (MQ=255)
           ggTGGCCGCATTCACCTGGGCGTTCGCCAAGGCTGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGGGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATTGTCGGGGCTGCAAAGCAGccccga                                                                                                                                           <  1:4503601/150‑3 (MQ=255)
            gTGGCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGGCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCTTAGTCGGGGCTGCAAAGCAGccccgac                                                                                                                                          <  2:3660968/150‑4 (MQ=255)
              ggCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGtt                                                                                                                                        <  2:4926764/150‑1 (MQ=255)
              ggCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGtt                                                                                                                                        <  2:2238093/150‑1 (MQ=255)
               gCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTg                                                                                                                                       <  2:1358853/150‑1 (MQ=255)
                ccGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCGGCGTTGCAGCCCCCAGTTgg                                                                                                                                      >  1:70741/1‑149 (MQ=255)
                 cGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCa                                                                                                                                     <  1:243269/150‑1 (MQ=255)
                 cGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCa                                                                                                                                     <  1:5024322/150‑1 (MQ=255)
                   cATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCGTTTCAGCCCCAAGTTGCAGt                                                                                                                                   >  1:5396666/1‑150 (MQ=255)
                    aTTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCAGGTTGCAGtt                                                                                                                                  >  2:4540318/1‑150 (MQ=255)
                     ttCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGGGTTTGATCCCCAAGTTGGCGTTg                                                                                                                                 >  1:2855252/1‑150 (MQ=255)
                        aCCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCa                                                                                                                              >  2:3230007/1‑150 (MQ=255)
                        aCCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCAAGTTGCAGTTGTCa                                                                                                                              >  1:3708184/1‑150 (MQ=255)
                         ccAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAg                                                                                                                             >  2:5560975/1‑150 (MQ=255)
                         ccAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCAAGTTGCAGTTGTCAg                                                                                                                             >  2:3589896/1‑150 (MQ=255)
                          cAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGAAGCCCCAAGTTGCAGTTGTAAGt                                                                                                                            >  1:324476/1‑150 (MQ=255)
                           aGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTc                                                                                                                           >  2:3240437/1‑150 (MQ=255)
                           aGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTc                                                                                                                           <  2:2866787/150‑1 (MQ=255)
                             ggCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAACATGACGGGCGCGGGAATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGt                                                                                                                         <  1:105294/150‑1 (MQ=255)
                               cGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTg                                                                                                                       <  1:4189144/150‑1 (MQ=255)
                               cGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCAGGTTGCCGTTGTCAGTCGTTg                                                                                                                       >  2:2441383/1‑150 (MQ=255)
                                gTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGGCGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGa                                                                                                                      <  2:4792743/150‑1 (MQ=255)
                                    gCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCCTAGTCGGGGCTGCGTTGCAGCCCCAGGTGGCCGGTGTCAGTCGGTGACGcg                                                                                                                  >  2:3296738/1‑149 (MQ=255)
                                    gCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCGTTGCAGCCCCAAGTTTCAGTTGTCAGTCGTTGTCGCt                                                                                                                  >  1:5427061/1‑150 (MQ=255)
                                      cAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGGTGCTTTTCAGCCCCCAGTTGCAGTTTTCAGTCGGTGAAGcagt                                                                                                                >  1:1189451/1‑147 (MQ=255)
                                      cAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCGTTTAAGCCCCAAGATTAAGGTGTCAAACGTTTACGCAGa                                                                                                                >  1:4855748/1‑150 (MQ=255)
                                             gTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCcg                                                                                                         >  2:2709178/1‑150 (MQ=255)
                                              tGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACTAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCcgc                                                                                                        <  2:5559794/150‑1 (MQ=255)
                                               gggggTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCcgcg                                                                                                       <  2:1190565/150‑1 (MQ=255)
                                               gggggTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCcgcg                                                                                                       <  2:4099364/150‑1 (MQ=255)
                                               gggggTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCcgcg                                                                                                       <  2:4098757/150‑1 (MQ=255)
                                                ggggTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTTCAGCCCCATTTTGCGGTTGTCAGTCGTTTACGCTGATAAACCgcgc                                                                                                      >  1:379618/1‑150 (MQ=255)
                                                   gTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCCCCCCCCAGTTGCAGTTGTCAGTCGTTGTCGCTGATAACCCgcgcccg                                                                                                   >  2:2033852/1‑147 (MQ=255)
                                                    tGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTg                                                                                                  >  2:4185748/1‑150 (MQ=255)
                                                      cAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCGTTGCAGCCCCATGGTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACAGtt                                                                                                >  2:4010543/1‑150 (MQ=255)
                                                         ccGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAAc                                                                                             <  2:5417127/150‑1 (MQ=255)
                                                           gggggCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCGGCGTTTGAGCCCCCAGTTTCAGGTGTGAGTCGTTGTAGCGGATAAACCCCCCACAGTTaaacc                                                                                           >  1:4715601/1‑147 (MQ=255)
                                                             gggCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTTCAGCCCCATGATGAAGTTGTCAGTCGTTGAAGCAGATAAAACCCCCACAGTTaaacccc                                                                                         >  1:3572432/1‑145 (MQ=255)
                                                                 gTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCGTTGCAGCCCCAAGTTGCAGTTGTCAGTCGTTGTCGCGGATAACCCCCCCACAGTTAACCCCCATcc                                                                                     >  1:4470423/1‑150 (MQ=255)
                                                                  tATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTccc                                                                                    <  2:1927307/150‑1 (MQ=255)
                                                                     cGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCCAAGCGCGGACATCGACCAAGCGGCCCTAGTCGGGGCGGCGGTGCAGCCCCCGGGTGCCGGTGTCAGGCGGGGCCGCGGAGAACCCGCGCCCGGGGAACCGCCGGCCcgcg                                                                                 >  2:5302754/1‑150 (MQ=255)
                                                                     cGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACGGTTAACCGCCGTCCcgcg                                                                                 >  2:3755563/1‑150 (MQ=255)
                                                                     cGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCAAGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCcgcg                                                                                 >  2:3741255/1‑150 (MQ=255)
                                                                      gcaggcagCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCgt                                                                                <  2:906856/150‑1 (MQ=255)
                                                                      gcaggcagCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCGGCGTTGCAGCCCCCAGGTGCAGTTGTCAGTCGTTGACGCGGATAAACCCCGCCCAGTTAACCCCCATccccccg                                                                                >  1:4919641/1‑146 (MQ=255)
                                                                       caggcagCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCgtg                                                                               <  1:3361217/150‑1 (MQ=255)
                                                                        aggcagCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCgtgt                                                                              <  2:4726435/150‑1 (MQ=255)
                                                                              caTGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAg                                                                        <  2:2577555/150‑1 (MQ=255)
                                                                               aTGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGGTTTGCAGCCCCAAGTTGCAGTTTTCAGTCGGTTAAGCTGATAAACCCCCCCCAGGTAAACCCCATCCCCCCTGTGGCCCgg                                                                       >  1:2970797/1‑150 (MQ=255)
                                                                                tGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGc                                                                      <  2:2550951/150‑1 (MQ=255)
                                                                                      gcgcGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGGGGGGTTTCCGCCCCCAGTTTCAGTTTTCAGTCGTTGTCGCTGATAAACCCCCCCCGGTTAAACCCCATCCCCCGGGTGGGCAGGCCGGAAt                                                                >  1:3057407/1‑150 (MQ=255)
                                                                                        gcgGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCCGTAATAc                                                              >  2:826799/1‑150 (MQ=255)
                                                                                          gggCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCCGCCCCAAGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAg                                                            >  2:2073936/1‑150 (MQ=255)
                                                                                           ggCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGc                                                           <  2:631599/150‑1 (MQ=255)
                                                                                            gCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCGTTGCAGCCCCCTGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAAACGCCGTCCCGCGGGTGGCCCGGCAGTAATACAGCg                                                          >  2:4474784/1‑150 (MQ=255)
                                                                                                gCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCCTGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCCCGGTTAACCGCCTTCCCGCGTGTGGCCAGGCCGTAATACAGCGGCAc                                                      >  2:2866499/1‑150 (MQ=255)
                                                                                                gCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCGTTGCAGCCCCATGGTGCAGTTGTCAGTCGTTGACGCTGCTAACCCGCGCACTGGTAACCGCCTTCCCGCGGGTGGCCAGGCCGTAAAACAGCGGCAc                                                      >  2:802089/1‑150 (MQ=255)
                                                                                                 ccTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCAcc                                                     <  2:97225/150‑1 (MQ=255)
                                                                                                   tGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTTCATTTTTCAGTAGTTGGCGCTGATAACCCGCGCACAGTTAACCGCCTTCCCCCGTGTGGCCAGGCCGTAATACAGCGGCACCGt                                                   >  1:5091041/1‑150 (MQ=255)
                                                                                                       ccAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCCTGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCCCTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACCGCGGCACCGTCAcc                                               >  2:5075217/1‑150 (MQ=255)
                                                                                                        cAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCa                                              <  2:4680583/150‑1 (MQ=255)
                                                                                                        cAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCa                                              <  2:1560491/150‑1 (MQ=255)
                                                                                                        cAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCa                                              <  2:5193516/150‑1 (MQ=255)
                                                                                                         aGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCAAGTTGCAGTTGTCAGTCGTTTACGCGGATAAACCCCGCACAGTTAAACCCCGTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCAcccg                                             >  2:1966689/1‑148 (MQ=255)
                                                                                                            aaaGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAAACCCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCaa                                          >  2:2311373/1‑150 (MQ=255)
                                                                                                                 gcgGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGa                                     <  2:221144/150‑1 (MQ=255)
                                                                                                                  cgGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGaa                                    <  2:4352919/150‑1 (MQ=255)
                                                                                                                   ggAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACCCTGATAAACCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATAAAGCGGCACCCTCACCAGCAAGCCCaaa                                   >  2:4846015/1‑149 (MQ=255)
                                                                                                                    gAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATa                                  <  1:4634589/150‑1 (MQ=255)
                                                                                                                    gAAATCGACCAAGCGGCCATAGTCGGGGCGGCGTTTCCGCCCCAAGTTGCAGTTTTCAGTCGTTGACGCTGATAACCCCCGCACGGTTAACCGCCGTCCCGCGGGTGGCCCGGCCGTAATACAGCGGCACCCTCACCACCACGCCGaaaa                                  >  2:2915686/1‑148 (MQ=255)
                                                                                                                       aTCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGcc                               <  2:1766735/150‑1 (MQ=255)
                                                                                                                         cGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCCCGCACGGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCAc                             >  2:3361389/1‑150 (MQ=255)
                                                                                                                         cGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCACATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCAc                             <  2:3433307/150‑1 (MQ=255)
                                                                                                                            ccAAGCGGCCATAGTCGGGGCTGCGTTGCAGCCCCAGGTTGCACTTGTCAGTCGTTGACCCTGATAACCCGCGCACCGGTAACCGCCATCCCCCGGGTGGCCCCGCCGGAAAACAGCGGCACCCGCACCAGCACGCCCAAAAGCCACGAc                          >  2:3002168/1‑150 (MQ=255)
                                                                                                                             cAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCTCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACa                         <  1:104867/150‑1 (MQ=255)
                                                                                                                                  ggCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCg                    <  2:2075647/150‑1 (MQ=255)
                                                                                                                                  ggCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCg                    <  1:2501460/150‑1 (MQ=255)
                                                                                                                                         catggggCTGCAAAGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCtt             >  2:1683042/4‑150 (MQ=255)
                                                                                                                                         catggggCTGCAAAGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCtt             >  1:3791757/4‑150 (MQ=255)
                                                                                                                                          tCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCATCAAGCCGAATAGCCACGACAGGTCGGCCCCTTc            >  1:890054/1‑150 (MQ=255)
                                                                                                                                          atggggCTGCAAAGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTc            >  1:2055682/3‑150 (MQ=255)
                                                                                                                                          atggggCTGCAAAGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTc            >  1:1638813/3‑150 (MQ=255)
                                                                                                                                            ggggCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTCGa          <  1:2033852/150‑1 (MQ=255)
                                                                                                                                            ggggCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTCGa          <  2:3322381/150‑1 (MQ=255)
                                                                                                                                                   ctTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTCGACAATGtt   <  2:4888190/150‑1 (MQ=255)
                                                                                                                                                    tttGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTCGACAATGTTc  <  1:700213/150‑1 (MQ=255)
                                                                                                                                                    tttGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTCGACAATGTTc  <  2:1111915/150‑1 (MQ=255)
                                                                                                                                                    tttGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTCGACAATGTTc  >  2:1878330/1‑150 (MQ=255)
                                                                                                                                                    tttGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTCGACAATGTTc  <  2:318397/150‑1 (MQ=255)
                                                                                                                                                    tttGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTCGACAATGTTc  <  1:4185625/150‑1 (MQ=255)
                                                                                                                                                    tttGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTCGACAATGTTc  <  1:1494857/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGTATCAGTCGGTGGCCGCATTCACCAGGGCGTTCGCCAAGGCAGTGGGGGTGCAGCCGGGGGCGTATCGCAGGCAGCATGACGGGCGCGGGCATGCCTGACCCAGCAAAGCGCGGAAATCGACCAAGCGGCCATAGTCGGGGCTGCTTTGCAGCCCCATGTTGCAGTTGTCAGTCGTTGACGCTGATAACCCGCGCACTGTTAACCGCCTTCCCGCGTGTGGCCAGGCAGTAATACAGCGGCACCGTCACCAGCAAGCCGAATAGCCACGACAGGTCGGCCCCTTCGACAATGTTC  >  NC_002947/4141604‑4141901

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