Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,863,038 G→C 5.4% R32P (CGC→CCC)  PP_4273 → MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,863,0380GC 5.4% 935.1 / 31.7 276R32P (CGC→CCC) PP_4273MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (125/135);  new base C (4/11);  total (129/147)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 287 total aligned reads displayed.

ATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGAC  >  NC_002947/4862889‑4863186
                                                                                                                                                     |                                                                                                                                                    
aTTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGgcg                                                                                                                                                      >  2:2131210/1‑150 (MQ=255)
  tACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGcgcg                                                                                                                                                    <  2:2342793/150‑1 (MQ=255)
    ccATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGc                                                                                                                                                  >  2:4839042/1‑150 (MQ=255)
    ccATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGc                                                                                                                                                  >  1:1747867/1‑150 (MQ=255)
    ccATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGc                                                                                                                                                  <  2:4795324/150‑1 (MQ=255)
     cATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGcc                                                                                                                                                 >  2:3146461/1‑150 (MQ=255)
       tGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGCGCCCGGCCaa                                                                                                                                               <  1:1311802/150‑2 (MQ=255)
       tGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGCGCCCGGCCaa                                                                                                                                               <  2:1848579/150‑2 (MQ=255)
        gCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGCGCCCGGCCaag                                                                                                                                              <  1:4595252/150‑3 (MQ=255)
        gCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGCGCCCGGCCaag                                                                                                                                              <  1:107282/150‑3 (MQ=255)
          cgcgCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGt                                                                                                                                            >  2:4426867/1‑150 (MQ=255)
          cgcgCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGCGCCCGGCCaagga                                                                                                                                            <  1:581991/150‑5 (MQ=255)
           gcgcAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTc                                                                                                                                           <  2:406121/150‑1 (MQ=255)
                   tAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGttt                                                                                                                                   >  2:2633438/1‑150 (MQ=255)
                    aGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTg                                                                                                                                  >  1:2521975/1‑150 (MQ=255)
                     gAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGc                                                                                                                                 <  2:4663061/150‑1 (MQ=255)
                      aaGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCa                                                                                                                                >  2:2186329/1‑150 (MQ=255)
                      aaGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCa                                                                                                                                <  1:1564106/150‑1 (MQ=255)
                           gAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGa                                                                                                                           >  1:1651570/1‑150 (MQ=255)
                              gTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGAtgt                                                                                                                        >  2:2234425/1‑150 (MQ=255)
                              gTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGAtgt                                                                                                                        <  1:990803/150‑1 (MQ=255)
                              gTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGAtgt                                                                                                                        >  2:1118351/1‑150 (MQ=255)
                                cAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTgg                                                                                                                      <  2:5127108/150‑1 (MQ=255)
                                cAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTgg                                                                                                                      <  1:2874493/150‑1 (MQ=255)
                                 aaCTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGTCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGc                                                                                                                     >  1:5358675/1‑150 (MQ=255)
                                    ttCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGa                                                                                                                  >  1:1434475/1‑150 (MQ=255)
                                     tcaagtcaagGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGAc                                                                                                                 <  2:4533830/150‑1 (MQ=255)
                                      caagtcaagGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGAcc                                                                                                                >  1:4483688/1‑150 (MQ=255)
                                         gtcaagGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGt                                                                                                             >  2:387264/1‑150 (MQ=255)
                                         gtcaagGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGt                                                                                                             <  2:888388/150‑1 (MQ=255)
                                           caagGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTg                                                                                                           <  2:3651311/150‑1 (MQ=255)
                                             agGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGgc                                                                                                         <  2:2510063/150‑1 (MQ=255)
                                             agGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGgc                                                                                                         >  1:4126388/1‑150 (MQ=255)
                                              gggCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGgcg                                                                                                        <  1:1600386/150‑1 (MQ=255)
                                                  ccGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGa                                                                                                    <  2:3177321/150‑1 (MQ=255)
                                                       aTGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGcc                                                                                               <  2:2556142/150‑1 (MQ=255)
                                                           aCATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTgg                                                                                           >  1:3904807/1‑150 (MQ=255)
                                                            cATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTggg                                                                                          <  2:3423829/150‑1 (MQ=255)
                                                               tGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCa                                                                                       <  2:4950039/150‑1 (MQ=255)
                                                               ggCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCa                                                                                       <  1:1559349/149‑1 (MQ=255)
                                                                      gTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCg                                                                                <  2:1672069/150‑1 (MQ=255)
                                                                       tCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCgg                                                                               >  2:1026504/1‑150 (MQ=255)
                                                                            aaaCGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGttt                                                                          >  1:4938906/1‑150 (MQ=255)
                                                                             aaCGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGtttt                                                                         <  2:1714112/150‑1 (MQ=255)
                                                                              aCGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTc                                                                        >  2:5602062/1‑150 (MQ=255)
                                                                                  aCTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTg                                                                    <  1:5488451/150‑1 (MQ=255)
                                                                                  aCTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTg                                                                    <  1:2235907/150‑1 (MQ=255)
                                                                                  aCTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTg                                                                    >  1:878557/1‑150 (MQ=255)
                                                                                  aCTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCGTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTg                                                                    <  2:3652658/150‑1 (MQ=255)
                                                                                   cTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTgg                                                                   <  2:5210143/150‑1 (MQ=255)
                                                                                    tGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGa                                                                  <  1:2684697/150‑1 (MQ=255)
                                                                                    gggCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGa                                                                  <  1:3086357/149‑1 (MQ=255)
                                                                                           ccTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTg                                                           <  1:3719183/150‑1 (MQ=255)
                                                                                            cTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTgg                                                          <  1:4316235/150‑1 (MQ=255)
                                                                                            cTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTgg                                                          >  1:3428474/1‑150 (MQ=255)
                                                                                             tGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTggg                                                         <  1:3841106/150‑1 (MQ=255)
                                                                                              gCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTgggg                                                        >  1:2016629/1‑150 (MQ=255)
                                                                                              gCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTgggg                                                        <  1:2633438/150‑1 (MQ=255)
                                                                                              gCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTgggg                                                        >  1:643817/1‑150 (MQ=255)
                                                                                               cATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTggggg                                                       <  2:4008690/150‑1 (MQ=255)
                                                                                                    acacTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATg                                                  >  1:1917362/1‑150 (MQ=255)
                                                                                                      acTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCt                                                >  1:1495419/1‑150 (MQ=255)
                                                                                                        tGCGTTACTACGAGAAGAAAGGCTTTCTCAAGTCCTTGGACGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGAGCTGATCGCCGTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCtgt                                              >  1:2381987/1‑149 (MQ=255)
                                                                                                        tGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCtgg                                              <  2:5372717/150‑1 (MQ=255)
                                                                                                         gCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCtggt                                             <  1:4047508/150‑1 (MQ=255)
                                                                                                          cGTTACTACGAGAAGAAAGGCTTGCTCGAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCtggtg                                            >  2:2870905/1‑150 (MQ=255)
                                                                                                          cGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCtggtg                                            <  1:467578/150‑1 (MQ=255)
                                                                                                             tactacGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGAc                                         <  2:1747867/150‑1 (MQ=255)
                                                                                                             tactacGAGAAGAAAGGCTTGCTCAAGTCCTGGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGAc                                         <  2:5282378/150‑1 (MQ=255)
                                                                                                                tacGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCATCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTg                                      >  1:1199547/1‑150 (MQ=255)
                                                                                                                tacGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTg                                      >  2:3928114/1‑150 (MQ=255)
                                                                                                                tacGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTg                                      >  1:5515464/1‑150 (MQ=255)
                                                                                                                   gagaAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAg                                   <  1:392755/150‑1 (MQ=255)
                                                                                                                        gaaAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGGTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGGTTTCGCTGGATGAGGTGGGGGCGATGCTGGTGGACCTGCTGGTCGa                              >  2:136038/1‑150 (MQ=255)
                                                                                                                         aaaGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGAc                             >  1:1049446/1‑150 (MQ=255)
                                                                                                                         aaaGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGAc                             >  1:4804693/1‑150 (MQ=255)
                                                                                                                           aGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCg                           >  2:4322066/1‑150 (MQ=255)
                                                                                                                            ggCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCgg                          <  1:1051951/150‑1 (MQ=255)
                                                                                                                             gCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGc                         <  2:3361474/150‑1 (MQ=255)
                                                                                                                              cTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCa                        <  2:1434475/150‑1 (MQ=255)
                                                                                                                              cTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCCGGTCGACCTGCAGGTCCACCAGCa                        <  1:318284/150‑1 (MQ=255)
                                                                                                                               ttGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAg                       >  1:257553/1‑150 (MQ=255)
                                                                                                                                tGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCCGGTCGACCTGCAGGTCCACCAGCatc                      <  1:2131210/150‑3 (MQ=255)
                                                                                                                                  cTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATg                    <  1:5220667/150‑1 (MQ=255)
                                                                                                                                   tCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGc                   >  2:1951009/1‑150 (MQ=255)
                                                                                                                                   tCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGc                   >  2:711714/1‑150 (MQ=255)
                                                                                                                                   tCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGc                   >  1:822989/1‑150 (MQ=255)
                                                                                                                                   tCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGc                   >  1:5465928/1‑150 (MQ=255)
                                                                                                                                   tCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGc                   >  1:3587621/1‑150 (MQ=255)
                                                                                                                                     aaGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTg                 >  2:1575295/1‑150 (MQ=255)
                                                                                                                                        tCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATg              >  2:854868/1‑150 (MQ=255)
                                                                                                                                         ccTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATgg             <  2:644921/150‑1 (MQ=255)
                                                                                                                                           ttGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGAc                             >  1:400871/1‑132 (MQ=255)
                                                                                                                                           ttGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGcc           >  1:5069484/1‑150 (MQ=255)
                                                                                                                                           ctGGCCGCGCCCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGAGGCCGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCCGATGCTGATGGcc           >  2:129299/2‑150 (MQ=255)
                                                                                                                                              gCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGAGGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAg        >  2:2967403/1‑150 (MQ=255)
                                                                                                                                                cGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGc      <  1:677677/150‑1 (MQ=255)
                                                                                                                                                cGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGc      <  2:105286/150‑1 (MQ=255)
                                                                                                                                                   gcgcGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGa   <  1:4426867/150‑1 (MQ=255)
                                                                                                                                                    cgcgGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGAc  <  1:4828739/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGAC  >  NC_002947/4862889‑4863186

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: