Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 719,056 C→T 6.8% R272Q (CGG→CAG)  kdpD ← fused sensory histidine kinase in two‑component regulatory system with KdpE: signal sensing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273719,0560CT 6.8% 899.5 / 33.2 251R272Q (CGG→CAG) kdpDfused sensory histidine kinase in two‑component regulatory system with KdpE: signal sensing protein
Reads supporting (aligned to +/- strand):  ref base C (115/119);  new base T (4/13);  total (119/132)
Fisher's exact test for biased strand distribution p-value = 4.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

TCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGGCGC  >  CP009273/718909‑719202
                                                                                                                                                   |                                                                                                                                                  
tCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAgccg                                                                                                                                                   >  2:329904/1‑149 (MQ=255)
   gAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAgccgcgc                                                                                                                                                >  1:318252/1‑149 (MQ=255)
       gtgtTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCt                                                                                                                                            >  2:1951606/1‑149 (MQ=255)
         gtTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTgc                                                                                                                                          >  2:3614177/1‑149 (MQ=255)
           tGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTgcgc                                                                                                                                        >  2:3652817/1‑149 (MQ=255)
            gCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGcgcg                                                                                                                                       >  2:1171754/1‑149 (MQ=255)
            gCCGTCTCCGCGCCCAGGTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGcgcg                                                                                                                                       <  2:3230421/149‑1 (MQ=255)
                tctcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGAcc                                                                                                                                   >  2:2864962/1‑149 (MQ=255)
                tctcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGAcc                                                                                                                                   >  1:2221527/1‑149 (MQ=255)
                    cgcgCCCAGGTCCTGCGCCAGACGTACGGCGCTGAGAATTGCCCGACGTTTGTTGTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGtt                                                                                                                               <  1:4642254/149‑1 (MQ=255)
                       gcCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGttttt                                                                                                                            >  1:3656822/1‑149 (MQ=255)
                       gcCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGttttt                                                                                                                            <  2:3787706/149‑1 (MQ=255)
                       gcCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCGACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGttttt                                                                                                                            >  1:1610727/1‑149 (MQ=255)
                          cAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGc                                                                                                                         >  2:2563530/1‑149 (MQ=255)
                          cAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGc                                                                                                                         >  2:171766/1‑149 (MQ=255)
                                cTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCgg                                                                                                                   <  2:1759850/149‑1 (MQ=255)
                                 tGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCCGCTGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGgt                                                                                                                  >  2:18177/1‑149 (MQ=255)
                                  gcgcCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGgtg                                                                                                                 >  1:2714037/1‑149 (MQ=255)
                                  gcgcCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGgtg                                                                                                                 >  2:2930867/1‑149 (MQ=255)
                                  gcgcCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGgtg                                                                                                                 <  2:1077635/149‑1 (MQ=255)
                                  gcgcCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGgtg                                                                                                                 <  1:218985/149‑1 (MQ=255)
                                   cgcCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGgtgt                                                                                                                <  2:4660865/149‑1 (MQ=255)
                                        gACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATgg                                                                                                           >  2:243550/1‑149 (MQ=255)
                                        gACGTAAGGAGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGTTGCAGGGCAGGGGTTTCAACATACACCGCGTTCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATgg                                                                                                           >  2:4294593/1‑149 (MQ=255)
                                         aCGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGc                                                                                                          >  2:1910237/1‑149 (MQ=255)
                                         aCGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGc                                                                                                          <  2:1829405/149‑1 (MQ=255)
                                         aCGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGc                                                                                                          >  2:606587/1‑149 (MQ=255)
                                            tAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGa                                                                                                       <  1:2682416/149‑1 (MQ=255)
                                            tAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGa                                                                                                       <  1:2980577/149‑1 (MQ=255)
                                            tAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGa                                                                                                       >  2:113278/1‑149 (MQ=255)
                                            tAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGa                                                                                                       <  2:3709309/149‑1 (MQ=255)
                                              aGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATg                                                                                                     >  2:476196/1‑149 (MQ=255)
                                               ggCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGc                                                                                                    <  1:3742151/149‑1 (MQ=255)
                                                  gcTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATa                                                                                                 >  2:1906394/1‑149 (MQ=255)
                                                     gagaATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATaaaa                                                                                              <  2:1469413/149‑1 (MQ=255)
                                                     gagaATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATaaaa                                                                                              >  1:4027726/1‑149 (MQ=255)
                                                     gagaATTGCCCGACGTTTTTGTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCGACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATaaaa                                                                                              <  1:4325245/149‑1 (MQ=255)
                                                       gaATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCTGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAgg                                                                                            <  2:3228234/149‑1 (MQ=255)
                                                        aaTTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGa                                                                                           <  1:3562460/149‑1 (MQ=255)
                                                         aTTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAt                                                                                          <  1:29274/149‑1 (MQ=255)
                                                          ttGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAtc                                                                                         <  2:1016156/149‑1 (MQ=255)
                                                            gCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAtcgc                                                                                       >  2:2091921/1‑149 (MQ=255)
                                                              ccGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAtcgcgt                                                                                     <  1:1916296/149‑1 (MQ=255)
                                                              ccGACGTTATTTGTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAtcgcgt                                                                                     <  1:1793142/149‑1 (MQ=255)
                                                                  cGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAtcgcgtcgcg                                                                                 >  2:4645833/1‑149 (MQ=255)
                                                                   gTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTcgcgc                                                                                <  1:4532801/149‑1 (MQ=255)
                                                                    ttttttttCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTcgcgcg                                                                               >  1:674103/1‑149 (MQ=255)
                                                                     tttttttCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCgt                                                                              <  1:2340426/149‑1 (MQ=255)
                                                                        ttttCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGtgtg                                                                           >  1:3629214/1‑149 (MQ=255)
                                                                          ttCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGcc                                                                         <  1:3332249/149‑1 (MQ=255)
                                                                           tCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCca                                                                        <  2:261576/149‑1 (MQ=255)
                                                                            ccGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCTGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCcac                                                                       <  2:74274/149‑1 (MQ=255)
                                                                              ggTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCacac                                                                     >  2:3693030/1‑149 (MQ=255)
                                                                               gTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCTCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCAACACt                                                                    >  2:1384360/1‑149 (MQ=255)
                                                                                tAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACttt                                                                  <  1:4525996/149‑1 (MQ=255)
                                                                                   gCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTtct                                                                >  1:1133421/1‑149 (MQ=255)
                                                                                    cGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTtctc                                                               <  1:2086661/149‑1 (MQ=255)
                                                                                     ggTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCTGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTctct                                                              <  1:1049953/149‑1 (MQ=255)
                                                                                     ggTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCTGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTctct                                                              <  1:3438790/149‑1 (MQ=255)
                                                                                     ggTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTctct                                                              <  2:1661793/149‑1 (MQ=255)
                                                                                      gTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCtt                                                             <  2:3468115/149‑1 (MQ=255)
                                                                                      gTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCtt                                                             <  2:3048267/149‑1 (MQ=255)
                                                                                       tgtagggcagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTc                                                            <  1:2588959/146‑1 (MQ=255)
                                                                                         cagggcagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGc                                                          <  1:3460654/149‑1 (MQ=255)
                                                                                           gggcagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCa                                                        >  1:2906292/1‑149 (MQ=255)
                                                                                           gggcagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCa                                                        >  1:3366937/1‑149 (MQ=255)
                                                                                           gggcagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCa                                                        >  2:3238487/1‑149 (MQ=255)
                                                                                              cagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGa                                                     <  2:3421649/149‑1 (MQ=255)
                                                                                                    tttCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCTGCGCCGCTGCGCGGACCAGTTTTGCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGcccc                                               <  2:4204845/149‑1 (MQ=255)
                                                                                                    tttCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGcccc                                               >  1:4584824/1‑149 (MQ=255)
                                                                                                     ttCAACAGACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGccccc                                              <  2:1575362/149‑1 (MQ=255)
                                                                                                        aaCATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGcc                                           >  1:950579/1‑149 (MQ=255)
                                                                                                         aCATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCa                                          <  1:1021756/149‑1 (MQ=255)
                                                                                                          cATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAg                                         <  1:4645788/149‑1 (MQ=255)
                                                                                                          cATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAg                                         >  1:624547/1‑149 (MQ=255)
                                                                                                            tACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGgc                                       <  1:1242293/149‑1 (MQ=255)
                                                                                                              cacCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGgcgc                                     >  2:1856597/1‑149 (MQ=255)
                                                                                                               acCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGgcgcg                                    >  2:2363184/1‑149 (MQ=255)
                                                                                                                  gcgTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCAt                                 <  1:2283253/149‑1 (MQ=255)
                                                                                                                  gcgTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCAt                                 <  1:4680561/149‑1 (MQ=255)
                                                                                                                   cgTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCAtt                                <  1:4379742/149‑1 (MQ=255)
                                                                                                                   cgTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCAtt                                <  1:1687170/149‑1 (MQ=255)
                                                                                                                         aGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCtca                          <  1:295422/149‑1 (MQ=255)
                                                                                                                          gACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCtcat                         >  2:2934087/1‑149 (MQ=255)
                                                                                                                           aCGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCtcatc                        <  1:2793383/149‑1 (MQ=255)
                                                                                                                             gCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCaa                      >  1:4263549/1‑149 (MQ=255)
                                                                                                                                aCCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAAcgc                   >  2:3402766/1‑149 (MQ=255)
                                                                                                                                 cccAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAAcgcg                  <  1:2122130/149‑1 (MQ=255)
                                                                                                                                  ccAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGa                 <  1:2930875/149‑1 (MQ=255)
                                                                                                                                    aGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATc               <  2:683468/149‑1 (MQ=255)
                                                                                                                                    aGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATc               >  2:826020/1‑149 (MQ=255)
                                                                                                                                      ccGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCgg             >  1:498181/1‑149 (MQ=255)
                                                                                                                                        gTGACGCCAGCTGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCa           <  2:3867925/149‑1 (MQ=255)
                                                                                                                                        gTGACGCCAGCTGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCa           <  2:4596164/149‑1 (MQ=255)
                                                                                                                                        gTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCa           >  1:401280/1‑149 (MQ=255)
                                                                                                                                          gACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGt         <  1:2564089/149‑1 (MQ=255)
                                                                                                                                           aCGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACATTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTa        <  2:4026151/149‑1 (MQ=255)
                                                                                                                                              ccAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACgg     <  2:2233535/149‑1 (MQ=255)
                                                                                                                                                 gccgcgccgcTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGgcgc  >  2:4357355/1‑149 (MQ=255)
                                                                                                                                                 gccgcgccgcTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGgcgc  >  2:2443242/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCAGAAAGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGGCGC  >  CP009273/718909‑719202

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: