Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 725,276 G→A 6.9% 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 6.9% 276.4 / 8.3 87Y587Y (TAC→TAT) kdpBpotassium translocating ATPase, subunit B
Reads supporting (aligned to +/- strand):  ref base G (42/39);  new base A (3/3);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

TAACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGC  >  NC_000913/725129‑725424
                                                                                                                                                   |                                                                                                                                                    
tAACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGt                                                                                                                                                     <  3:598391/149‑1 (MQ=255)
 aaCGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTa                                                                                                                                                    <  3:554779/149‑1 (MQ=255)
    gCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTAttt                                                                                                                                                 >  2:163862/1‑149 (MQ=255)
       aGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGc                                                                                                                                              <  1:163862/149‑1 (MQ=255)
       aGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGc                                                                                                                                              <  4:166251/149‑1 (MQ=255)
          ggAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAATACTTCGCcat                                                                                                                                           <  1:204403/149‑2 (MQ=255)
              tCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCg                                                                                                                                       >  1:53174/1‑149 (MQ=255)
                aaaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGtt                                                                                                                                     >  2:531300/1‑149 (MQ=255)
                aaaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGtt                                                                                                                                     >  1:557751/1‑149 (MQ=255)
                 aaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTg                                                                                                                                    >  1:158391/1‑149 (MQ=255)
                 aaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTg                                                                                                                                    <  1:349468/149‑1 (MQ=255)
                   aaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGc                                                                                                                                  >  2:297913/1‑149 (MQ=255)
                     gACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCaa                                                                                                                                >  2:450905/1‑149 (MQ=255)
                        gATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGc                                                                                                                             >  1:42444/1‑149 (MQ=255)
                        gATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGACGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGc                                                                                                                             >  3:148015/1‑149 (MQ=255)
                            aTCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGaa                                                                                                                         <  2:158391/149‑1 (MQ=255)
                             tCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAg                                                                                                                        >  3:102033/1‑149 (MQ=255)
                               aaGGCGTTGAAAATCACCGCACTGTGGACTGCGGAGTCGGGCGAATGCACGCACATGATGTTAATCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGCAggt                                                                                                                      >  2:643319/1‑149 (MQ=255)
                               aaGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAggt                                                                                                                      >  3:201657/1‑149 (MQ=255)
                                 ggCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGtt                                                                                                                                     >  4:32389/1‑132 (MQ=255)
                                 ggCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGtt                                                                                                                                     <  3:32389/132‑1 (MQ=255)
                                     ttGAAAATCACAGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAg                                                                                                                >  1:558307/1‑149 (MQ=255)
                                       gAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCg                                                                                                              >  2:191104/1‑149 (MQ=255)
                                        aaaaTCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGa                                                                                                             <  3:78786/149‑1 (MQ=255)
                                           aTCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTGGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGcc                                                                                                          >  1:523186/1‑149 (MQ=255)
                                              aCCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACg                                                                                                       >  1:652799/1‑149 (MQ=255)
                                               ccGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACgg                                                                                                      <  1:572187/149‑1 (MQ=255)
                                               ccGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACgg                                                                                                      <  2:101633/149‑1 (MQ=255)
                                               ccGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACgg                                                                                                      <  1:188841/149‑1 (MQ=255)
                                                   aCTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGtca                                                                                                  <  3:78807/149‑1 (MQ=255)
                                                     tGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGtcatc                                                                                                >  4:394618/1‑149 (MQ=255)
                                                          aTTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGACAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCAt                                                                                           <  4:401298/149‑1 (MQ=255)
                                                           ttGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATc                                                                                          >  1:602764/1‑149 (MQ=255)
                                                           ttGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATc                                                                                          >  1:195186/1‑149 (MQ=255)
                                                           ttGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATc                                                                                          >  2:87633/1‑149 (MQ=255)
                                                            tGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCt                                                                                         <  1:297913/149‑1 (MQ=255)
                                                                gAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGttt                                                                                     <  2:53174/149‑1 (MQ=255)
                                                                     gggCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCaa                                                                                >  4:323903/1‑149 (MQ=255)
                                                                      ggCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAAt                                                                               <  4:247889/149‑1 (MQ=255)
                                                                      ggCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAAt                                                                               <  1:610069/149‑1 (MQ=255)
                                                                       gCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAAtg                                                                              <  1:87633/149‑1 (MQ=255)
                                                                        cGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAAtgt                                                                             <  3:394618/149‑1 (MQ=255)
                                                                          aaTGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGc                                                                           <  1:690278/149‑1 (MQ=255)
                                                                           aTGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGca                                                                          >  4:247720/1‑149 (MQ=255)
                                                                            tGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGcac                                                                         <  4:63779/149‑1 (MQ=255)
                                                                               aGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGcaccac                                                                      >  4:123055/1‑149 (MQ=255)
                                                                                   acaTGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCg                                                                  >  1:546388/1‑149 (MQ=255)
                                                                                   acaTGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCg                                                                  >  2:45564/1‑149 (MQ=255)
                                                                                    caTGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGa                                                                 <  2:195186/149‑1 (MQ=255)
                                                                                     atgatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAGTGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGAt                                                                >  4:87032/1‑149 (MQ=255)
                                                                                     atgatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGTTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGAt                                                                >  1:247527/1‑149 (MQ=255)
                                                                                     atgatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGAt                                                                <  4:102033/149‑1 (MQ=255)
                                                                                      tgatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGc                                                                           >  1:321349/1‑137 (MQ=255)
                                                                                      tgatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGc                                                                           >  3:479330/1‑137 (MQ=255)
                                                                                      tgatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGc                                                                           <  2:321349/137‑1 (MQ=255)
                                                                                      tgatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGc                                                                           <  4:479330/137‑1 (MQ=255)
                                                                                       gatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCa                                                              >  4:454206/1‑149 (MQ=255)
                                                                                       gatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCa                                                              >  4:73258/1‑149 (MQ=255)
                                                                                        atgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCaa                                                             <  2:669899/149‑1 (MQ=255)
                                                                                        atgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCaa                                                             >  4:76138/1‑149 (MQ=255)
                                                                                         tgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAAc                                                            <  2:42444/149‑1 (MQ=255)
                                                                                         tgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAAc                                                            <  4:201657/149‑1 (MQ=255)
                                                                                           ttCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGa                                                                 <  4:489271/142‑1 (MQ=255)
                                                                                           ttCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGa                                                                 >  3:489271/1‑142 (MQ=255)
                                                                                                   aTTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGt                                                  >  4:291934/1‑149 (MQ=255)
                                                                                                        aCTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTagag                                             <  3:178660/149‑1 (MQ=255)
                                                                                                         cTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGttt                                                                                     <  1:220802/108‑1 (MQ=255)
                                                                                                         cTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGttt                                                                                     >  2:220802/1‑108 (MQ=255)
                                                                                                            cGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCga                                         <  2:602764/149‑1 (MQ=255)
                                                                                                            cGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCga                                         <  1:531300/149‑1 (MQ=255)
                                                                                                            cGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCga                                         <  2:557751/149‑1 (MQ=255)
                                                                                                            cGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCga                                         <  3:454206/149‑1 (MQ=255)
                                                                                                              gATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCgaga                                       <  2:637269/149‑1 (MQ=255)
                                                                                                                 aCGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCg                                    >  4:97992/1‑149 (MQ=255)
                                                                                                                   gTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGAc                                  <  2:198085/149‑1 (MQ=255)
                                                                                                                    ttGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGAcc                                 <  1:559177/149‑1 (MQ=255)
                                                                                                                                ccgccgGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGcc                     >  4:390061/1‑149 (MQ=255)
                                                                                                                                        ataatGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGCTCGACCATATTGCCCGCCTCTTTc               >  2:491799/1‑147 (MQ=255)
                                                                                                                                        ataatGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTc               <  1:491799/147‑1 (MQ=255)
                                                                                                                                            tGGCGAAATACTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTaga                                              >  1:648494/1‑112 (MQ=255)
                                                                                                                                            tGGCGAAATACTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTaga                                              <  2:145406/112‑1 (MQ=255)
                                                                                                                                            tGGCGAAATACTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTaga                                              <  2:648494/112‑1 (MQ=255)
                                                                                                                                            tGGCGAAATACTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTaga                                              >  1:145406/1‑112 (MQ=255)
                                                                                                                                              gCGAAATACTTCGCCACATCGTTGGCAATGCCGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTg       >  3:596592/1‑149 (MQ=255)
                                                                                                                                                gAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTggg     >  1:456412/1‑149 (MQ=255)
                                                                                                                                                gAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTggg     >  2:327078/1‑149 (MQ=255)
                                                                                                                                                 aaGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGt    >  2:572869/1‑149 (MQ=255)
                                                                                                                                                   gTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGc  <  3:291934/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TAACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGC  >  NC_000913/725129‑725424

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