Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 725,276 G→A 5.2% Y587Y (TAC→TAT kdpB ← potassium translocating ATPase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913725,2760GA 5.2% 521.6 / 11.9 157Y587Y (TAC→TATkdpBpotassium translocating ATPase, subunit B
Reads supporting (aligned to +/- strand):  ref base G (65/83);  new base A (2/6);  total (68/89)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

Only 100 of 157 total aligned reads displayed.

AACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGT  >  NC_000913/725130‑725422
                                                                                                                                                  |                                                                                                                                                  
aaCGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTa                                                                                                                                                  >  3:1315709/1‑149 (MQ=255)
 aCGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTAt                                                                                                                                                 <  3:472038/149‑1 (MQ=255)
  cGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTAtt                                                                                                                                                <  3:277860/149‑1 (MQ=255)
    ccAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTc                                                                                                                                              <  1:1142854/149‑1 (MQ=255)
    ccAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTc                                                                                                                                              <  2:829266/149‑1 (MQ=255)
    ccAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTc                                                                                                                                              <  4:1286140/149‑1 (MQ=255)
     cAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCg                                                                                                                                             >  4:929719/1‑149 (MQ=255)
       ggggAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGcc                                                                                                                                           <  3:331387/149‑1 (MQ=255)
        gggAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAATACTTCGCca                                                                                                                                          <  4:623241/149‑1 (MQ=255)
         ggAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAATACTTCGCcat                                                                                                                                         <  4:733272/149‑2 (MQ=255)
         ggAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAATACTTCGCcat                                                                                                                                         <  2:405410/149‑2 (MQ=255)
             tCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCg                                                                                                                                     >  3:1424038/1‑149 (MQ=255)
             tCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCg                                                                                                                                     >  4:800404/1‑149 (MQ=255)
               aaaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGtt                                                                                                                                   >  3:723391/1‑149 (MQ=255)
               aaaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGtt                                                                                                                                   <  1:820601/149‑1 (MQ=255)
               aaaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGtt                                                                                                                                   <  1:517192/149‑1 (MQ=255)
                aaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTg                                                                                                                                  <  4:966490/149‑1 (MQ=255)
                aaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTg                                                                                                                                  <  3:1011004/149‑1 (MQ=255)
                aaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTg                                                                                                                                  <  2:269106/149‑1 (MQ=255)
                aaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTg                                                                                                                                  <  1:1151187/149‑1 (MQ=255)
                 aaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTgg                                                                                                                                 <  1:669314/149‑1 (MQ=255)
                 aaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTgg                                                                                                                                 <  3:203673/149‑1 (MQ=255)
                 aaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTgg                                                                                                                                 <  4:10282/149‑1 (MQ=255)
                 aaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTgg                                                                                                                                 <  2:334571/149‑1 (MQ=255)
                 aaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTgg                                                                                                                                 >  1:1460255/1‑149 (MQ=255)
                  aaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCca                                                                                                                                          <  2:1284554/139‑1 (MQ=255)
                    gACGATCATCAAGGCTTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCaa                                                                                                                              <  2:737097/149‑1 (MQ=255)
                       gATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGc                                                                                                                           >  1:1326414/1‑149 (MQ=255)
                       gATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGc                                                                                                                           >  4:325875/1‑149 (MQ=255)
                          aaTCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGa                                                                                                                        >  2:923329/1‑149 (MQ=255)
                            tCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTg                                                                                                                                  >  4:1100364/1‑137 (MQ=255)
                               aGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAggtg                                                                                                                   <  3:1170324/149‑1 (MQ=255)
                                   gTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCa                                                                                                               <  3:725790/149‑1 (MQ=255)
                                       aaaaTCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACAGCGTTGGCAATGCTGAAGGTGGTCAGCGa                                                                                                           <  3:97984/149‑1 (MQ=255)
                                        aaaTCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAg                                                                                                          <  1:923329/149‑1 (MQ=255)
                                         aaTCACCGCACTGAGGATTGCGGAGTCGGGTGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGc                                                                                                         >  4:480859/1‑149 (MQ=255)
                                           tCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGc                                                                                                                                <  4:386412/124‑1 (MQ=255)
                                              ccGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGTATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACgg                                                                                                    <  3:896812/149‑1 (MQ=255)
                                                gCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGt                                                                                                  <  1:473241/149‑1 (MQ=255)
                                                 cACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGtc                                                                                                 <  2:216290/149‑1 (MQ=255)
                                                    tGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGtcatc                                                                                              <  2:1326414/149‑1 (MQ=255)
                                                     gAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGtcatca                                                                                             <  1:503192/149‑1 (MQ=255)
                                                         aTTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCAt                                                                                         >  2:1416882/1‑149 (MQ=255)
                                                             cGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGt                                                                                     <  3:815938/149‑1 (MQ=255)
                                                             cGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGt                                                                                     <  4:1297741/149‑1 (MQ=255)
                                                              ggAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTgg                                                                                                                                 >  3:365232/1‑104 (MQ=255)
                                                              ggAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGtt                                                                                    <  4:925034/149‑1 (MQ=255)
                                                              ggAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAATACTTCGCcat                                                                                                                                         <  3:624853/96‑2 (MQ=255)
                                                               gAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGttt                                                                                   <  3:1401344/149‑1 (MQ=255)
                                                                 gTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGc                                                                                 <  2:448790/149‑1 (MQ=255)
                                                                 gTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGc                                                                                 <  2:1244326/149‑1 (MQ=255)
                                                                    ggacGAATGCATCCATATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAAAGCCGCCGGAAGAATGGCGAAGTATTTCGCCACATCCGTGGCAATGCTGAAGGTGGTCAACGAGCCACAGATCATCAGCATCTGTTTGCCaa                                                                              >  4:274586/4‑149 (MQ=255)
                                                                       cGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAAtgt                                                                           <  3:1213296/149‑1 (MQ=255)
                                                                       cGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAAtgt                                                                           <  3:553541/149‑1 (MQ=255)
                                                                          aTGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGca                                                                        >  2:199132/1‑149 (MQ=255)
                                                                          aTGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGca                                                                        >  1:1102551/1‑149 (MQ=255)
                                                                              aGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGcaccac                                                                    >  1:1357598/1‑149 (MQ=255)
                                                                                gCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCt                                                                  >  4:49688/1‑149 (MQ=255)
                                                                                gCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCt                                                                  >  4:359833/1‑149 (MQ=255)
                                                                                     tgatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATc                                                             >  1:145723/1‑149 (MQ=255)
                                                                                           tCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTg                                                       <  3:1022980/149‑1 (MQ=255)
                                                                                             aGCGCCTTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGt                                                     >  1:400712/1‑149 (MQ=255)
                                                                                              gcgcATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACgg                                                                                                    >  3:943870/1‑101 (MQ=255)
                                                                                              gcgcATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTc                                                    <  1:121178/149‑1 (MQ=255)
                                                                                              gcgcATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTc                                                    >  3:151663/1‑149 (MQ=255)
                                                                                              gcgcATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTc                                                    <  1:275804/149‑1 (MQ=255)
                                                                                              gcgcATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAAGAATGGCGAACTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTc                                                    >  4:339701/1‑149 (MQ=255)
                                                                                                  aTTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAAc                                                          <  2:1314548/139‑1 (MQ=255)
                                                                                                  aTTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGt                                                <  1:543802/149‑1 (MQ=255)
                                                                                                  aTTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGt                                                >  4:264284/1‑149 (MQ=255)
                                                                                                   tttAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGtt                                               >  2:469587/1‑149 (MQ=255)
                                                                                                      aaCTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTaga                                            <  1:1386104/149‑1 (MQ=255)
                                                                                                       aCTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTagag                                           >  1:690547/1‑149 (MQ=255)
                                                                                                        cTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGt                                          >  4:603592/1‑149 (MQ=255)
                                                                                                        cTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGt                                          <  1:1309366/149‑1 (MQ=255)
                                                                                                             gATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCgaga                                     <  1:199132/149‑1 (MQ=255)
                                                                                                               tACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGt                                          <  3:665379/142‑1 (MQ=255)
                                                                                                               tACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGt                                          >  4:665379/1‑142 (MQ=255)
                                                                                                                  gTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGAc                                <  2:1369402/149‑1 (MQ=255)
                                                                                                                  gTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGAc                                <  3:459728/149‑1 (MQ=255)
                                                                                                                   ttGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTc                                         <  4:475278/139‑1 (MQ=255)
                                                                                                                      ccGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCt                                                                  >  2:1025050/1‑111 (MQ=255)
                                                                                                                      ccGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCt                                                                  <  1:1025050/111‑1 (MQ=255)
                                                                                                                       cgcgAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTg                                                                                  >  4:1249457/1‑94 (MQ=255)
                                                                                                                           aaTGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAAc                                                          >  4:1375628/1‑114 (MQ=255)
                                                                                                                           aaTGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAAc                                                          <  3:1375628/114‑1 (MQ=255)
                                                                                                                                 gccgGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGc                        >  3:1027462/1‑142 (MQ=255)
                                                                                                                                 gccgGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCctc                 >  4:54760/1‑149 (MQ=255)
                                                                                                                                   cgGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCtt               <  3:264284/149‑1 (MQ=255)
                                                                                                                                    ggAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCttt              <  4:723391/149‑1 (MQ=255)
                                                                                                                                     gAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTc             <  3:800404/149‑1 (MQ=255)
                                                                                                                                     gAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTc             >  2:767068/1‑149 (MQ=255)
                                                                                                                                       ataatGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTcgc           >  1:835943/1‑149 (MQ=255)
                                                                                                                                          atGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTcgccgc        >  2:1111550/1‑149 (MQ=255)
                                                                                                                                            ggCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCt      <  1:1300451/149‑1 (MQ=255)
                                                                                                                                               gAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTc                                                    >  2:17555/1‑100 (MQ=255)
                                                                                                                                               gAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTggg   >  1:1240578/1‑149 (MQ=255)
                                                                                                                                               gAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTggg   >  3:376233/1‑149 (MQ=255)
                                                                                                                                               gAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTggg   >  4:232559/1‑149 (MQ=255)
                                                                                                                                                aaGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGt  <  3:213661/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGT  >  NC_000913/725130‑725422

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: