Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 871,072 C→A 7.6% V94V (GTC→GTA yliF → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273871,0720CA 7.6% 1134.0 / 50.9 315V94V (GTC→GTAyliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (121/170);  new base A (6/18);  total (127/188)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 322 total aligned reads displayed.

ACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCT  >  CP009273/870923‑871221
                                                                                                                                                     |                                                                                                                                                     
aCTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTc                                                                                                                                                       <  1:3297095/150‑1 (MQ=255)
 cTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCg                                                                                                                                                      <  2:691411/150‑1 (MQ=255)
  tGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGAGGTCgg                                                                                                                                                     >  1:117137/1‑150 (MQ=255)
       ttGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCtt                                                                                                                                                <  1:3876632/150‑1 (MQ=255)
           ggTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCt                                                                                                                                            >  2:149253/1‑150 (MQ=255)
              ggCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGAt                                                                                                                                         >  2:2480730/1‑150 (MQ=255)
               gCAAATATGGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGAtt                                                                                                                                        <  2:3175221/150‑1 (MQ=255)
                 aaaTATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTc                                                                                                                                      >  1:3527329/1‑150 (MQ=255)
                 aaaTATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTc                                                                                                                                      >  1:2938243/1‑150 (MQ=255)
                 aaaTATTGCAGAATATCGGCGAAGCGATGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGATGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTAATTTGTACGATGTCGGTCCTTATCTGATTTc                                                                                                                                      >  2:1402663/1‑150 (MQ=255)
                 aaaTATTGCAGAAGATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTc                                                                                                                                      <  2:3947359/150‑1 (MQ=255)
                  aaTATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTcc                                                                                                                                     <  1:4002623/150‑1 (MQ=255)
                     atTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGt                                                                                                                                  <  2:5088451/150‑1 (MQ=255)
                     atTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGt                                                                                                                                  <  1:4866597/150‑1 (MQ=255)
                       tgatGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGa                                                                                                                                >  1:636686/4‑150 (MQ=255)
                       tGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGa                                                                                                                                >  1:3180329/1‑150 (MQ=255)
                           gAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGac                                                                                                                            >  1:4409929/1‑149 (MQ=255)
                           gAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGaa                                                                                                                            >  1:956127/1‑150 (MQ=255)
                            aaTATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAAt                                                                                                                           <  2:4621732/150‑1 (MQ=255)
                            aaTATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAAt                                                                                                                           <  2:2268754/150‑1 (MQ=255)
                            aaTATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAAt                                                                                                                           >  2:369772/1‑150 (MQ=255)
                            aaTATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAAt                                                                                                                           <  1:1125569/150‑1 (MQ=255)
                                tCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTAt                                                                                                                       <  2:3308573/150‑1 (MQ=255)
                                    cGTAGCGTTGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAg                                                                                                                   <  1:4225723/150‑1 (MQ=255)
                                       aGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTa                                                                                                                >  1:2430923/1‑150 (MQ=255)
                                          gTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAg                                                                                                             <  2:4348531/150‑1 (MQ=255)
                                          gTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGATTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAg                                                                                                             >  2:3507570/1‑150 (MQ=255)
                                           tGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAgg                                                                                                            >  1:5177772/1‑150 (MQ=255)
                                            ggCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGa                                                                                                           <  1:5040401/150‑1 (MQ=255)
                                                aTTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTc                                                                                                       <  2:938260/150‑1 (MQ=255)
                                                aTTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTc                                                                                                       >  1:1670699/1‑150 (MQ=255)
                                                aTTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTc                                                                                                       <  1:3202547/150‑1 (MQ=255)
                                                aTTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTc                                                                                                       <  1:2803738/150‑1 (MQ=255)
                                                  tGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCTGGTAAAGGAGTTCga                                                                                                     >  1:3723442/1‑150 (MQ=255)
                                                   gATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCgag                                                                                                    <  2:2247179/150‑1 (MQ=255)
                                                       ttATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGaaaa                                                                                                <  1:1620695/150‑1 (MQ=255)
                                                         aTCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAat                                                                                              >  2:3625929/1‑150 (MQ=255)
                                                         aTCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAat                                                                                              >  2:1343426/1‑150 (MQ=255)
                                                           cAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAatta                                                                                            <  2:3152063/150‑1 (MQ=255)
                                                            aaCACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAattat                                                                                           <  1:4596074/150‑1 (MQ=255)
                                                            aaCACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAagtat                                                                                           >  1:3077058/1‑150 (MQ=255)
                                                            aaCACGCTAAATCCCACGGAGGAACAGCGGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAattat                                                                                           >  1:144456/1‑150 (MQ=255)
                                                             acacGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAattatt                                                                                          >  2:2303578/1‑150 (MQ=255)
                                                             acacGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAattatt                                                                                          <  2:3839661/150‑1 (MQ=255)
                                                             acacGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAattatt                                                                                          <  2:290933/150‑1 (MQ=255)
                                                             acacGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAattatt                                                                                          <  2:3083162/150‑1 (MQ=255)
                                                             acacGCTAAATCCCACGGAGGAACAGCGGTGGGTTGGTTTGCGCACAGCTTACGCCGACTCGGGTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGGTCGAGAAAAattatt                                                                                          >  2:1733580/1‑150 (MQ=255)
                                                              cacGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATtg                                                                                         >  2:3368932/1‑150 (MQ=255)
                                                              cacGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATtg                                                                                         >  2:1311488/1‑150 (MQ=255)
                                                               acGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATtgt                                                                                        <  1:4661560/150‑1 (MQ=255)
                                                                   naaaTCCCACGGAGGAACAGCTGTTGGTTGTTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCCTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGctg                                                                                    >  2:2399512/2‑148 (MQ=255)
                                                                    aaaTCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGa                                                                                   >  2:2524555/1‑150 (MQ=255)
                                                                    aaaTCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGa                                                                                   >  2:2041676/1‑150 (MQ=255)
                                                                     aaTCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGat                                                                                  <  2:1238560/150‑1 (MQ=255)
                                                                     aaTCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGat                                                                                  <  2:1922561/150‑1 (MQ=255)
                                                                      ctCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGata                                                                                 <  1:1385933/149‑1 (MQ=255)
                                                                        cccACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGAtatt                                                                               <  1:4749583/150‑1 (MQ=255)
                                                                         ccACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGAtatta                                                                              <  2:88012/150‑1 (MQ=255)
                                                                          cACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCTGTCCTTATCTGATTTCCAGTGACGAATGTATTAAGGTAAAGTAGTTAGAGAAAAATTATTGTGAAGAtattat                                                                             >  2:2715451/1‑150 (MQ=255)
                                                                           aCGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATg                                                                            <  1:1852538/150‑1 (MQ=255)
                                                                                 gaACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGtt                                                                      <  1:2854655/150‑1 (MQ=255)
                                                                                  aaCAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTtg                                                                     <  2:4934011/150‑1 (MQ=255)
                                                                                    cAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTtgtg                                                                   <  1:2509186/150‑1 (MQ=255)
                                                                                    cAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTtgtg                                                                   <  1:2910923/150‑1 (MQ=255)
                                                                                      gCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGAAGATATTATGCAGGTTGTGaa                                                                 >  2:3935536/1‑150 (MQ=255)
                                                                                           tggttggtTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATc                                                            <  1:2058334/150‑1 (MQ=255)
                                                                                            ggttggttTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCg                                                           >  1:4532203/1‑150 (MQ=255)
                                                                                               tggttTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACa                                                        >  2:612671/1‑150 (MQ=255)
                                                                                                 gttTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGACAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATg                                                      >  2:4029044/1‑150 (MQ=255)
                                                                                                 gttTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATg                                                      >  1:836775/1‑150 (MQ=255)
                                                                                                 gttTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATg                                                      <  1:1343284/150‑1 (MQ=255)
                                                                                                 gttTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATg                                                      >  2:1158560/1‑150 (MQ=255)
                                                                                                  tttGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGt                                                     <  2:4147649/150‑1 (MQ=255)
                                                                                                  tttGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGt                                                     <  2:4545897/150‑1 (MQ=255)
                                                                                                   ttGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCTGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATAGACATGTc                                                    <  2:343825/150‑1 (MQ=255)
                                                                                                       gcACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCaaaa                                                >  1:260166/1‑150 (MQ=255)
                                                                                                        cacaGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCaaaaa                                               >  1:1965829/1‑150 (MQ=255)
                                                                                                         acaGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATCTTATGCAGGTTGTGAAGTATCGACATGTCAAAAAt                                              >  1:1522537/1‑150 (MQ=255)
                                                                                                               tACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAggg                                        <  2:3353712/150‑1 (MQ=255)
                                                                                                               tACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAggg                                        <  2:3713176/150‑1 (MQ=255)
                                                                                                                 cGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGtt                                      <  2:2155805/150‑1 (MQ=255)
                                                                                                                 cGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGtt                                      <  2:5710/150‑1 (MQ=255)
                                                                                                                  gCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGttt                                     >  2:4971733/1‑150 (MQ=255)
                                                                                                                    cGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTAt                                   >  2:4662198/1‑150 (MQ=255)
                                                                                                                       cTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTAtctc                                <  2:4773272/150‑1 (MQ=255)
                                                                                                                               cccccTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACg                        >  1:2015841/1‑150 (MQ=255)
                                                                                                                               cccccTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACg                        <  1:4688873/150‑1 (MQ=255)
                                                                                                                                    tctTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTaaa                   <  2:2642124/150‑1 (MQ=255)
                                                                                                                                     ctTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTaaaa                  >  1:1321617/1‑150 (MQ=255)
                                                                                                                                     ctTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTaaaa                  >  1:3332557/1‑150 (MQ=255)
                                                                                                                                      ttATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAAc                 <  2:4443978/150‑1 (MQ=255)
                                                                                                                                            gTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGt           >  1:2212940/1‑150 (MQ=255)
                                                                                                                                             tACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTc          >  2:1332202/1‑150 (MQ=255)
                                                                                                                                              aCGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCt         <  1:271016/150‑1 (MQ=255)
                                                                                                                                              aCGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCt         <  1:3395633/150‑1 (MQ=255)
                                                                                                                                                gATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTAt       >  2:1954753/1‑150 (MQ=255)
                                                                                                                                                  tGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTa     <  2:5194859/150‑1 (MQ=255)
                                                                                                                                                   gTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTAc    <  2:1075267/150‑1 (MQ=255)
                                                                                                                                                   gTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTAc    >  1:4045148/1‑150 (MQ=255)
                                                                                                                                                     cGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACct  <  1:346075/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCT  >  CP009273/870923‑871221

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: