Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 958,417 A→T 12.1% N323Y (AAC→TAC)  rpsA → 30S ribosomal subunit protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273958,4170AT12.1% 1040.9 / 106.4 327N323Y (AAC→TAC) rpsA30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base A (124/163);  new base T (23/17);  total (147/180)
Fisher's exact test for biased strand distribution p-value = 9.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 332 total aligned reads displayed.

GTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAG  >  CP009273/958268‑958566
                                                                                                                                                     |                                                                                                                                                     
   aTCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACg                                                                                                                                                    <  2:3695822/150‑1 (MQ=255)
      cGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTg                                                                                                                                                 >  1:1309896/1‑150 (MQ=255)
      cGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTg                                                                                                                                                 >  2:2088778/1‑150 (MQ=255)
      cGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTg                                                                                                                                                 >  1:3007085/1‑150 (MQ=255)
         aaGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCg                                                                                                                                              <  2:3651903/150‑1 (MQ=255)
         aaGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCg                                                                                                                                              >  2:4611096/1‑150 (MQ=255)
               ccAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAg                                                                                                                                        <  2:2211030/150‑1 (MQ=255)
                cAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTTACGTTGGCGATGTAGt                                                                                                                                       >  1:3871714/1‑150 (MQ=255)
                cAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGt                                                                                                                                       <  2:1450773/150‑1 (MQ=255)
                 aaaCTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTg                                                                                                                                      >  2:3416199/1‑150 (MQ=255)
                   actgactgGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGa                                                                                                                                    >  2:3736073/1‑150 (MQ=255)
                    ctgactgGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGaa                                                                                                                                   <  2:4647315/150‑1 (MQ=255)
                        ctgGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTa                                                                                                                               <  1:2752250/150‑1 (MQ=255)
                         tgGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTAt                                                                                                                              <  2:4752237/150‑1 (MQ=255)
                          ggTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATg                                                                                                                             <  1:4898690/150‑1 (MQ=255)
                            tCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGt                                                                                                                           >  1:1955761/1‑150 (MQ=255)
                                 tGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTgg                                                                                                                      <  1:3257857/150‑1 (MQ=255)
                                  gACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGa                                                                                                                     <  2:5070587/150‑1 (MQ=255)
                                   aCCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGat                                                                                                                    >  1:1615741/1‑150 (MQ=255)
                                    ccAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTTACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGata                                                                                                                   <  1:1856653/150‑1 (MQ=255)
                                     cAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGatat                                                                                                                  <  1:5101149/150‑1 (MQ=255)
                                         cTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACCCGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcgac                                                                                                              <  2:4446736/150‑1 (MQ=255)
                                         cTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcgac                                                                                                              <  1:2131339/150‑1 (MQ=255)
                                          tGACCGACGACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcgacg                                                                                                             <  1:1491332/150‑1 (MQ=255)
                                                aCTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAc                                                                                                       <  1:4788167/150‑1 (MQ=255)
                                                aCTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAc                                                                                                       <  1:2243017/150‑1 (MQ=255)
                                                   aCGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtc                                                                                                    <  2:3342503/150‑1 (MQ=255)
                                                   aCGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtc                                                                                                    <  1:3504559/150‑1 (MQ=255)
                                                     ggCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtcgt                                                                                                  >  2:2992178/1‑150 (MQ=255)
                                                      gctgctTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtcgtc                                                                                                 <  2:3970451/150‑1 (MQ=255)
                                                          ctTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTAt                                                                                             <  2:3795019/150‑1 (MQ=255)
                                                          ctTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTAt                                                                                             <  2:4310357/150‑1 (MQ=255)
                                                          ctTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTAt                                                                                             <  2:437323/150‑1 (MQ=255)
                                                           ttCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTAtc                                                                                            <  1:938280/150‑1 (MQ=255)
                                                             cGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTAtctc                                                                                          >  2:2724981/1‑150 (MQ=255)
                                                              gTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTcc                                                                                         <  1:4505575/150‑1 (MQ=255)
                                                              gTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTcc                                                                                         >  2:3612334/1‑150 (MQ=255)
                                                               ttGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTccc                                                                                        >  1:4089148/1‑150 (MQ=255)
                                                               ttGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTccc                                                                                        <  2:3282269/150‑1 (MQ=255)
                                                                tGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCt                                                                                       <  1:3395678/150‑1 (MQ=255)
                                                                 gAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTg                                                                                      <  1:4374472/150‑1 (MQ=255)
                                                                   anTCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTggg                                                                                    >  1:4172727/1‑150 (MQ=255)
                                                                   aaTCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTggg                                                                                    >  2:5125118/1‑150 (MQ=255)
                                                                     tCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTc                                                                                  >  2:495351/1‑150 (MQ=255)
                                                                      cGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTAGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCt                                                                                 <  1:1151846/150‑1 (MQ=255)
                                                                        aagaagGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAACACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGa                                                                               <  1:3075565/150‑1 (MQ=255)
                                                                        aagaagGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGa                                                                               <  1:2595048/150‑1 (MQ=255)
                                                                        aagaagGCGTTGAAGGCATGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGAAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGAATCTCCCTGGGTCTGa                                                                               >  1:1004816/1‑150 (MQ=255)
                                                                         agaagGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGaa                                                                              >  2:1780091/1‑150 (MQ=255)
                                                                            agGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGACCGTCGTCGTATCTCCCTGGGTCTGAAACa                                                                           <  2:2181324/150‑1 (MQ=255)
                                                                              gCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGt                                                                         <  1:57519/150‑1 (MQ=255)
                                                                               cGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTTACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTg                                                                        >  1:668067/1‑150 (MQ=255)
                                                                                  tGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCaa                                                                     >  1:1330201/1‑150 (MQ=255)
                                                                                  tGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCaa                                                                     >  2:177007/1‑150 (MQ=255)
                                                                                   gAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTTACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCaaa                                                                    >  2:4173682/1‑150 (MQ=255)
                                                                                     aGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGc                                                                  <  2:3316264/150‑1 (MQ=255)
                                                                                     aGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGc                                                                  <  1:2407265/150‑1 (MQ=255)
                                                                                       gCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTa                                                                <  1:540315/150‑1 (MQ=255)
                                                                                        ccTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTaa                                                               <  1:3757981/150‑1 (MQ=255)
                                                                                         cTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTTACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAAc                                                              <  2:76099/150‑1 (MQ=255)
                                                                                          tGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAAcc                                                             >  1:2158435/1‑150 (MQ=255)
                                                                                           ggTACACGTTTCAGAAATGGACTGGAACAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAAGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAAccc                                                            <  1:3528288/150‑1 (MQ=255)
                                                                                            gTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAACCAGTGCAAAGCTAACCCg                                                           <  2:3596377/150‑1 (MQ=255)
                                                                                             tACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGt                                                          <  2:873833/150‑1 (MQ=255)
                                                                                             tACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCCGAAACAGTGCAAAGCTAACCCGt                                                          >  2:4809809/1‑150 (MQ=255)
                                                                                                  gTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGgcag                                                     >  2:4679582/1‑150 (MQ=255)
                                                                                                  gTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGgcag                                                     >  2:4981799/1‑150 (MQ=255)
                                                                                                    ttCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTAGGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGgcagca                                                   >  1:1839644/1‑150 (MQ=255)
                                                                                                     tCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGcagcag                                                  >  1:1367017/1‑150 (MQ=255)
                                                                                                       cGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGtt                                                <  1:1655219/150‑1 (MQ=255)
                                                                                                        gAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTc                                               >  1:3213590/1‑150 (MQ=255)
                                                                                                          aaTGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTcgc                                             >  2:4579463/1‑150 (MQ=255)
                                                                                                             ggactggacCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGa                                          >  1:734354/1‑150 (MQ=255)
                                                                                                              gactggacCAACAAAAACATCCACCCGTCCAAAGTTGTTTACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGaa                                         >  1:3334737/1‑150 (MQ=255)
                                                                                                                  ggacCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAAccc                                     <  1:4344345/150‑1 (MQ=255)
                                                                                                                  ggacCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGGTATGGTTCTGGATATCGACGCAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGGGCAAAGCGAACCCGTGGCAGCAGGTCGCGGAAcccc                                     >  2:3290206/1‑150 (MQ=255)
                                                                                                                    acCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCcac                                   <  2:4877052/150‑1 (MQ=255)
                                                                                                                     ccAACAAACACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCcaca                                  <  2:2313744/150‑1 (MQ=255)
                                                                                                                     ccAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAACCCcaca                                  >  2:2880110/1‑150 (MQ=255)
                                                                                                                      caacaaAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCacaa                                 >  2:615079/1‑150 (MQ=255)
                                                                                                                      caacaaAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCacaa                                 <  2:4473545/150‑1 (MQ=255)
                                                                                                                        acaaAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCacaaca                               <  2:168060/150‑1 (MQ=255)
                                                                                                                           aaaaCATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAgg                            <  2:107450/150‑1 (MQ=255)
                                                                                                                           aaaaCATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAgg                            <  1:1188684/150‑1 (MQ=255)
                                                                                                                            aaaCATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAggg                           <  2:1309862/150‑1 (MQ=255)
                                                                                                                                 tCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGAcc                      <  2:3463623/150‑1 (MQ=255)
                                                                                                                                 tCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGAcc                      <  2:1826978/150‑1 (MQ=255)
                                                                                                                                 tCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGAcc                      <  1:4244256/150‑1 (MQ=255)
                                                                                                                                  ccACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCg                     >  1:2941032/1‑150 (MQ=255)
                                                                                                                                   cACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCgt                    >  2:2027009/1‑150 (MQ=255)
                                                                                                                                    aCCCGTCCAAAGTTGTTTACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAATAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCgtg                   >  1:1020105/1‑150 (MQ=255)
                                                                                                                                     cccGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCgtgt                  >  2:3494465/1‑150 (MQ=255)
                                                                                                                                     cccGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCgtgt                  >  2:1648532/1‑150 (MQ=255)
                                                                                                                                      ccGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGtt                 >  2:3088176/1‑150 (MQ=255)
                                                                                                                                      ccGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGtt                 <  1:498663/150‑1 (MQ=255)
                                                                                                                                      ccGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGtt                 <  1:2459852/150‑1 (MQ=255)
                                                                                                                                        gTCCAAAGTTGTTTACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGa               <  1:508873/150‑1 (MQ=255)
                                                                                                                                           cAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAgg            <  1:4832046/150‑1 (MQ=255)
                                                                                                                                               gttgttAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTaaa        >  1:3645361/1‑150 (MQ=255)
                                                                                                                                                    tAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCaa   >  2:4239343/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAG  >  CP009273/958268‑958566

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: