Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,677,363 C→A 10.2% F75L (TTC→TTA rstB → sensory histidine kinase of RstAB two‑component system, low Mg‑responsive via PhoQP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,677,3630CA10.2% 736.2 / 39.4 216F75L (TTC→TTArstBsensory histidine kinase of RstAB two‑component system, low Mg‑responsive via PhoQP
Reads supporting (aligned to +/- strand):  ref base C (98/96);  new base A (13/9);  total (111/105)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 218 total aligned reads displayed.

GTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCTGGCAGTTGGTCCT  >  CP009273/1677216‑1677510
                                                                                                                                                   |                                                                                                                                                   
   cAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGAt                                                                                                                                                  >  2:1735185/1‑148 (MQ=255)
   cAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGAt                                                                                                                                                  >  1:2335175/1‑148 (MQ=255)
   aaaaTTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGAt                                                                                                                                                  <  2:2341457/147‑1 (MQ=255)
    aaaTTTACCGCCGAACGCGCGGGCAAACAGTCGCTGTATTATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATc                                                                                                                                                 <  1:299659/148‑1 (MQ=255)
    aaaTTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATc                                                                                                                                                 <  1:1106632/148‑1 (MQ=255)
        ttACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCg                                                                                                                                             <  2:1239393/148‑1 (MQ=255)
         tACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCgt                                                                                                                                            <  2:353371/148‑1 (MQ=255)
            cgccgAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTc                                                                                                                                         >  1:763519/1‑148 (MQ=255)
             gccgAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTAGATCTGCGTGTCg                                                                                                                                        <  1:1221260/148‑1 (MQ=255)
                     cgcgGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTg                                                                                                                                <  2:460370/148‑1 (MQ=255)
                     cgcgGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTg                                                                                                                                <  2:497665/148‑1 (MQ=255)
                      gcgGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGa                                                                                                                               <  1:2412284/148‑1 (MQ=255)
                         ggCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCCCGGAGTa                                                                                                                            >  2:2157854/1‑148 (MQ=255)
                          gCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTaa                                                                                                                           >  2:2858917/1‑148 (MQ=255)
                            aaaCAGTCGCTGGGTGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAAt                                                                                                                         >  1:2608490/1‑148 (MQ=255)
                             aaCAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATa                                                                                                                        <  1:1893922/148‑1 (MQ=255)
                              aCAGTCGCTGGATGATTTGATGAACAGTTAGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACNACTGGGGTAAAAATCTGAAAGAGATTGATTTAAATATCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATAc                                                                                                                       >  1:1681707/1‑148 (MQ=255)
                               cAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATAcc                                                                                                                      <  2:1037409/148‑1 (MQ=255)
                               cAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATAcc                                                                                                                      >  2:2396938/1‑148 (MQ=255)
                               cAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATAcc                                                                                                                      <  1:3013568/148‑1 (MQ=255)
                                aGTCGCTGGATTATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCa                                                                                                                     <  1:1757449/148‑1 (MQ=255)
                                aGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCa                                                                                                                     <  2:1338052/148‑1 (MQ=255)
                                 gTCGCTGGATGATTTGATTAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCAt                                                                                                                    <  2:1814189/148‑1 (MQ=255)
                                 gTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCAt                                                                                                                    <  2:3042906/148‑1 (MQ=255)
                                 gTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCAt                                                                                                                    <  2:141527/148‑1 (MQ=255)
                                     cTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTAGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtg                                                                                                                >  1:2817165/1‑148 (MQ=255)
                                       ggATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTAGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgat                                                                                                              >  1:655780/1‑148 (MQ=255)
                                         atgatTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgatga                                                                                                            >  2:2899348/1‑148 (MQ=255)
                                         atgatTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgatga                                                                                                            >  2:1248405/1‑148 (MQ=255)
                                          tgatTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTgatgat                                                                                                           >  1:565732/1‑148 (MQ=255)
                                           gatTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGata                                                                                                          >  1:1627864/1‑148 (MQ=255)
                                           gatTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGata                                                                                                          >  2:615621/1‑148 (MQ=255)
                                            atTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGatat                                                                                                         <  1:2042134/148‑1 (MQ=255)
                                            atTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGatat                                                                                                         >  2:1863729/1‑148 (MQ=255)
                                              ttGATGAACAGTTTGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                                                                       >  1:487768/1‑148 (MQ=255)
                                              ttGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                                                                       <  1:575271/148‑1 (MQ=255)
                                                 atgaACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCa                                                                                                    >  1:2978252/1‑148 (MQ=255)
                                                 atgaACAGTTAGCTTTATCTGATGCGCATCGAATTGCGTGAGATCCCCCCACACGACTGGGTTAAAACTCTGATAGAGATTGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCa                                                                                                    <  2:2524413/148‑1 (MQ=255)
                                                  tgaACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCAt                                                                                                   >  1:827938/1‑148 (MQ=255)
                                                   gaACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTAGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATg                                                                                                  <  2:2374553/148‑1 (MQ=255)
                                                    aaCTGTTCGCTGTATCTGATGCGCAGCGATTTGCGTGAGATCCCCCCCCACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGc                                                                                                 <  2:3199945/148‑1 (MQ=255)
                                                     aCAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                                                                <  1:2977037/148‑1 (MQ=255)
                                                     aCAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                                                                >  1:387649/1‑148 (MQ=255)
                                                     aCAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                                                                <  1:2682401/148‑1 (MQ=255)
                                                        gTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCg                                                                                             <  2:469837/148‑1 (MQ=255)
                                                         ttCGCTGTATCTGATGCGCAGCTAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                                                                            <  1:1026210/148‑1 (MQ=255)
                                                         ttCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                                                                            <  2:1592929/148‑1 (MQ=255)
                                                            gCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTg                                                                                         >  2:2981787/1‑148 (MQ=255)
                                                               gTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGt                                                                                      >  2:1624017/1‑148 (MQ=255)
                                                                  tCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggc                                                                                   >  1:18606/1‑148 (MQ=255)
                                                                   cTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcg                                                                                  >  1:1180265/1‑148 (MQ=255)
                                                                    tGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTAGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcgg                                                                                 <  1:1018514/148‑1 (MQ=255)
                                                                        gcgcAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGaa                                                                             <  1:2163928/148‑1 (MQ=255)
                                                                        gcgcAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACAGACTGCGTGGCGGCGaa                                                                             >  1:2062927/1‑148 (MQ=255)
                                                                          gcAGCGAATGGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAAt                                                                           <  2:1998088/148‑1 (MQ=255)
                                                                            aGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTg                                                                         <  1:2894560/148‑1 (MQ=255)
                                                                                aaTTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTAGATCTGCGTGTCGAGCCACTCAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGACATTGTCGc                                                                     >  1:1705351/1‑148 (MQ=255)
                                                                                 aTTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGcc                                                                    <  2:2703559/148‑1 (MQ=255)
                                                                                     cGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTgg                                                                <  2:763519/148‑1 (MQ=255)
                                                                                      gTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGa                                                               <  2:2658195/148‑1 (MQ=255)
                                                                                        gagaTCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACg                                                             <  1:1350218/148‑1 (MQ=255)
                                                                                         agaTCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGa                                                            <  2:90604/148‑1 (MQ=255)
                                                                                         agaTCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGa                                                            <  2:505293/148‑1 (MQ=255)
                                                                                         agaTCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGa                                                            >  2:1567029/1‑148 (MQ=255)
                                                                                         agaTCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGa                                                            >  2:1401571/1‑148 (MQ=255)
                                                                                          gaTCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGAt                                                           <  2:1188837/148‑1 (MQ=255)
                                                                                           aTCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATc                                                          >  2:3062401/1‑148 (MQ=255)
                                                                                            tCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCa                                                         >  1:2809047/1‑148 (MQ=255)
                                                                                            tCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCa                                                         <  1:1991854/148‑1 (MQ=255)
                                                                                            tCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCa                                                         <  1:1991890/148‑1 (MQ=255)
                                                                                              cccccACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGt                                                       >  2:1001308/1‑148 (MQ=255)
                                                                                                cccACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTac                                                     <  1:2565723/148‑1 (MQ=255)
                                                                                                   acacGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACg                                                  <  2:558909/148‑1 (MQ=255)
                                                                                                     acGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGtt                                                <  1:1269194/148‑1 (MQ=255)
                                                                                                          tGGGGTAAAACTCTGAAAGAGATGGATTTAACTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAATTACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGATATTGTCGCCCTGGACGATCAGTACACGTTTTTGc                                           <  2:3098807/148‑1 (MQ=255)
                                                                                                            gggTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAg                                         >  2:2965637/1‑148 (MQ=255)
                                                                                                            gggTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAg                                         >  1:3097250/1‑148 (MQ=255)
                                                                                                             ggTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGc                                        >  1:1528378/1‑148 (MQ=255)
                                                                                                               tAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGt                                      >  1:2117399/1‑148 (MQ=255)
                                                                                                               tAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGt                                      >  1:2177771/1‑148 (MQ=255)
                                                                                                                 aaaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTAt                                    >  2:273169/1‑148 (MQ=255)
                                                                                                                         aaaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCa                            >  2:3091873/1‑148 (MQ=255)
                                                                                                                         aaaGAGATGGATTTAAATCTCTCTTTAGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCa                            >  1:2512120/1‑148 (MQ=255)
                                                                                                                          aaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAg                           >  2:315597/1‑148 (MQ=255)
                                                                                                                           agagATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGc                          >  2:2046284/1‑148 (MQ=255)
                                                                                                                                tGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTa                     >  2:902451/1‑148 (MQ=255)
                                                                                                                                 ggATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTAc                    >  1:2613822/1‑148 (MQ=255)
                                                                                                                                 ggATTTAAATCTCTCTTTAGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTAc                    >  1:1128263/1‑148 (MQ=255)
                                                                                                                                     ttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGc                >  2:826010/1‑148 (MQ=255)
                                                                                                                                     ttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGc                >  2:1850010/1‑148 (MQ=255)
                                                                                                                                      tAAATCTCTCTTTAGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCt               <  2:2817203/148‑1 (MQ=255)
                                                                                                                                       aaaTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCTg              >  2:1462025/1‑148 (MQ=255)
                                                                                                                                         aTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCTGGc            >  2:686995/1‑148 (MQ=255)
                                                                                                                                          tctctcTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCTGGCa           >  2:1589609/1‑148 (MQ=255)
                                                                                                                                             ctctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCTGGCAggt        >  2:1689787/1‑146 (MQ=255)
                                                                                                                                             ctctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCTGGCAGtt        >  1:358185/1‑148 (MQ=255)
                                                                                                                                              tctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCTGGCAGTTg       <  2:2795038/148‑1 (MQ=255)
                                                                                                                                                 ttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCTGGCAGTTGGTc    <  1:1850010/148‑1 (MQ=255)
                                                                                                                                                 ttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCTGGCAGTTGGTc    <  1:1579172/148‑1 (MQ=255)
                                                                                                                                                  tCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCTGGCAGTTGGTcc   <  1:316626/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTATCCCGCGCAGCCACTACGTGCTGGCAGTTGGTCCT  >  CP009273/1677216‑1677510

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