Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 520,254 G→T 9.9% R513L (CGC→CTC)  rhsD → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273520,2540GT 9.9% 1005.1 / 55.8 272R513L (CGC→CTC) rhsDRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base G (138/107);  new base T (15/12);  total (153/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 296 total aligned reads displayed.

GGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCC  >  CP009273/520105‑520403
                                                                                                                                                     |                                                                                                                                                     
  aCGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTAcgctaccgct                                                                                                                                                     >  2:1268679/1‑150 (MQ=255)
   cgggcgggAGACGAAATTGTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTAcgctaccgcta                                                                                                                                                    <  2:1032250/150‑1 (MQ=255)
    gggcgggAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTAcgctaccgctac                                                                                                                                                   >  1:2720160/1‑150 (MQ=255)
    gggcgggAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTAcgctaccgctac                                                                                                                                                   >  1:3553737/1‑150 (MQ=255)
     ggcgggAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACg                                                                                                                                                  >  2:3114299/1‑150 (MQ=255)
       cgggAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGAt                                                                                                                                                >  2:568579/1‑150 (MQ=255)
       cgggAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGAt                                                                                                                                                >  2:3980063/1‑150 (MQ=255)
       cgggAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGAt                                                                                                                                                <  1:4083493/150‑1 (MQ=255)
          gagaCGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGAc                                                                                                                                             >  1:2567244/1‑150 (MQ=255)
           agaCGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGAcg                                                                                                                                            >  2:1107730/1‑150 (MQ=255)
           agaCGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGAcg                                                                                                                                            >  2:190508/1‑150 (MQ=255)
           agaCGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGAcg                                                                                                                                            >  2:4151401/1‑150 (MQ=255)
             aCGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGAcgcg                                                                                                                                          >  2:1075490/1‑150 (MQ=255)
             aCGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGCAACCAGAGCCGCCCGCAATAAGATGAACACGGCAGGCAGGTCTCGGAGACATCGCTCACAGGCGAGACAGTACGCTACCGCTACGATGActcg                                                                                                                                          >  2:3000724/1‑150 (MQ=255)
               gAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGca                                                                                                                                        <  1:1488890/150‑1 (MQ=255)
                aaaTTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGcac                                                                                                                                       <  1:4491961/150‑1 (MQ=255)
                aaaTTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCAGGGACTATGATGAACCGGGCAGGCTGGTATCGGATACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGcac                                                                                                                                       >  1:3204216/1‑150 (MQ=255)
                  aTTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAg                                                                                                                                     >  1:4565871/1‑150 (MQ=255)
                  aTTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAg                                                                                                                                     <  1:106610/150‑1 (MQ=255)
                      tACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAg                                                                                                                                 <  1:3115056/150‑1 (MQ=255)
                       aCTATACCGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGt                                                                                                                                <  1:1061810/150‑1 (MQ=255)
                        cTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGtt                                                                                                                               <  1:2950959/150‑1 (MQ=255)
                         tataACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTa                                                                                                                              <  1:3472569/150‑1 (MQ=255)
                         tataACGACGGGAACCAGCTGACGGATGTGGTGTCCCCGGACTGGCAGGATAGACGCCGGGAATATGATCAACCGTGCAGGCAGGTCTCGGAGACATCGCGCAGCGGGGAGACACTACGCTACCGCTTCGATGACGCGCACAGTGAGTTa                                                                                                                              >  2:1790389/1‑150 (MQ=255)
                          ataACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCTCTACGATGACGCGCACAGTGAGTTAc                                                                                                                             <  1:1283338/150‑1 (MQ=255)
                             acgacgGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCgg                                                                                                                          <  1:3826596/150‑1 (MQ=255)
                              cgacgGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGc                                                                                                                         >  2:781625/1‑150 (MQ=255)
                              cgacgGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGc                                                                                                                         >  2:3933427/1‑150 (MQ=255)
                               gacgGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGcg                                                                                                                        >  1:2531391/1‑150 (MQ=255)
                               gacgGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGcg                                                                                                                        >  2:3742871/1‑150 (MQ=255)
                               gacgGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGcg                                                                                                                        >  2:3742962/1‑150 (MQ=255)
                               gacgGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGcg                                                                                                                        >  1:3078429/1‑150 (MQ=255)
                                 cgGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGcgac                                                                                                                      <  1:3321449/150‑1 (MQ=255)
                                  gggAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGcgacg                                                                                                                     <  1:2831986/150‑1 (MQ=255)
                                  gggAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGcgacg                                                                                                                     >  1:777521/1‑150 (MQ=255)
                                    gAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCgacgac                                                                                                                   >  1:2038254/1‑150 (MQ=255)
                                     aaCCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACa                                                                                                                  <  1:3656881/150‑1 (MQ=255)
                                       ccAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAAc                                                                                                                >  1:140760/1‑150 (MQ=255)
                                        cAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCTCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACg                                                                                                               <  2:2693221/150‑1 (MQ=255)
                                        cAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACg                                                                                                               >  1:3120567/1‑150 (MQ=255)
                                            tGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATg                                                                                                           <  2:96974/150‑1 (MQ=255)
                                            tGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATg                                                                                                           <  1:398516/150‑1 (MQ=255)
                                              aCGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCTCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCg                                                                                                         >  2:3129860/1‑150 (MQ=255)
                                              aCGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCTCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCg                                                                                                         <  1:2342029/150‑1 (MQ=255)
                                              aCGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCg                                                                                                         <  1:1328533/150‑1 (MQ=255)
                                               cggcggTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGa                                                                                                        <  2:694591/150‑1 (MQ=255)
                                                    gtggtgTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGgc                                                                                                   >  1:2230964/1‑150 (MQ=255)
                                                      ggtgTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGgcag                                                                                                 >  1:2883727/1‑150 (MQ=255)
                                                        tgtCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCTCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGgcagca                                                                                               <  1:3425129/150‑1 (MQ=255)
                                                           ccccGGACGGGCTGGAGAGCCGCCGGGACTATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCTCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCAccc                                                                                            >  2:1569987/1‑150 (MQ=255)
                                                             ccGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCgg                                                                                          >  2:3767985/1‑150 (MQ=255)
                                                             ccGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCgg                                                                                          >  2:2563039/1‑150 (MQ=255)
                                                              cGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGc                                                                                         >  2:3244053/1‑150 (MQ=255)
                                                                  cGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGAt                                                                                     >  1:3779637/1‑150 (MQ=255)
                                                                     gCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGAc                                                                                  <  1:3303271/150‑1 (MQ=255)
                                                                       tGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCTCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCt                                                                                >  2:1771890/1‑150 (MQ=255)
                                                                        ggAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTg                                                                               >  2:1691656/1‑150 (MQ=255)
                                                                         gagaGCCGCCGGGAATATGATGCACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTgg                                                                              <  2:117571/150‑1 (MQ=255)
                                                                         gagaGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTgg                                                                              <  2:2429122/150‑1 (MQ=255)
                                                                         gagaGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTgg                                                                              >  2:2855732/1‑150 (MQ=255)
                                                                          agagCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGa                                                                             <  2:1490425/150‑1 (MQ=255)
                                                                           gagCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAg                                                                            <  2:2920740/150‑1 (MQ=255)
                                                                                  cgGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTAc                                                                     <  2:1886276/150‑1 (MQ=255)
                                                                                  cgGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTAc                                                                     >  2:143370/1‑150 (MQ=255)
                                                                                   gggAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACg                                                                    >  1:2633219/1‑150 (MQ=255)
                                                                                   gggAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACg                                                                    >  2:1014012/1‑150 (MQ=255)
                                                                                    ggAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACgg                                                                   >  2:3435298/1‑150 (MQ=255)
                                                                                     gAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCTCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACggg                                                                  >  1:1876142/1‑150 (MQ=255)
                                                                                       atatGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCa                                                                <  1:2206012/150‑1 (MQ=255)
                                                                                       atatGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCa                                                                <  1:1760125/150‑1 (MQ=255)
                                                                                          tgatgaACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGtt                                                             <  1:4182563/150‑1 (MQ=255)
                                                                                          tgatgaACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGtt                                                             <  1:67943/150‑1 (MQ=255)
                                                                                          tgatgaACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGtt                                                             >  2:2320889/1‑150 (MQ=255)
                                                                                              gaACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTg                                                         <  2:3539907/150‑1 (MQ=255)
                                                                                                aCCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGc                                                       >  2:4401982/1‑150 (MQ=255)
                                                                                                aCCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGc                                                       >  2:3836262/1‑150 (MQ=255)
                                                                                                      cAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCAc                                                 <  2:2650543/150‑1 (MQ=255)
                                                                                                       aGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGTGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCAcc                                                >  1:4044058/1‑150 (MQ=255)
                                                                                                          cTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGAc                                             >  1:2407347/1‑150 (MQ=255)
                                                                                                               aTCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTc                                        <  2:1314673/150‑1 (MQ=255)
                                                                                                                 cGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCgg                                      >  2:2235433/1‑150 (MQ=255)
                                                                                                                     gaCATCGCGCCGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTa                                  <  2:2494184/150‑1 (MQ=255)
                                                                                                                       cATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTAcc                                >  2:4580495/1‑150 (MQ=255)
                                                                                                                           gcgcAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACATTGAGTTACAGGAGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGaa                            >  1:1556158/1‑149 (MQ=255)
                                                                                                                           gcgcAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGAc                            <  2:2961719/150‑1 (MQ=255)
                                                                                                                            cgcAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGAcc                           <  2:2057768/150‑1 (MQ=255)
                                                                                                                             gcAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGAccc                          >  1:1168428/1‑150 (MQ=255)
                                                                                                                             gcAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGAccc                          >  2:2741548/1‑150 (MQ=255)
                                                                                                                                 cGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTa                      >  2:509618/1‑150 (MQ=255)
                                                                                                                                  ggggAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTAt                     >  1:1371376/1‑150 (MQ=255)
                                                                                                                                   gggAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATg                    >  2:663249/1‑150 (MQ=255)
                                                                                                                                      agaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAAt                 >  1:3938441/1‑150 (MQ=255)
                                                                                                                                       gaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATa                >  2:4265664/1‑150 (MQ=255)
                                                                                                                                       gaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATa                >  2:2692710/1‑150 (MQ=255)
                                                                                                                                         cAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACg              <  2:1143860/150‑1 (MQ=255)
                                                                                                                                         cAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACg              <  2:2107130/150‑1 (MQ=255)
                                                                                                                                             aCGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGACTACGACCg          >  1:1724279/1‑150 (MQ=255)
                                                                                                                                             aCGCTACCGCTACGATGACGCGCACAGGGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCg          >  2:2164462/1‑150 (MQ=255)
                                                                                                                                               gctaccgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCt        <  1:4303190/150‑1 (MQ=255)
                                                                                                                                                ctaccgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCtt       <  1:3807760/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCC  >  CP009273/520105‑520403

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: