Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,707,321 A→T 7.3% intergenic (‑197/+1) lepA ← / ← rseC back‑translocating elongation factor EF4, GTPase/SoxR iron‑sulfur cluster reduction factor component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,707,3210AT 7.3% 287.7 / 9.9 96intergenic (‑197/+1)lepA/rseCback‑translocating elongation factor EF4, GTPase/SoxR iron‑sulfur cluster reduction factor component
Reads supporting (aligned to +/- strand):  ref base A (52/36);  new base T (2/5);  total (55/41)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

GTCGCTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGG  >  NC_000913/2707173‑2707469
                                                                                                                                                    |                                                                                                                                                    
gTCGCTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGaaa                                                                                                                                                      <  2:557984/149‑1 (MQ=255)
    cTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCAc                                                                                                                                                  <  2:1579619/149‑1 (MQ=255)
    cTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCAc                                                                                                                                                  <  2:1585674/149‑1 (MQ=255)
     tGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACt                                                                                                                                                 >  1:246009/1‑149 (MQ=255)
     tGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACt                                                                                                                                                 >  2:1534528/1‑149 (MQ=255)
     tGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACt                                                                                                                                                 <  1:793677/149‑1 (MQ=255)
      gTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTg                                                                                                                                                <  2:1575899/149‑1 (MQ=255)
          tGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtt                                                                                                                                            <  1:651396/149‑2 (MQ=255)
           gAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTcg                                                                                                                                           >  1:480197/1‑149 (MQ=255)
            agattAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCATAATTCAACAAAGCCAGTGAAATCACTGGCTcgc                                                                                                                                          >  2:1277430/4‑149 (MQ=255)
            aGCTTAACGTCTGTATTAATAGAAAAGCCGCATTCTACTCTCCAACATTGAGGCGAGGAATTGTTCATACCGTATGGATTGTGGTTTCTGGTATCGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttg                                                                                                                                          <  1:789587/149‑4 (MQ=255)
             gCTTACCTTCTTTATTAATAGAACCGGCGCTTTCTACACTACTAAATTGAGGCGAGGAAATGTTCATACCGTTTGGTTTGTGGTATCTTTATAGGTCCTCGCATTTGTTATGAAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgt                                                                                                                                         <  2:53724/149‑5 (MQ=255)
               ttAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgttg                                                                                                                                       <  1:308995/149‑7 (MQ=255)
                tAACGTCTGTATTAATAGTAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCTTATGGTTTGTGGTATCTGTAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTTAAATCACTGGCTCGCGTCt                                                                                                                                      <  2:398641/149‑1 (MQ=255)
                 aaCGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCtt                                                                                                                                     >  1:1402126/1‑149 (MQ=255)
                 aaCGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCtt                                                                                                                                     >  2:1527615/1‑149 (MQ=255)
                 aaCGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGCTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCtt                                                                                                                                     >  2:442077/1‑149 (MQ=255)
                      cTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGaa                                                                                                                                <  1:525057/149‑1 (MQ=255)
                       tGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAg                                                                                                                               >  1:1376513/1‑149 (MQ=255)
                       tGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAg                                                                                                                               >  2:1295089/1‑149 (MQ=255)
                        gTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGTAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGa                                                                                                                              >  2:866816/1‑149 (MQ=255)
                         tATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGAt                                                                                                                             >  2:1584791/1‑149 (MQ=255)
                         tATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGAt                                                                                                                             >  2:1145595/1‑149 (MQ=255)
                           ttAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGt                                                                                                                           <  2:422417/149‑1 (MQ=255)
                             aaTAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGttt                                                                                                                         <  1:375634/149‑1 (MQ=255)
                              aTAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTAGCGTCTTCCGAAGATGTTTc                                                                                                                        >  2:686038/1‑149 (MQ=255)
                               tAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCa                                                                                                                       >  2:421263/1‑149 (MQ=255)
                                 gAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCTTACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCaaa                                                                                                                     >  2:304407/1‑149 (MQ=255)
                                 gAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCaaa                                                                                                                     >  2:1406179/1‑149 (MQ=255)
                                  aaaCGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAAt                                                                                                                    <  2:126102/149‑1 (MQ=255)
                                  aaaCGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAAt                                                                                                                    <  2:1523682/149‑1 (MQ=255)
                                       ccgcATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGcac                                                                                                               >  2:1100912/1‑149 (MQ=255)
                                       ccgcATTCTACACAACAACATTGAGGCGAGGAAATGTTCATACCGTATCGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAACTCACTGGCTCGCGTCTTCAGAAGATGTTTCAAATCGcac                                                                                                               >  1:1369321/1‑149 (MQ=255)
                                               tacactacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGATATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAg                                                                                                       >  1:609593/1‑149 (MQ=255)
                                                acactacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAgg                                                                                                      <  2:1136876/149‑1 (MQ=255)
                                                 cactacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAggc                                                                                                     <  2:654643/149‑1 (MQ=255)
                                                   ctacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAggcgg                                                                                                   <  1:1295089/149‑1 (MQ=255)
                                                   ctacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAggcgg                                                                                                   <  2:1132441/149‑1 (MQ=255)
                                                    tacaACATAGAGGCGAGGAAAGGTTCATACCGTATGGTTTGTGGTATCTGGTAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGAAGTTTAAAATCGCAAATGGCCTggagga                                                                                                  >  1:802942/1‑148 (MQ=255)
                                                        acaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATg                                                                                              <  1:33788/149‑1 (MQ=255)
                                                         caTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGc                                                                                             <  2:985914/149‑1 (MQ=255)
                                                           ttGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGCAAACGTCCTCGCATTTGTTATGCAAAATGAAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTATCAAATCCCAGAAGGCCATGCGGCAATGCCa                                                                                           >  2:1444059/1‑149 (MQ=255)
                                                            tGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCAc                                                                                          <  2:596735/149‑1 (MQ=255)
                                                             gAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACg                                                                                         <  1:813106/149‑1 (MQ=255)
                                                              aGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGc                                                                                        <  1:1594483/149‑1 (MQ=255)
                                                                gCGAGGATATGTTCTTACCGTTTGGATTGTGGTCTCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAAGCACTGTCTCGCGTCTTCCGAAGATGTTTCAAATCTCACAAGGCCAGGCGGCAATGCCACGCtt                                                                                      <  2:54724/149‑1 (MQ=255)
                                                                gCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCCCGCtt                                                                                      >  1:284338/1‑149 (MQ=255)
                                                                gCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCtt                                                                                      >  1:552211/1‑149 (MQ=255)
                                                                 cGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTa                                                                                     >  1:1325499/1‑149 (MQ=255)
                                                                 cGAGGAAATGTTCATACCGGATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTa                                                                                     <  1:1018497/149‑1 (MQ=255)
                                                                   aGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTGCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAg                                                                                   <  1:1265752/149‑1 (MQ=255)
                                                                    ggAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAga                                                                                  <  2:246009/149‑1 (MQ=255)
                                                                     aaaaTGTACATACCGTATGGATTGTGGTATCTGGAAACGTCCTCCAAATTGTTATGCAAAAAGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGACGGCAAAGCCACCATAAAgat                                                                                 >  2:98147/2‑149 (MQ=255)
                                                                       aaTGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgatga                                                                               <  2:671590/149‑1 (MQ=255)
                                                                        aTGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgatgat                                                                              <  1:1100912/149‑1 (MQ=255)
                                                                        aTGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgatgat                                                                              <  1:470470/149‑1 (MQ=255)
                                                                           ttCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCCCGCTTAAGATGATCgg                                                                           >  2:341085/1‑149 (MQ=255)
                                                                           ttCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCgg                                                                           >  1:1105893/1‑149 (MQ=255)
                                                                           ttCATACCGAATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCgg                                                                           <  2:583403/149‑1 (MQ=255)
                                                                             cATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGtt                                                                         >  2:1195711/1‑149 (MQ=255)
                                                                             cATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGtt                                                                         >  1:1388419/1‑149 (MQ=255)
                                                                             cATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGtt                                                                         >  1:904241/1‑149 (MQ=255)
                                                                                   gTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTGCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGACACGCTTAAGATGATCGGTTGCCAtt                                                                   >  2:680981/1‑149 (MQ=255)
                                                                                        gATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGcc                                                              >  2:417508/1‑149 (MQ=255)
                                                                                         aTTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGccc                                                             >  1:490004/1‑149 (MQ=255)
                                                                                           tgtgGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCgg                                                           >  1:467004/1‑149 (MQ=255)
                                                                                            gtgGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCggg                                                          >  2:1363198/1‑149 (MQ=255)
                                                                                            gtgGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCggg                                                          <  1:882014/149‑1 (MQ=255)
                                                                                              ggTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACACGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGGTGCCATACTGCCCGGgca                                                        >  2:1028824/1‑149 (MQ=255)
                                                                                               gTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgcag                                                       >  1:1331145/1‑149 (MQ=255)
                                                                                               gTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTACGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgcag                                                       >  1:284733/1‑149 (MQ=255)
                                                                                               gTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGACAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgcag                                                       <  2:975992/149‑1 (MQ=255)
                                                                                                   cTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCaaa                                                   <  1:304407/149‑1 (MQ=255)
                                                                                                    tGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACAGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCAAAAGGCCAGGAGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAAc                                                  >  1:530177/1‑149 (MQ=255)
                                                                                                      gAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACtt                                                >  2:850602/1‑149 (MQ=255)
                                                                                                         aCGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTAcg                                             <  1:1584791/149‑1 (MQ=255)
                                                                                                          cGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTAcgc                                            <  2:1390839/149‑1 (MQ=255)
                                                                                                          cGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTAcgc                                            <  1:1145595/149‑1 (MQ=255)
                                                                                                              cTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCCCGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGt                                        >  2:1265334/1‑149 (MQ=255)
                                                                                                               tCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTa                                       >  1:993451/1‑149 (MQ=255)
                                                                                                               tCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAACGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTa                                       >  2:1073687/1‑149 (MQ=255)
                                                                                                                 gCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCCCTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCAATTCTGCCCGGGCAGCAACCTTACGCGAGTAGc                                     >  1:1287350/1‑149 (MQ=255)
                                                                                                                  cATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGcc                                    <  2:480197/149‑1 (MQ=255)
                                                                                                                       gTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAAAGCCACGCATAAGATGATCGGTTGCCATTCTGCCAGGGCAGCAAACTTACGCGAGTAGCagcgcg                               >  2:78532/1‑149 (MQ=255)
                                                                                                                           tGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAATCTTACGCGAGTAGCCGCGCGCAAt                           <  1:866816/149‑1 (MQ=255)
                                                                                                                            gaaaaaTGCAACAAAGCAAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGTGCAGCAAACTTACGCGAGTAGCCGCGCGCAAta                          >  2:896704/3‑148 (MQ=255)
                                                                                                                              aaaaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAg                        >  1:843707/1‑149 (MQ=255)
                                                                                                                                 aTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGaa                     <  1:1534528/149‑1 (MQ=255)
                                                                                                                                    caacaaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGCTGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAcccc                  >  1:949842/1‑149 (MQ=255)
                                                                                                                                    caacaaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAAccc                  >  2:101719/1‑149 (MQ=255)
                                                                                                                                     aacaaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCg                 >  2:136499/1‑149 (MQ=255)
                                                                                                                                     aacaaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGCCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAAcccc                 >  1:1571606/1‑148 (MQ=255)
                                                                                                                                       gcgaGCCAGTGATTTCACTGGCTTGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAAt                           >  1:1265766/4‑137 (MQ=255)
                                                                                                                                       gcgaGCCAGTGATTTCACTGGCTTGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAAt                           <  2:1265766/134‑1 (MQ=255)
                                                                                                                                           gCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTATGCCCGGGCAGCAAACTTACGCGAGTAGCCGAGCGCAATCAGGAACCCGCCAATa           >  2:1404817/1‑149 (MQ=255)
                                                                                                                                            ccAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATac          >  1:1187143/1‑149 (MQ=255)
                                                                                                                                               gTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATaccac       >  1:490232/1‑149 (MQ=255)
                                                                                                                                                  aaaTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGa    >  1:548282/1‑149 (MQ=255)
                                                                                                                                                    aTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCCTTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAgg  >  1:542214/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCGCTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGG  >  NC_000913/2707173‑2707469

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