Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,295,024 T→A 11.8% intergenic (+34/+36) slyB → / ← nhaA‑I outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,295,0240TA11.8% 289.0 / 26.4 102intergenic (+34/+36)slyB/nhaA‑Iouter membrane lipoprotein/Na(+)/H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (43/47);  new base A (7/5);  total (50/52)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

Only 100 of 110 total aligned reads displayed.

AAGATGATGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGC  >  NC_002947/1294875‑1295172
                                                                                                                                                     |                                                                                                                                                     
      atgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGcc                                                                                                                                                 <  2:287001/150‑1 (MQ=255)
          cagcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTc                                                                                                                                            <  1:1158276/150‑1 (MQ=255)
          cagcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCatgtc                                                                                                                                            <  2:2277285/150‑5 (MQ=255)
           agcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCg                                                                                                                                           >  2:574202/1‑150 (MQ=255)
           agcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCatgtcg                                                                                                                                           >  2:911596/1‑145 (MQ=255)
             caTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCggg                                                                                                                                         <  1:1587082/150‑1 (MQ=255)
               tGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGt                                                                                                                                       <  2:83880/150‑1 (MQ=255)
               tGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCATGTCGGGGt                                                                                                                                       <  1:911596/150‑1 (MQ=255)
                  gcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGGAGCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGtttt                                                                                                                                    >  1:2126712/1‑150 (MQ=255)
                  gcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCATGTCGGGGtttt                                                                                                                                    >  2:893108/1‑150 (MQ=39)
                   cgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGttttt                                                                                                                                   <  2:752078/150‑1 (MQ=255)
                    gcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGCGTGGGTGACCGCGTGCGCATCATGACCGTTGACGGTACCAGCCGCGTTTCGCACTGACCGTAGCGCAGTAAAACAAAACCCCGACCTGGCTGCTGCCAGTCGCCGttcttt                                                                                                                                  <  2:1543184/150‑1 (MQ=255)
                      cTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCa                                                                                                                                >  2:2249934/1‑150 (MQ=255)
                      cTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCAAACATTGCTGCTGCCAGTCGGGGTTTTTTCa                                                                                                                                >  2:302121/1‑150 (MQ=255)
                      cTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGAACGTCGTCCAGCAATAACAACCTCCCATATGCCTCCTGCCAGCCGGGGTCTTTTCa                                                                                                                                >  2:1453902/1‑150 (MQ=255)
                        aTGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCATGTCGGGGttttttatta                                                                                                                              <  1:893108/150‑5 (MQ=38)
                         tgtgCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGc                                                                                                                             <  1:1289842/150‑1 (MQ=255)
                                      cAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCa                                                                                                                <  1:2249934/150‑1 (MQ=255)
                                        ggAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAgg                                                                                                              >  1:1397255/1‑150 (MQ=255)
                                           gaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCaa                                                                                                           >  1:1619512/1‑150 (MQ=255)
                                            aaCGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAAt                                                                                                          >  1:896594/1‑150 (MQ=255)
                                              cGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGc                                                                                                        >  2:1158460/1‑150 (MQ=255)
                                               gAAATCTTCCGTGTTGGCGAACGCGCACGGATTCTGGCTGTCGGCCGTTCCAGGTGCGGTTCGCAGTTAGCGTTGTGCAGGTATAAAAAACCCCGGCATGGCTGCTGCCAGTCGGGGTTTTTTTATGCCTCAGCCCAGTCATGCAATGCt                                                                                                       >  2:1466840/1‑150 (MQ=255)
                                                aaaTCTTCCGTGTTGGCGACCGAGTGCGCATCATGATTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTg                                                                                                      <  1:1957640/150‑1 (MQ=255)
                                                aaaTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTg                                                                                                      >  2:796318/1‑150 (MQ=255)
                                                  aTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGt                                                                                                    <  1:941202/150‑1 (MQ=255)
                                                  aTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGCGTTTGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGCCATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGt                                                                                                    <  1:1453902/150‑1 (MQ=255)
                                                   tCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGtt                                                                                                   <  1:1552760/150‑1 (MQ=255)
                                                     ttCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCATGTCggg                                                                                                                                         <  2:760448/110‑1 (MQ=38)
                                                     ttCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCATGTCggg                                                                                                                                         >  1:760448/1‑110 (MQ=38)
                                                      tCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTg                                                                                                <  1:796318/150‑1 (MQ=255)
                                                       ccGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGc                                                                                               <  2:395182/150‑1 (MQ=255)
                                                            ttGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTg                                                                                          >  1:815704/1‑150 (MQ=255)
                                                                 gACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGgctgct                                                                                                                                                    >  2:893191/1‑88 (MQ=255)
                                                                 gACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGgctgct                                                                                                                                                    <  1:893191/88‑1 (MQ=255)
                                                                     gAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCg                                                                                 <  1:1873069/150‑1 (MQ=255)
                                                                     gAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCg                                                                                 <  2:1619512/150‑1 (MQ=255)
                                                                         gcgcGTCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGa                                                                             <  2:410588/150‑1 (MQ=255)
                                                                          cgcATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGAc                                                                            <  1:1531950/150‑1 (MQ=255)
                                                                           gcATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAAGAGGTTTCTCCGCCATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACg                                                                           <  2:324226/150‑1 (MQ=255)
                                                                             atcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGc                                                                         >  1:1904805/1‑150 (MQ=255)
                                                                              tcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTgg                                                                                         <  1:1211521/133‑1 (MQ=255)
                                                                              tcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTgg                                                                                         >  2:1211521/1‑133 (MQ=255)
                                                                                atGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTa                                                                      >  2:284666/1‑150 (MQ=255)
                                                                                    cTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTTGTTTTGCTGCTGGCCATTGCGGTGAAGGCATAggac                                                                  >  2:2230845/1‑147 (MQ=255)
                                                                                     tGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGCCGGGGTTTTTTCATGCCTCAGCCAAGTAAGGGAGTACGTTATTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGa                                                                 >  1:2377122/1‑150 (MQ=255)
                                                                                      gTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGAt                                                                <  1:574202/150‑1 (MQ=255)
                                                                                        aGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCa                                                              <  1:1177640/150‑1 (MQ=255)
                                                                                          aCGNTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAcc                                                            >  2:1514233/1‑150 (MQ=255)
                                                                                          aCGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAcc                                                            >  1:374699/1‑150 (MQ=255)
                                                                                          aCGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAcc                                                            <  2:374699/150‑1 (MQ=255)
                                                                                            ggTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTTCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgc                                                          <  1:627140/150‑1 (MQ=255)
                                                                                             gTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgca                                                         <  1:775715/150‑1 (MQ=255)
                                                                                                  aGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGAGACGGCGTAGCCGATCACCGCAGCaaa                                                    >  2:1270703/1‑150 (MQ=255)
                                                                                                    tCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTTACGGCGTAGCCGATCACCGCAGCAAAga                                                  >  1:1318293/1‑150 (MQ=255)
                                                                                                    tCGCGTTTCGCACTGACCGTAGCGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTg                                                                                                      >  2:1489680/1‑98 (MQ=255)
                                                                                                    tCGCGTTTCGCACTGACCGTAGCGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTg                                                                                                      <  1:1489680/98‑1 (MQ=255)
                                                                                                       cgTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacg                                               >  1:1096877/1‑150 (MQ=255)
                                                                                                         tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaa                                             <  1:1372222/150‑1 (MQ=255)
                                                                                                         tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaa                                                 <  2:2118223/146‑1 (MQ=255)
                                                                                                         tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaa                                                 >  1:2118200/1‑146 (MQ=255)
                                                                                                         tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaa                                                 >  1:2118223/1‑146 (MQ=255)
                                                                                                          ttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGa                                                                             <  1:2248475/117‑1 (MQ=255)
                                                                                                          ttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGa                                                                             >  2:2248475/1‑117 (MQ=255)
                                                                                                            cGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGTCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCa                                          <  2:1374408/150‑1 (MQ=255)
                                                                                                              cACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAgg                                        >  1:459437/1‑149 (MQ=255)
                                                                                                              cACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGt                                        >  1:2063006/1‑150 (MQ=255)
                                                                                                                  gACCGTAGTGCAGTAATAAAAAACCCCGACATGGCtgctgc                                                                                                                                                  <  2:1338544/41‑1 (MQ=255)
                                                                                                                  gACCGTAGTGCAGTAATAAAAAACCCCGACATGGCtgctgc                                                                                                                                                  >  1:1338544/1‑41 (MQ=255)
                                                                                                                    ccGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAAt                                  <  1:284666/150‑1 (MQ=255)
                                                                                                                     cGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCgg                                                                                >  2:2321653/1‑103 (MQ=255)
                                                                                                                     cGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCgg                                                                                <  1:2321653/103‑1 (MQ=255)
                                                                                                                          tGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCAt                            <  2:435339/150‑1 (MQ=255)
                                                                                                                            cAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTc                                                                                                                                            <  2:720764/36‑1 (MQ=255)
                                                                                                                            cAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTc                                                                                                                                            >  1:720764/1‑36 (MQ=255)
                                                                                                                            cAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGc                          <  1:1466840/150‑1 (MQ=255)
                                                                                                                             aGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAg                                         >  1:2407233/1‑134 (MQ=255)
                                                                                                                             aGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCg                         <  1:1834285/150‑1 (MQ=255)
                                                                                                                             aGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTCGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAg                                         <  2:2407233/134‑1 (MQ=255)
                                                                                                                              gTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGATGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCgg                        <  1:629918/150‑1 (MQ=255)
                                                                                                                               taataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTATTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCCTGCggc                       >  2:1697666/1‑149 (MQ=255)
                                                                                                                               aaataaAAAACCCCGACATGGCTGCTGTCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGAAGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGt                       <  1:2014573/149‑1 (MQ=255)
                                                                                                                                aataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTc                      >  2:304573/1‑150 (MQ=255)
                                                                                                                                 ataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGc                                                                    <  1:1535025/103‑1 (MQ=255)
                                                                                                                                 ataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGc                                                                    >  2:1535025/1‑103 (MQ=255)
                                                                                                                                    aaaaaCCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCCCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATg                  <  1:302121/150‑1 (MQ=255)
                                                                                                                                       aaCCCCGACTGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCg               >  1:385675/1‑150 (MQ=255)
                                                                                                                                        aCCCCGACTGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                            >  2:562045/1‑120 (MQ=38)
                                                                                                                                        aCCCCGACTGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                            <  1:562045/120‑1 (MQ=38)
                                                                                                                                         ccccGACTGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc             >  2:2338918/1‑150 (MQ=255)
                                                                                                                                          cccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGttt                                                                                                  >  2:117656/1‑64 (MQ=255)
                                                                                                                                          cccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGttt                                                                                                  <  1:117656/64‑1 (MQ=255)
                                                                                                                                           ccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCgg                                                                                >  2:138460/1‑81 (MQ=255)
                                                                                                                                           ccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCgg                                                                                <  1:138460/81‑1 (MQ=255)
                                                                                                                                             gACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTAc         >  1:1956421/1‑150 (MQ=255)
                                                                                                                                              aCATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACt        <  1:1158460/150‑1 (MQ=255)
                                                                                                                                                 tGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCAc     <  2:534055/150‑1 (MQ=255)
                                                                                                                                                   gctgctGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGcc                               >  2:2309076/1‑122 (MQ=255)
                                                                                                                                                   gctgctGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGcc                               <  1:2309076/122‑1 (MQ=255)
                                                                                                                                                    ctgctgCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGc  <  2:2377122/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGATGATGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGC  >  NC_002947/1294875‑1295172

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