Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,687,926 A→G 5.9% Q485R (CAG→CGG)  recJ → ssDNA‑specific exonuclease RecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,687,9260AG 5.9% 1030.6 / 31.7 322Q485R (CAG→CGG) recJssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base A (152/151);  new base G (6/13);  total (158/164)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 325 total aligned reads displayed.

GGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAGTGTGGC  >  NC_002947/1687777‑1688075
                                                                                                                                                     |                                                                                                                                                     
  cGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGc                                                                                                                                                     <  1:326205/150‑1 (MQ=255)
  cGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGc                                                                                                                                                     <  1:2089315/150‑1 (MQ=255)
    tttGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCCATGGGCCGGCat                                                                                                                                                   <  1:5490431/150‑2 (MQ=255)
    tttGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGCAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCAc                                                                                                                                                   >  2:3913974/1‑150 (MQ=255)
      tGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAACTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCTGGCCCCATGGGCCGGCatgg                                                                                                                                                 <  1:3871045/150‑4 (MQ=255)
        ccGAGGCGTTCGACGAAGAGGTGCGGCGTCAGTTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCTGGCCCCATGGGCCGGCa                                                                                                                                                    >  1:4456346/1‑145 (MQ=255)
        ccGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTcc                                                                                                                                               >  2:2038515/1‑150 (MQ=255)
        ccGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTcc                                                                                                                                               <  2:2666090/150‑1 (MQ=255)
        ccGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTcc                                                                                                                                               >  2:1633924/1‑150 (MQ=255)
        ccGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTcc                                                                                                                                               >  1:2289859/1‑150 (MQ=255)
          gAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCt                                                                                                                                             >  1:5062302/1‑150 (MQ=255)
            ggCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGa                                                                                                                                           <  2:6299568/150‑1 (MQ=255)
                ttCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCg                                                                                                                                       >  2:1071085/1‑150 (MQ=255)
                 tCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGc                                                                                                                                      <  1:3714992/150‑1 (MQ=255)
                  cgacgaAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCt                                                                                                                                     <  1:2679574/150‑1 (MQ=255)
                   gacgaAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTg                                                                                                                                    >  2:1554977/1‑150 (MQ=255)
                      gaAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTc                                                                                                                                 >  2:339952/1‑150 (MQ=255)
                       aaGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTcc                                                                                                                                <  2:6752268/150‑1 (MQ=255)
                         gagGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCAt                                                                                                                              >  2:5199265/1‑150 (MQ=255)
                         gagGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCAt                                                                                                                              <  1:793460/150‑1 (MQ=255)
                           ggTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATgg                                                                                                                            <  1:18093/150‑1 (MQ=255)
                             tGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCg                                                                                                                          >  2:801544/1‑150 (MQ=255)
                               cggcgTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCgtg                                                                                                                        <  1:6352395/150‑1 (MQ=255)
                                ggcgTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGACAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCgtgt                                                                                                                       >  1:3740722/1‑150 (MQ=255)
                                    tCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGGTGTCCGACGGCAGCCTGGCGGTAGAGGAGGTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCa                                                                                                                   <  1:3582598/150‑1 (MQ=255)
                                      aGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGt                                                                                                                 >  1:6869494/1‑150 (MQ=255)
                                        cTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTtg                                                                                                               <  2:3206509/150‑1 (MQ=255)
                                         tGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTtgg                                                                                                              <  2:6146943/150‑1 (MQ=255)
                                          gcgcGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTtggt                                                                                                             <  2:1705618/150‑1 (MQ=255)
                                           cgcgAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTtggtg                                                                                                            <  1:5933210/150‑1 (MQ=255)
                                            gcgAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTtggtgg                                                                                                           >  1:4072358/1‑150 (MQ=255)
                                              gaagaaGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGag                                                                                                         >  2:6121128/1‑150 (MQ=255)
                                              gaagaaGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGag                                                                                                         >  1:6424747/1‑150 (MQ=255)
                                              gaagaaGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGag                                                                                                         <  2:2292741/150‑1 (MQ=255)
                                               aagaagACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGagc                                                                                                        >  2:5734814/1‑150 (MQ=255)
                                                       cTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGtgtg                                                                                                >  1:1568945/1‑150 (MQ=255)
                                                       cTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGtgtg                                                                                                <  2:3055774/150‑1 (MQ=255)
                                                         gACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGt                                                                                              >  2:5114161/1‑150 (MQ=255)
                                                         gACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGt                                                                                              >  1:4619823/1‑150 (MQ=255)
                                                         gACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGt                                                                                              >  2:4818423/1‑150 (MQ=255)
                                                            tGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTgg                                                                                           >  2:2780829/1‑150 (MQ=255)
                                                            tGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTgg                                                                                           >  1:2468206/1‑150 (MQ=255)
                                                            tGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTgg                                                                                           >  2:134334/1‑150 (MQ=255)
                                                             ggCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGc                                                                                          <  1:4142741/150‑1 (MQ=255)
                                                             ggCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGc                                                                                          >  1:6575316/1‑150 (MQ=255)
                                                              gccgccTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCg                                                                                         >  1:1645627/1‑150 (MQ=255)
                                                               ccgccTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTGGGTGGAGCAGCGTGTGGTTGGCGa                                                                                        >  2:5602312/1‑150 (MQ=255)
                                                                cgccTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGaa                                                                                       <  2:3508580/150‑1 (MQ=255)
                                                                   cTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACgg                                                                                    <  2:3310534/150‑1 (MQ=255)
                                                                   cTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACCTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACgg                                                                                    >  1:4244844/1‑150 (MQ=255)
                                                                   cTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACATGGCCAATGCCCTGCGCCATGCCGGCCAATGGGGCCAGCACTTCCCTGAGCCGCTGTTCAATGGCGTGTTACAGTTGGTGGAGCAGCGTGTGGTTGGCGAACgg                                                                                    >  1:257382/1‑150 (MQ=255)
                                                                   cTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTATACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGACGCTGTTCCATGGAGTGTTCCAGTTGGTGGAGCAGAGTGTGGTTGGCGAACgg                                                                                    >  2:271229/1‑150 (MQ=255)
                                                                       tgtCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCAcc                                                                                >  1:1057905/1‑150 (MQ=255)
                                                                         tCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCAGGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCGGCGGGTGGTTGGCGAACGGCACCTg                                                                              >  1:2278750/1‑150 (MQ=255)
                                                                           cGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGaa                                                                            >  2:3138597/1‑150 (MQ=255)
                                                                            gACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGaaa                                                                           <  2:3877470/150‑1 (MQ=255)
                                                                            gACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGaaa                                                                           >  1:2712091/1‑150 (MQ=255)
                                                                            gACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGaaa                                                                           >  1:1096957/1‑150 (MQ=255)
                                                                                 cAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAgtggt                                                                      <  1:6796780/150‑1 (MQ=255)
                                                                                  aGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAgtggtg                                                                     <  1:2156435/150‑1 (MQ=255)
                                                                                   gCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGc                                                                    <  1:3024374/150‑1 (MQ=255)
                                                                                      tGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCa                                                                 >  2:6388602/1‑150 (MQ=255)
                                                                                      tGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCa                                                                 >  1:128657/1‑150 (MQ=255)
                                                                                      tGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGGTGGTGGAGCAGCGGGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCa                                                                 >  2:4152831/1‑150 (MQ=255)
                                                                                       ggCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCaa                                                                >  1:4707898/1‑150 (MQ=255)
                                                                                         cGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAaga                                                              <  1:6015531/150‑1 (MQ=255)
                                                                                           gTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGc                                                            <  2:5971375/150‑1 (MQ=255)
                                                                                           gTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGc                                                            <  1:4009142/150‑1 (MQ=255)
                                                                                               aggagTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAAt                                                        <  2:3700590/150‑1 (MQ=255)
                                                                                                 gagTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAAtgt                                                      <  2:6428672/150‑1 (MQ=255)
                                                                                                 gagTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAAtgt                                                      <  2:2928523/150‑1 (MQ=255)
                                                                                                 gagTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAAtgt                                                      <  2:1231758/150‑1 (MQ=255)
                                                                                                   gTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTgg                                                    <  2:592735/150‑1 (MQ=255)
                                                                                                   gTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTgg                                                    <  2:4961920/150‑1 (MQ=255)
                                                                                                         ccTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGt                                              >  2:3614204/1‑150 (MQ=255)
                                                                                                                  ggCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGa                                     <  1:339952/150‑1 (MQ=255)
                                                                                                                   gCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGAt                                    >  2:6512153/1‑150 (MQ=255)
                                                                                                                    ccAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGGTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCGGAAAGTGGTGCTCAAGTGCGAAGGTGGCGCGGTGCGGCTGGATg                                   >  2:617624/1‑150 (MQ=255)
                                                                                                                     cAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATgg                                  >  1:4799305/1‑150 (MQ=255)
                                                                                                                      aaGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGc                                 >  2:343910/1‑150 (MQ=255)
                                                                                                                      aaGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGc                                 <  1:5199265/150‑1 (MQ=255)
                                                                                                                        ggCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCAt                               <  1:3154533/150‑1 (MQ=255)
                                                                                                                         gCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCAtt                              >  2:923082/1‑150 (MQ=255)
                                                                                                                           ccTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGc                            >  1:4852544/1‑150 (MQ=255)
                                                                                                                            cTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTGGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGcc                           >  1:2884999/1‑150 (MQ=255)
                                                                                                                               cgcCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTc                        <  2:467803/150‑1 (MQ=255)
                                                                                                                                gcCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCg                       >  2:2263010/1‑150 (MQ=255)
                                                                                                                                 ccATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCgg                      <  1:986214/150‑1 (MQ=255)
                                                                                                                                    tGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCAt                   >  1:1029014/1‑150 (MQ=255)
                                                                                                                                    tGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCAt                   >  1:4514462/1‑150 (MQ=255)
                                                                                                                                     ctggccccATGGGCCGGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTg                 >  2:4465560/6‑150 (MQ=255)
                                                                                                                                      ccGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTg                 <  2:450792/150‑1 (MQ=255)
                                                                                                                                              tGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGaa         >  1:6237219/1‑150 (MQ=255)
                                                                                                                                                gggCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAgt       <  2:4898484/150‑1 (MQ=255)
                                                                                                                                                  gCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAgtgt     <  1:1710706/150‑1 (MQ=255)
                                                                                                                                                  gCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAgtgt     >  2:5383851/1‑150 (MQ=255)
                                                                                                                                                   ccAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAGtgtg    <  1:5945419/150‑1 (MQ=255)
                                                                                                                                                   ccAGCACGTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAGtgtg    <  1:4639723/150‑1 (MQ=255)
                                                                                                                                                    cAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAGTGTgg   >  2:5452924/1‑150 (MQ=255)
                                                                                                                                                     aGCACTACCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGACTGTGACTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAGTGTGGc  >  1:5657966/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAGTGTGGC  >  NC_002947/1687777‑1688075

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: