Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 159,121 A→T 10.6% E177V (GAA→GTA)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273159,1210AT10.6% 1028.7 / 89.0 329E177V (GAA→GTA) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (122/172);  new base T (13/22);  total (135/194)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 335 total aligned reads displayed.

CGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTATTTTTACCTGGCGTCGGAGAAATTCAGC  >  CP009273/158973‑159270
                                                                                                                                                    |                                                                                                                                                     
cgcAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGaa                                                                                                                                                      <  2:4306876/150‑1 (MQ=255)
    gCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGgc                                                                                                                                                  <  2:608366/150‑1 (MQ=255)
      ttGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGgcgc                                                                                                                                                >  2:1241743/1‑150 (MQ=255)
       tGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCt                                                                                                                                               <  2:1374937/150‑1 (MQ=255)
        gCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTc                                                                                                                                              <  1:3608231/150‑1 (MQ=255)
            gCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGttt                                                                                                                                          <  2:3900483/150‑1 (MQ=255)
               gATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGTAGGGCGCTCGTTTCCg                                                                                                                                       <  2:2707382/150‑1 (MQ=255)
                aTTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCgg                                                                                                                                      <  1:956696/150‑1 (MQ=255)
                 tttGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGt                                                                                                                                     <  1:3256398/150‑1 (MQ=255)
                  ttggcgttggcgCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGtt                                                                                                                                    <  2:650860/150‑1 (MQ=255)
                   tggcgttggcgCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTg                                                                                                                                   >  2:1783623/1‑150 (MQ=255)
                      cgttggcgCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAAc                                                                                                                                >  1:1781368/1‑150 (MQ=255)
                       gttggcgCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACg                                                                                                                               <  1:3131384/150‑1 (MQ=255)
                         tggcgCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGcc                                                                                                                             <  1:1312100/150‑1 (MQ=255)
                           gcgcTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGt                                                                                                                           >  2:4412716/1‑150 (MQ=255)
                           gcgcTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGt                                                                                                                           >  1:3676317/1‑150 (MQ=255)
                              cTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTAt                                                                                                                        >  1:1063515/1‑150 (MQ=255)
                              cTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTAt                                                                                                                        >  2:2239820/1‑150 (MQ=255)
                              cTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTAt                                                                                                                        >  2:3363881/1‑150 (MQ=255)
                              cTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTAt                                                                                                                        <  2:399490/150‑1 (MQ=255)
                                   aCTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTAcc                                                                                                                   <  2:2422106/150‑1 (MQ=255)
                                   aCTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTAcc                                                                                                                   >  1:4484310/1‑150 (MQ=255)
                                      cGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCt                                                                                                                >  1:2250037/1‑150 (MQ=255)
                                          gtgCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGccc                                                                                                            >  2:3204750/1‑150 (MQ=255)
                                          gtgCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGccc                                                                                                            <  1:3855250/150‑1 (MQ=255)
                                           tgCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGTAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCcg                                                                                                           <  1:2345120/150‑1 (MQ=255)
                                             caacaaGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCcgcg                                                                                                         >  2:3030666/1‑150 (MQ=255)
                                             caacaaGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCcgcg                                                                                                         <  1:401869/150‑1 (MQ=255)
                                             caacaaGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCcgcg                                                                                                         <  1:1204662/150‑1 (MQ=255)
                                                    gTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGc                                                                                                  >  1:2915392/1‑150 (MQ=255)
                                                      cTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGt                                                                                                <  1:1361359/150‑1 (MQ=255)
                                                        gCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGttt                                                                                              <  2:2580656/150‑1 (MQ=255)
                                                          gTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTg                                                                                            >  1:3547096/1‑150 (MQ=255)
                                                          gTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTg                                                                                            >  2:1604560/1‑150 (MQ=255)
                                                           tgatgaCCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGa                                                                                           <  1:1486400/150‑1 (MQ=255)
                                                            gatgaCCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGAc                                                                                          <  1:3766814/150‑1 (MQ=255)
                                                            gatgaCCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGAc                                                                                          >  1:1299361/1‑150 (MQ=255)
                                                              tgaCCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGTAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGa                                                                                        >  2:1109727/1‑150 (MQ=255)
                                                              tgaCCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGa                                                                                        <  2:3086667/150‑1 (MQ=255)
                                                                aCCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATg                                                                                      >  2:3709700/1‑150 (MQ=255)
                                                                aCCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATg                                                                                      <  2:648659/150‑1 (MQ=255)
                                                                 ccTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGc                                                                                     >  1:1559759/1‑150 (MQ=255)
                                                                 ccTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGc                                                                                     <  2:169111/150‑1 (MQ=255)
                                                                  cTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGcc                                                                                    <  2:2003119/150‑1 (MQ=255)
                                                                  cTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGcc                                                                                    >  1:1925105/1‑150 (MQ=255)
                                                                   ttAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGTAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCg                                                                                   >  1:2563178/1‑150 (MQ=255)
                                                                       aCTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGACGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGACCGCCGTTATGTACCGCTGCCCGCGCATCAGCGTTTTGTCGATGCCGTTGc                                                                               >  2:3100100/1‑150 (MQ=255)
                                                                       aCTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGc                                                                               <  2:850079/150‑1 (MQ=255)
                                                                       aCTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGc                                                                               >  2:4153888/1‑150 (MQ=255)
                                                                        ctgctgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGTAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCg                                                                              >  1:4225901/1‑150 (MQ=255)
                                                                        ctgctgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCg                                                                              <  1:4044707/150‑1 (MQ=255)
                                                                         tgctgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCgg                                                                             >  2:4067242/1‑150 (MQ=255)
                                                                         tgctgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCgg                                                                             >  1:2463708/1‑150 (MQ=255)
                                                                          gctgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAATCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGt                                                                            <  1:298919/150‑1 (MQ=255)
                                                                          gctgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGt                                                                            <  2:532527/150‑1 (MQ=255)
                                                                           ctgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCTCATCAGCGTTTTGACGATGCCGTTGCGGTa                                                                           <  1:4285999/150‑1 (MQ=255)
                                                                           ctgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTa                                                                           >  2:2847344/1‑150 (MQ=255)
                                                                           ctgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTa                                                                           >  2:2323179/1‑150 (MQ=255)
                                                                            tgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGTGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAg                                                                          >  1:2159699/1‑150 (MQ=255)
                                                                             gATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGc                                                                         >  2:3249836/1‑150 (MQ=255)
                                                                              aTTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGcc                                                                        <  1:2094764/150‑1 (MQ=255)
                                                                               ttATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGACGCCGTTGCGGTAGCCa                                                                       <  1:550138/150‑1 (MQ=255)
                                                                                tATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCAc                                                                      >  1:1282468/1‑150 (MQ=255)
                                                                                    tCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCt                                                                  <  2:4611498/150‑1 (MQ=255)
                                                                                     cGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGTAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTg                                                                 >  1:1462096/1‑150 (MQ=255)
                                                                                     cGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTg                                                                 <  1:3439856/150‑1 (MQ=255)
                                                                                      ggCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGa                                                                <  1:90603/150‑1 (MQ=255)
                                                                                      ggCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGa                                                                <  1:4544843/150‑1 (MQ=255)
                                                                                            ccTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGTAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAAtgct                                                          <  2:1370850/150‑1 (MQ=255)
                                                                                            ccTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCACTGAAAtgct                                                          >  1:651989/1‑150 (MQ=255)
                                                                                                gACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGt                                                      >  2:3584155/1‑150 (MQ=255)
                                                                                                  cAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCa                                                    <  2:3466283/150‑1 (MQ=255)
                                                                                                   aaCGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAg                                                   <  1:3366855/150‑1 (MQ=255)
                                                                                                     cGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGa                                                 >  2:378703/1‑150 (MQ=255)
                                                                                                      gACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGaa                                                >  2:1581789/1‑150 (MQ=255)
                                                                                                      gACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGaa                                                >  2:685376/1‑150 (MQ=255)
                                                                                                        ccGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAg                                              <  2:4010217/150‑1 (MQ=255)
                                                                                                         cGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGTAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGc                                             <  1:1998123/150‑1 (MQ=255)
                                                                                                         cGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGc                                             <  2:2477057/150‑1 (MQ=255)
                                                                                                          gCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCg                                            <  1:181223/150‑1 (MQ=255)
                                                                                                           cTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCgg                                           >  1:994406/1‑150 (MQ=255)
                                                                                                              gcagcaAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATc                                        <  2:2719736/150‑1 (MQ=255)
                                                                                                              gcagcaAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATc                                        <  1:3894233/150‑1 (MQ=255)
                                                                                                                 gcaAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCAtt                                     <  2:2815847/150‑1 (MQ=255)
                                                                                                                   aaaTGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGAGGATCATTAc                                   >  2:1840466/1‑150 (MQ=255)
                                                                                                                    aaTGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACt                                  <  1:3957501/150‑1 (MQ=255)
                                                                                                                    aaTGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACt                                  <  2:4474064/150‑1 (MQ=255)
                                                                                                                         tgcCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTAt                             <  1:4369421/150‑1 (MQ=255)
                                                                                                                           ccAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTAttt                           <  2:3648414/150‑1 (MQ=255)
                                                                                                                            cAGAAGCGCCTGTCGTCATCTCAGTAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTAtttt                          >  2:1949840/1‑150 (MQ=255)
                                                                                                                            cAGAAGCGCCTGTCGTCATCTCAGTAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTAtttt                          <  2:2394415/150‑1 (MQ=255)
                                                                                                                            cAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTAtttt                          >  1:826618/1‑150 (MQ=255)
                                                                                                                                 gcgcCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTATTTTTACCt                     <  2:1436548/150‑1 (MQ=255)
                                                                                                                                  cgcCTGTCGTCATCTCAGTAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTATTTTTACCTg                    <  1:2041998/150‑1 (MQ=255)
                                                                                                                                         cgtcATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTATTTTTACCTGGCGTCgg             <  2:1028377/150‑1 (MQ=255)
                                                                                                                                          gtcATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTATTTTTACCTGGCGTCGga            <  2:2299098/150‑1 (MQ=255)
                                                                                                                                          gtcATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTATTTTTACCTGGCGTCGga            >  2:940443/1‑150 (MQ=255)
                                                                                                                                             aTCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTATTTTTACCTGGCGTCGGAGaa         >  2:4228926/1‑150 (MQ=255)
                                                                                                                                                tcAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTATTTTTACCTGGCGTCGGAGAAAtt      >  2:564435/1‑150 (MQ=255)
                                                                                                                                                   gAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTATTTTTACCTGGCGTCGGAGAAATTCAg   <  2:1507065/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTACTGTTATTTTTACCTGGCGTCGGAGAAATTCAGC  >  CP009273/158973‑159270

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: