Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,640,268 A→C 19.7% K39Q (AAG→CAG)  ydfV → Qin prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,640,2680AC19.7% 766.2 / 221.1 362K39Q (AAG→CAG) ydfVQin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (168/123);  new base C (41/30);  total (209/153)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

Only 100 of 367 total aligned reads displayed.

ATGCTACCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTG  >  CP009273/1640119‑1640417
                                                                                                                                                     |                                                                                                                                                     
aTGCTACCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTa                                                                                                                                                       >  1:2093760/1‑150 (MQ=255)
 tGCTACCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTaa                                                                                                                                                      <  2:1178438/150‑1 (MQ=255)
    tACCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTc                                                                                                                                                   >  2:5018977/1‑150 (MQ=255)
     aCCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCTCTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTcc                                                                                                                                                  <  1:3108310/150‑1 (MQ=255)
        cATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTTTGCAAGGAGCAGCCATTGCTAAGTCCATc                                                                                                                                               >  2:3164628/1‑150 (MQ=255)
        cATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATc                                                                                                                                               >  2:2288251/1‑150 (MQ=255)
         aTGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATcc                                                                                                                                              <  2:3077333/150‑1 (MQ=255)
          tGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCt                                                                                                                                             <  2:2434312/150‑1 (MQ=255)
          tGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCt                                                                                                                                             <  2:2022974/150‑1 (MQ=255)
          tGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCt                                                                                                                                             >  1:4191730/1‑150 (MQ=255)
           gTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTg                                                                                                                                            <  1:3455989/150‑1 (MQ=255)
              ttACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTAt                                                                                                                                         >  1:2409709/1‑150 (MQ=255)
                  acCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATtgtg                                                                                                                                     <  1:2453062/150‑1 (MQ=255)
                     tctTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAg                                                                                                                                  >  1:2249987/1‑150 (MQ=255)
                        tGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTCTTGTGCAGGTc                                                                                                                               >  2:851762/1‑150 (MQ=255)
                          tAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAg                                                                                                                             >  2:1778217/1‑150 (MQ=255)
                          tAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAg                                                                                                                             >  2:1777466/1‑150 (MQ=255)
                           aaTTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGc                                                                                                                            >  1:212484/1‑150 (MQ=255)
                           aaTTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGc                                                                                                                            >  1:1753111/1‑150 (MQ=255)
                             ttACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTc                                                                                                                          >  1:2827084/1‑150 (MQ=255)
                              tACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCg                                                                                                                         >  1:1445150/1‑150 (MQ=255)
                                  aaTGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGtttt                                                                                                                     <  2:4291522/150‑1 (MQ=255)
                                  aaTGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGAGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGtttt                                                                                                                     >  1:3636730/1‑150 (MQ=255)
                                   aTGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGttttt                                                                                                                    >  2:462621/1‑150 (MQ=255)
                                   aTGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGttttt                                                                                                                    >  1:3091115/1‑150 (MQ=255)
                                     gTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCTCGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTaa                                                                                                                  <  2:3027775/150‑1 (MQ=255)
                                       cATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTGTTTATag                                                                                                                >  2:3985950/1‑150 (MQ=255)
                                        aTTACAAGTATCGCACTTCAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAaga                                                                                                               <  2:4985829/150‑1 (MQ=255)
                                        aTTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAaga                                                                                                               <  1:1979183/150‑1 (MQ=255)
                                        aTTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAaga                                                                                                               >  2:1275009/1‑150 (MQ=255)
                                         ttACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAagtg                                                                                                              >  2:310570/1‑148 (MQ=255)
                                            cAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTcc                                                                                                           <  1:673283/150‑1 (MQ=255)
                                             taGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCg                                                                                                          <  1:4269004/149‑1 (MQ=255)
                                             aaGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTCTTGTGCAGGTCAGCTCGTTTTTATAGAGTCCg                                                                                                          >  2:4550056/1‑150 (MQ=255)
                                             aaGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCg                                                                                                          >  1:2580419/1‑150 (MQ=255)
                                             aaGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCg                                                                                                          <  2:4799878/150‑1 (MQ=255)
                                              aGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCgg                                                                                                         >  1:254990/1‑150 (MQ=255)
                                                tATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGcc                                                                                                       <  1:1737852/150‑1 (MQ=255)
                                                   cGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCcatc                                                                                                    >  2:4254323/1‑150 (MQ=255)
                                                     cACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCcatcat                                                                                                  >  1:4843317/1‑150 (MQ=255)
                                                          caacaTGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTAc                                                                                             <  2:1570895/150‑1 (MQ=255)
                                                           aacaTGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACt                                                                                            >  1:512816/1‑150 (MQ=255)
                                                           aacaTGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACt                                                                                            <  1:524358/150‑1 (MQ=255)
                                                            acaTGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTg                                                                                           <  1:3355871/150‑1 (MQ=255)
                                                              aTGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGt                                                                                         >  2:4837102/1‑150 (MQ=255)
                                                                gCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTAc                                                                                       >  1:4012958/1‑150 (MQ=255)
                                                                  ttAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCTTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAg                                                                                     <  2:1983434/150‑1 (MQ=255)
                                                                  ttAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTTCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAg                                                                                     >  2:1292396/1‑150 (MQ=255)
                                                                   tAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGa                                                                                    >  1:4935001/1‑150 (MQ=255)
                                                                     gggCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGaca                                                                                  <  2:1466428/150‑1 (MQ=255)
                                                                       gCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACAcc                                                                                <  1:984863/150‑1 (MQ=255)
                                                                        cAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCa                                                                               <  2:3706765/150‑1 (MQ=255)
                                                                         aaGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCat                                                                              >  2:233414/1‑150 (MQ=255)
                                                                           gTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCatat                                                                            >  1:2320985/1‑150 (MQ=255)
                                                                             cACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATAc                                                                          >  1:3919775/1‑150 (MQ=255)
                                                                              aCGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACt                                                                         >  1:4082882/1‑150 (MQ=255)
                                                                                gaaggaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACttt                                                                       >  1:3574819/1‑150 (MQ=255)
                                                                                gaaggaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACttt                                                                       >  1:1225351/1‑150 (MQ=255)
                                                                                gaaggaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACttt                                                                       >  1:3193349/1‑150 (MQ=255)
                                                                                    gaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGa                                                                   <  1:3821401/150‑1 (MQ=255)
                                                                                             aaGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGc                                                          >  2:1947824/1‑150 (MQ=255)
                                                                                              agtagtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCt                                                         >  1:1012687/1‑150 (MQ=255)
                                                                                               gtagtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTa                                                        >  1:1214954/1‑150 (MQ=255)
                                                                                                tagtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTAc                                                       >  1:1816096/1‑150 (MQ=255)
                                                                                                  gtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACat                                                     <  2:2609163/150‑1 (MQ=255)
                                                                                                    cGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCCGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATCTACTTTGTGACGGTCAGGCTACatat                                                   >  1:266716/1‑150 (MQ=255)
                                                                                                     gTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATg                                                  >  2:460167/1‑150 (MQ=255)
                                                                                                        aGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGcac                                               >  2:1633413/1‑150 (MQ=255)
                                                                                                           aCGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCAcaac                                            <  1:1375289/150‑1 (MQ=255)
                                                                                                           aCGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTACAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCAcaac                                            <  2:5096898/150‑1 (MQ=255)
                                                                                                             ggTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCAcaactc                                          >  2:4051160/1‑150 (MQ=255)
                                                                                                                 aTCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCAcaactcaact                                      >  1:1663117/1‑150 (MQ=255)
                                                                                                                  tCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACtt                                     >  1:2922910/1‑150 (MQ=255)
                                                                                                                       gTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCa                                >  1:834580/1‑150 (MQ=255)
                                                                                                                        tCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTGCATATGCACAACTCAACTTATTCAt                               <  1:4254323/150‑1 (MQ=255)
                                                                                                                         ccGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATc                              <  2:7498/150‑1 (MQ=255)
                                                                                                                         ccGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATc                              >  2:4406542/1‑150 (MQ=255)
                                                                                                                         ccGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATc                              >  1:1380872/1‑150 (MQ=255)
                                                                                                                               tgtgCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAtttt                        <  2:2232294/150‑1 (MQ=255)
                                                                                                                                   ccAGGATCATCGATTTATAAGTCAATCCTGTATTGTGCATGTCAGATCGTTTTTAAAGAGTCCTGCCATCATCTTTCTGGTACAGACTCAATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATATATTTTTTGc                    <  2:2785823/150‑1 (MQ=255)
                                                                                                                                     aGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCtt                  <  1:1919205/150‑1 (MQ=255)
                                                                                                                                      ggAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCttt                 <  2:3530074/150‑1 (MQ=255)
                                                                                                                                       gAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTa                <  2:3419564/150‑1 (MQ=255)
                                                                                                                                        agcagcCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAg               <  1:4102812/150‑1 (MQ=255)
                                                                                                                                        agcagcCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAg               <  2:1753169/150‑1 (MQ=255)
                                                                                                                                         gcagcCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGc              <  1:778149/150‑1 (MQ=255)
                                                                                                                                         gcagcCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGc              <  1:1391168/150‑1 (MQ=255)
                                                                                                                                         gcagcCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGc              <  2:4167304/150‑1 (MQ=255)
                                                                                                                                          cagcCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCa             <  1:15624/150‑1 (MQ=255)
                                                                                                                                           agcCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCAt            >  2:2578147/1‑150 (MQ=255)
                                                                                                                                             ccATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGt          >  1:4583614/1‑150 (MQ=255)
                                                                                                                                               aTTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCa        >  2:3380196/1‑150 (MQ=255)
                                                                                                                                                ttGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAg       >  2:3278746/1‑150 (MQ=255)
                                                                                                                                                ttGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAg       >  2:653358/1‑150 (MQ=255)
                                                                                                                                                ttGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAg       >  1:4368018/1‑150 (MQ=255)
                                                                                                                                                ttGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAg       >  2:3820312/1‑150 (MQ=255)
                                                                                                                                                 tGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAgt      >  2:4371478/1‑150 (MQ=255)
                                                                                                                                                  gCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAgtg     >  2:5216644/1‑150 (MQ=255)
                                                                                                                                                    tAAGTCCATCCTGTATTGTGCTGGTCAGATAGTTTTTAAAGAGTCCGTCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTAAATATGCACTACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGtt   <  2:2031102/150‑1 (MQ=255)
                                                                                                                                                     aaGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTg  <  2:5146336/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGCTACCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTG  >  CP009273/1640119‑1640417

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: