Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 871,072 C→A 10.5% 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,0720CA10.5% 1271.1 / 87.7 362V94V (GTC→GTAyliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (137/187);  new base A (15/23);  total (152/210)
Fisher's exact test for biased strand distribution p-value = 8.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 370 total aligned reads displayed.

ACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCT  >  CP009273/870923‑871221
                                                                                                                                                     |                                                                                                                                                     
aCTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTc                                                                                                                                                       <  1:4844230/150‑1 (MQ=255)
aCTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTc                                                                                                                                                       >  2:2027812/1‑150 (MQ=255)
aCTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTc                                                                                                                                                       <  2:2261674/150‑1 (MQ=255)
 cTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCAAAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCg                                                                                                                                                      <  2:4386185/150‑1 (MQ=255)
    cTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTc                                                                                                                                                   >  2:1079568/1‑150 (MQ=255)
     ttttGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTcc                                                                                                                                                  <  2:3395663/150‑1 (MQ=255)
      tttGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCt                                                                                                                                                 <  2:3325622/150‑1 (MQ=255)
         gCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTAt                                                                                                                                              <  1:4562287/150‑1 (MQ=255)
           ggTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCt                                                                                                                                            <  1:3175449/150‑1 (MQ=255)
            gTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTg                                                                                                                                           <  2:3234247/150‑1 (MQ=255)
             tGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGa                                                                                                                                          >  2:89926/1‑150 (MQ=255)
               gCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGAtt                                                                                                                                        >  1:1353382/1‑150 (MQ=255)
                   atatTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGACCAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCa                                                                                                                                    >  1:827251/1‑150 (MQ=255)
                   atatTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCa                                                                                                                                    <  1:1849572/150‑1 (MQ=255)
                   atatTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCa                                                                                                                                    <  2:3551612/150‑1 (MQ=255)
                    tatTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAg                                                                                                                                   >  1:4300716/1‑150 (MQ=255)
                    tatTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAg                                                                                                                                   <  2:924800/150‑1 (MQ=255)
                    tatTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAg                                                                                                                                   <  1:807575/150‑1 (MQ=255)
                      ttGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTg                                                                                                                                 >  1:208444/1‑150 (MQ=255)
                      ttGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTg                                                                                                                                 >  1:3018687/1‑150 (MQ=255)
                       tGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGa                                                                                                                                >  2:2644267/1‑150 (MQ=255)
                        gCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGAc                                                                                                                               >  1:584066/1‑150 (MQ=255)
                          tGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGa                                                                                                                             >  1:2406874/1‑150 (MQ=255)
                             atatCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCCACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATg                                                                                                                          >  1:4981662/1‑150 (MQ=255)
                              tatCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGt                                                                                                                         <  2:3264918/150‑1 (MQ=255)
                               atCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTa                                                                                                                        <  2:1557043/150‑1 (MQ=255)
                               atCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTa                                                                                                                        >  1:3808002/1‑150 (MQ=255)
                                tCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTAt                                                                                                                       <  2:2908275/150‑1 (MQ=255)
                                 cGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTAtt                                                                                                                      <  1:1166953/150‑1 (MQ=255)
                                  ggCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTc                                                                                                                     >  1:2128847/1‑150 (MQ=255)
                                          gTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAg                                                                                                             >  1:1931265/1‑150 (MQ=255)
                                           tGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAgg                                                                                                            >  1:906405/1‑150 (MQ=255)
                                                aTTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTc                                                                                                       <  2:5146280/150‑1 (MQ=255)
                                                    aTCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCgaga                                                                                                   <  1:1238272/150‑1 (MQ=255)
                                                    aTCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCgaga                                                                                                   >  1:1802041/1‑150 (MQ=255)
                                                     tCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGaa                                                                                                  <  2:1768297/150‑1 (MQ=255)
                                                       ttATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTACAGGAGTTCGAGaaaa                                                                                                >  2:76647/1‑150 (MQ=255)
                                                       ttATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGaaaa                                                                                                >  2:4290083/1‑150 (MQ=255)
                                                        tATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGaaaaa                                                                                               <  2:694963/150‑1 (MQ=255)
                                                          tCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAatt                                                                                             >  2:5490879/1‑150 (MQ=255)
                                                          tCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAatt                                                                                             <  1:3979223/150‑1 (MQ=255)
                                                           cAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAatta                                                                                            <  2:1325128/150‑1 (MQ=255)
                                                             acacGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATCCAGGTAAAGGAGTTCCAGAAAAattatt                                                                                          >  1:1690108/1‑150 (MQ=255)
                                                              cacGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAATTTATtg                                                                                         >  2:4497431/1‑150 (MQ=255)
                                                              cacGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATtg                                                                                         >  2:1664801/1‑150 (MQ=255)
                                                               acGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATtgt                                                                                        <  1:2047778/150‑1 (MQ=255)
                                                                 gCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTTTGc                                                                                      >  2:3831159/1‑150 (MQ=255)
                                                                 gCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGc                                                                                      <  2:5089884/150‑1 (MQ=255)
                                                                     aaTCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGat                                                                                  <  1:1576259/150‑1 (MQ=255)
                                                                          cACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGAtattat                                                                             <  2:147481/150‑1 (MQ=255)
                                                                          cACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGAtattat                                                                             <  1:2773794/150‑1 (MQ=255)
                                                                            cGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGc                                                                           >  2:4839446/1‑150 (MQ=255)
                                                                               aggaACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAgg                                                                        >  1:3422094/1‑150 (MQ=255)
                                                                               aggaACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAgg                                                                        >  1:3087677/1‑150 (MQ=255)
                                                                                ggaACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGt                                                                       <  1:4788301/150‑1 (MQ=255)
                                                                                  aaCAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTtg                                                                     <  2:400201/150‑1 (MQ=255)
                                                                                  aaCAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGACGGTCCTTATCTGAATTCCAGTGACGAATGTATTAAGGTAAAGGAGTTCGAGAAAAATTAATGTGCAGATATTATGCAGGTtg                                                                     >  1:2815369/1‑150 (MQ=255)
                                                                                    cAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTtgtg                                                                   >  1:3598729/1‑150 (MQ=255)
                                                                                     aGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGa                                                                  >  2:2476069/1‑150 (MQ=255)
                                                                                        tGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGCTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGt                                                               >  2:1211051/1‑150 (MQ=255)
                                                                                                ggttTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACAt                                                       <  1:2367410/150‑1 (MQ=255)
                                                                                                ggttTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACAt                                                       <  1:4660617/150‑1 (MQ=255)
                                                                                                ggttTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACAt                                                       <  1:2581230/150‑1 (MQ=255)
                                                                                                ggttTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACAt                                                       <  2:5430296/150‑1 (MQ=255)
                                                                                                  tttGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGt                                                     <  2:480884/150‑1 (MQ=255)
                                                                                                   ttGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTc                                                    <  2:3784502/150‑1 (MQ=255)
                                                                                                    tGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCa                                                   <  1:652880/150‑1 (MQ=255)
                                                                                                    tGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCa                                                   <  1:4444920/150‑1 (MQ=255)
                                                                                                    tGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCa                                                   <  2:526679/150‑1 (MQ=255)
                                                                                                      cgcACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCaaa                                                 <  2:4461387/150‑1 (MQ=255)
                                                                                                        cacaGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCaaaaa                                               >  2:2821433/1‑150 (MQ=255)
                                                                                                         acaGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAAt                                              >  2:4632675/1‑150 (MQ=255)
                                                                                                           aGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATAc                                            >  1:552907/1‑150 (MQ=255)
                                                                                                           aGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATAc                                            >  1:3493051/1‑150 (MQ=255)
                                                                                                           aGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATAc                                            <  1:531534/150‑1 (MQ=255)
                                                                                                                  gCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTAGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGttt                                     <  2:3922037/150‑1 (MQ=255)
                                                                                                                   ccGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTCATGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATCCAGGGTTTa                                    >  1:4259433/1‑150 (MQ=255)
                                                                                                                   ccGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTa                                    >  1:1418612/1‑150 (MQ=255)
                                                                                                                    cGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTAt                                   <  1:1079568/150‑1 (MQ=255)
                                                                                                                       cTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTAtctc                                <  2:827251/150‑1 (MQ=255)
                                                                                                                       cTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTAtctc                                <  1:5186766/150‑1 (MQ=255)
                                                                                                                         cGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCtt                              >  1:3893996/1‑150 (MQ=255)
                                                                                                                         cGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCtt                              >  1:2734493/1‑150 (MQ=255)
                                                                                                                         cGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCtt                              >  1:781790/1‑150 (MQ=255)
                                                                                                                         cGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCtt                              >  1:1051736/1‑150 (MQ=255)
                                                                                                                          ggTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCttt                             >  1:4899525/1‑150 (MQ=255)
                                                                                                                             ttCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGa                          >  1:5144892/1‑150 (MQ=255)
                                                                                                                               cccccTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACg                        >  1:4424766/1‑150 (MQ=255)
                                                                                                                                  ccTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAAAACAGGGTTTATCTCTTTTGACGGTa                     >  2:1849171/1‑150 (MQ=255)
                                                                                                                                   ctctTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTaa                    <  2:3692984/150‑1 (MQ=255)
                                                                                                                                             tACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTc          <  1:3354534/150‑1 (MQ=255)
                                                                                                                                             gaCGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTc          <  1:3855675/149‑1 (MQ=255)
                                                                                                                                              aCGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCt         >  1:1625682/1‑150 (MQ=255)
                                                                                                                                                gATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTAt       >  1:323924/1‑150 (MQ=255)
                                                                                                                                                gATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTAt       >  2:2370659/1‑150 (MQ=255)
                                                                                                                                                  tGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTa     <  1:2929469/150‑1 (MQ=255)
                                                                                                                                                  tGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTa     >  1:22445/1‑150 (MQ=255)
                                                                                                                                                   gTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTAc    <  2:5075211/150‑1 (MQ=255)
                                                                                                                                                   gTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTAc    <  2:1306115/150‑1 (MQ=255)
                                                                                                                                                     cGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACct  <  1:1712146/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: