Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 725,276 G→A 7.8% Y587Y (TAC→TAT kdpB ← K(+) transporting P‑type ATPase subunit KdpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913725,2760GA 7.8% 273.3 / 7.8 89Y587Y (TAC→TATkdpBK(+) transporting P‑type ATPase subunit KdpB
Reads supporting (aligned to +/- strand):  ref base G (40/42);  new base A (5/2);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TTTAACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGC  >  NC_000913/725127‑725424
                                                                                                                                                     |                                                                                                                                                    
tttAACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAg                                                                                                                                                      <  1:1476580/150‑1 (MQ=255)
tttAACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAg                                                                                                                                                      <  2:866480/150‑1 (MQ=255)
     cGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTAttt                                                                                                                                                 <  1:69819/150‑1 (MQ=255)
     cGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCACATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTAttt                                                                                                                                                 >  1:1409831/1‑150 (MQ=255)
       ccAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCg                                                                                                                                               >  1:81370/1‑150 (MQ=255)
           gggAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAATACTTCGCcat                                                                                                                                           <  1:289157/150‑2 (MQ=255)
            ggAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCGCTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAATACTTCGCcatt                                                                                                                                          <  1:333462/150‑3 (MQ=255)
              aaTCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATc                                                                                                                                        >  1:747156/1‑150 (MQ=255)
               aTCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCg                                                                                                                                       >  1:826974/1‑150 (MQ=255)
                 cAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGtt                                                                                                                                     >  2:1470206/1‑150 (MQ=255)
                 cAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGtt                                                                                                                                     <  1:905822/150‑1 (MQ=255)
                 cAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGtt                                                                                                                                     >  2:905822/1‑150 (MQ=255)
                  aaaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGt                                                                                                                                      >  2:961046/1‑148 (MQ=255)
                  aaaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGt                                                                                                                                      <  1:961046/148‑1 (MQ=255)
                  aaaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTg                                                                                                                                    <  2:289633/150‑1 (MQ=255)
                   aaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTgg                                                                                                                                   <  2:344883/150‑1 (MQ=255)
                    aaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGc                                                                                                                                  <  2:126446/150‑1 (MQ=255)
                    aaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGc                                                                                                                                  <  2:479093/150‑1 (MQ=255)
                    aaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGc                                                                                                                                  <  2:1122717/150‑1 (MQ=255)
                         cGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGc                                                                                                                             >  2:1071552/1‑150 (MQ=255)
                         cGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGc                                                                                                                             <  2:1610634/150‑1 (MQ=255)
                            tAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGa                                                                                                                          <  1:87980/150‑1 (MQ=255)
                             aaTCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGaa                                                                                                                         >  2:1262496/1‑150 (MQ=255)
                              aTCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAg                                                                                                                        <  2:1059029/150‑1 (MQ=255)
                               tCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAgg                                                                                                                       <  1:1233097/150‑1 (MQ=255)
                                 aaGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAggtg                                                                                                                     >  1:92699/1‑150 (MQ=255)
                                     cGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCa                                                                                                                 >  1:1294045/1‑150 (MQ=255)
                                        tGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCCGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCg                                                                                                              >  2:1475569/1‑150 (MQ=255)
                                         gAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGa                                                                                                             >  2:1079626/1‑150 (MQ=255)
                                            aaTCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGcc                                                                                                          <  2:67795/150‑1 (MQ=255)
                                               cACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACg                                                                                                       >  2:134251/1‑150 (MQ=255)
                                               cACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACg                                                                                                       <  2:555534/150‑1 (MQ=255)
                                                 ccGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGGTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGGGGGCAGCGAGCCACggg                                                                                                     <  1:57847/150‑1 (MQ=255)
                                                   gCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGCATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGTAATGCTGAAGGTGGTCAGCGAGCCACGGGtc                                                                                                   >  2:1490650/1‑150 (MQ=255)
                                                   gCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGtc                                                                                                   <  2:1048093/150‑1 (MQ=255)
                                                               gCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGt                                                                                       >  1:338005/1‑150 (MQ=255)
                                                                cGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGtt                                                                                      >  1:856266/1‑150 (MQ=255)
                                                                cGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGtt                                                                                      >  1:1448142/1‑150 (MQ=255)
                                                                  gAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTg                                                                                    <  1:134251/150‑1 (MQ=255)
                                                                    gTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGcc                                                                                  >  2:342305/1‑150 (MQ=255)
                                                                      cGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCaa                                                                                <  2:175362/150‑1 (MQ=255)
                                                                        ggCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAAtg                                                                              >  1:1157422/1‑150 (MQ=255)
                                                                        ggCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAAtg                                                                              <  2:1294045/150‑1 (MQ=255)
                                                                        ggCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAAtg                                                                              >  1:398557/1‑150 (MQ=255)
                                                                           gAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGc                                                                           >  1:1251377/1‑150 (MQ=255)
                                                                            aaTGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGca                                                                          <  1:1209597/150‑1 (MQ=255)
                                                                              tGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGAGcacc                                                                        >  1:1066428/1‑150 (MQ=255)
                                                                               gCAGGCACATGATGTTCAGCGCATTTAACTGAGGATACGTTGCCGCGAATGCCGCCGGAATAATGTCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGcacca                                                                       >  2:255303/1‑150 (MQ=255)
                                                                                     acaTGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGa                                                                 <  2:44696/150‑1 (MQ=255)
                                                                                      caTGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGAt                                                                <  1:55673/150‑1 (MQ=255)
                                                                                       atgatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATc                                                               >  1:1339310/1‑150 (MQ=255)
                                                                                       atgatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATc                                                               <  1:819488/150‑1 (MQ=255)
                                                                                            gTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACtt                                                          >  2:531616/1‑150 (MQ=255)
                                                                                               cAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGt                                                       >  2:990381/1‑150 (MQ=255)
                                                                                                 gcgcATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCg                                                     <  2:856266/150‑1 (MQ=255)
                                                                                                  cgcATTTAACTGCGGATACGTTGCCGCGAATGCTGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCgg                                                    >  1:462753/1‑150 (MQ=255)
                                                                                                   gcATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCggg                                                   >  1:1328872/1‑150 (MQ=255)
                                                                                                    cATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGt                                                  <  1:873708/150‑1 (MQ=255)
                                                                                                    cATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGt                                                  >  1:80584/1‑150 (MQ=255)
                                                                                                       ttAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTag                                               <  1:1071552/150‑1 (MQ=255)
                                                                                                        tAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTaga                                              <  2:1601233/150‑1 (MQ=255)
                                                                                                          aCTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGt                                            <  1:1262496/150‑1 (MQ=255)
                                                                                                           cTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGcc                                                                                  <  1:1144557/111‑1 (MQ=255)
                                                                                                           cTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGcc                                                                                  >  2:1144557/1‑111 (MQ=255)
                                                                                                            tGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCg                                          <  2:92699/150‑1 (MQ=255)
                                                                                                             gCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCga                                         <  2:1409831/150‑1 (MQ=255)
                                                                                                              cGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCgag                                        <  2:747156/150‑1 (MQ=255)
                                                                                                               ggATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCgaga                                       >  2:122115/1‑150 (MQ=255)
                                                                                                                gATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGAt                                      <  2:503150/150‑1 (MQ=255)
                                                                                                                     gTTGCCGCGAATGCCGCCGGAAGAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTTGTCAGCGATCCACGTGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGAcc                                 <  2:1121113/150‑1 (MQ=255)
                                                                                                                      ttGCCGCGAATGCTGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCa                                <  2:462753/150‑1 (MQ=255)
                                                                                                                      ttGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCa                                >  1:801410/1‑150 (MQ=255)
                                                                                                                        gcagcgAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACAGCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGGTTGCCAATGTGCACCACCTCGATAAACTTGGTCGGGTTAGAGTCGAGCGCGACCata                              >  2:844575/4‑150 (MQ=255)
                                                                                                                         ccGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCatat                             >  2:438821/1‑150 (MQ=255)
                                                                                                                            cgAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGc                          <  1:990381/150‑1 (MQ=255)
                                                                                                                              aaTGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGccc                        <  2:451971/150‑1 (MQ=255)
                                                                                                                              aaTGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCGGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGccc                        <  2:920754/150‑1 (MQ=255)
                                                                                                                                   cgccgGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCctc                   >  1:1444242/1‑150 (MQ=255)
                                                                                                                                     ccgGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATACTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCCCCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGCCCATATTGCCCGCCTCtt                 >  2:1333074/1‑150 (MQ=255)
                                                                                                                                        acgatgtGGCGAAATACTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTcg              >  1:901870/7‑150 (MQ=255)
                                                                                                                                          ataatGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCCCCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTcgcc            >  1:1327043/1‑150 (MQ=255)
                                                                                                                                           atgtGGCGAAATACTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTcgccg           >  1:527037/4‑150 (MQ=255)
                                                                                                                                           atgtGGCGAAATACTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCAACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTcgccg           >  2:1032343/4‑150 (MQ=255)
                                                                                                                                            tgtGGCGAAATACTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTcgccgc          >  1:1370873/3‑150 (MQ=255)
                                                                                                                                             atGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGAc                                  <  2:1413376/125‑1 (MQ=255)
                                                                                                                                             atGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGAc                                  >  1:1413376/1‑125 (MQ=255)
                                                                                                                                             atGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCgccgcc         <  2:801410/150‑1 (MQ=255)
                                                                                                                                               ggCGAAATACTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTg       >  1:607790/1‑150 (MQ=255)
                                                                                                                                                 cGAGGTATTTCGCCACATCGTTGGCAATGCTGAAGGGGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTggg     <  1:1542593/150‑1 (MQ=255)
                                                                                                                                                 cGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTggg     <  2:1018405/150‑1 (MQ=255)
                                                                                                                                                    aGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGc  >  1:1611475/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTAACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGC  >  NC_000913/725127‑725424

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