Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,924,395 C→A 5.4% D85Y (GAT→TAT)  yebF ← extracellular Colicin M immunity family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,924,3950CA 5.4% 882.6 / 15.9 230D85Y (GAT→TAT) yebFextracellular Colicin M immunity family protein
Reads supporting (aligned to +/- strand):  ref base C (114/103);  new base A (6/7);  total (120/110)
Fisher's exact test for biased strand distribution p-value = 7.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

CTCGCGCCGCATCCGCCAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTCCGGCAGCATCCAGATCTTCA  >  CP009273/1924247‑1924543
                                                                                                                                                    |                                                                                                                                                    
cTCGCGCCGCATCCGCCAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTatc                                                                                                                                                      <  2:3319512/149‑1 (MQ=255)
cTCGCGCCGCATCCGCCAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTatc                                                                                                                                                      >  1:2225593/1‑149 (MQ=255)
cTCGCGCCGCATCCGCCAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTatc                                                                                                                                                      <  2:1269217/149‑1 (MQ=255)
   gcgcCGCATCCGCCAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTatcatc                                                                                                                                                   <  1:3647165/149‑1 (MQ=255)
       cgcATCCGCCAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTa                                                                                                                                               >  1:40093/1‑149 (MQ=255)
       cgcATCCGCCAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTa                                                                                                                                               <  2:1641575/149‑1 (MQ=255)
        gcATCCGCCAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTAc                                                                                                                                              <  2:3243796/149‑1 (MQ=255)
          aTCCGCCAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTAccc                                                                                                                                            <  2:900968/149‑1 (MQ=255)
             cGCCAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGa                                                                                                                                         <  1:2900041/149‑1 (MQ=255)
               ccAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAAt                                                                                                                                       <  2:2290569/149‑1 (MQ=255)
                cAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATg                                                                                                                                      >  1:2040324/1‑149 (MQ=255)
                 aaCAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGt                                                                                                                                     <  1:2011939/149‑1 (MQ=255)
                    agacAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACACTGAATGACCt                                                                                                                                  >  2:868798/4‑149 (MQ=255)
                          tCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAAc                                                                                                                            >  2:2328574/1‑149 (MQ=255)
                          tCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAAc                                                                                                                            <  1:3277773/149‑1 (MQ=255)
                          tCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAAc                                                                                                                            >  2:2328593/1‑149 (MQ=255)
                            gAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACATGGTAATGAATATCGGCACTTTTACCACGCACCGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCacacag                                                                                                                          >  1:2792516/1‑147 (MQ=255)
                             aaTTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGa                                                                                                                         <  1:4102719/149‑1 (MQ=255)
                             aaTTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGa                                                                                                                         >  2:3434607/1‑149 (MQ=255)
                               ttATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGAcc                                                                                                                       >  1:1922866/1‑149 (MQ=255)
                                  tCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACACTGAATGTCCTGCAAACTGACCCaa                                                                                                                    >  2:1667333/1‑149 (MQ=255)
                                    aGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGc                                                                                                                  >  1:1240709/1‑149 (MQ=255)
                                    aGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGc                                                                                                                  >  1:1239202/1‑149 (MQ=255)
                                     gAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGcc                                                                                                                 <  2:987979/149‑1 (MQ=255)
                                      aaGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCa                                                                                                                >  2:2858035/1‑149 (MQ=255)
                                      aaGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCa                                                                                                                >  1:3024418/1‑149 (MQ=255)
                                      aaGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCa                                                                                                                <  2:995326/149‑1 (MQ=255)
                                      aaGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCa                                                                                                                <  1:2778032/149‑1 (MQ=255)
                                      aaGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCa                                                                                                                >  2:518768/1‑149 (MQ=255)
                                          cATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACggg                                                                                                            >  2:2207732/1‑149 (MQ=255)
                                          cATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACggg                                                                                                            <  2:71620/149‑1 (MQ=255)
                                          cATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACggg                                                                                                            >  1:1845116/1‑149 (MQ=255)
                                            ttAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGAt                                                                                                          <  1:148262/149‑1 (MQ=255)
                                            ttAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGAt                                                                                                          <  1:2906070/149‑1 (MQ=255)
                                             tAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATc                                                                                                         <  2:4095252/149‑1 (MQ=255)
                                               aCGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATAATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCgg                                                                                                       >  1:2198363/1‑149 (MQ=255)
                                                 gccgcTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGcc                                                                                                     <  1:3858198/149‑1 (MQ=255)
                                                   cgcTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTg                                                                                                   >  2:3549257/1‑149 (MQ=255)
                                                    gcTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGa                                                                                                  >  1:907993/1‑149 (MQ=255)
                                                       gATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCg                                                                                               >  1:1254842/1‑149 (MQ=255)
                                                        atatTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGa                                                                                              >  2:638894/1‑149 (MQ=255)
                                                        atatTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGa                                                                                              >  2:571644/1‑149 (MQ=255)
                                                         tatTCCGCCAGTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATAATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGAc                                                                                             <  2:3822506/149‑1 (MQ=255)
                                                             ccGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAAtt                                                                                         >  2:1550460/1‑149 (MQ=255)
                                                             ccGCCAGTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTATTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAAtt                                                                                         <  2:2216354/149‑1 (MQ=255)
                                                              cGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAAttt                                                                                        <  2:2507880/149‑1 (MQ=255)
                                                                ccATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAAttttt                                                                                      <  1:4057738/149‑1 (MQ=255)
                                                                 cATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATAATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAAtttttt                                                                                     <  1:818915/149‑1 (MQ=255)
                                                                  aTTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTg                                                                                    >  1:220371/1‑149 (MQ=255)
                                                                   ttCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGa                                                                                   >  1:2270721/1‑149 (MQ=255)
                                                                    tCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGat                                                                                  <  2:1724739/149‑1 (MQ=255)
                                                                       cGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGatcat                                                                               <  1:2307955/149‑1 (MQ=255)
                                                                        gCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATAATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGatcatc                                                                              <  2:2275281/149‑1 (MQ=255)
                                                                           ttGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGc                                                                           >  1:2620519/1‑149 (MQ=255)
                                                                           ttGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGc                                                                           >  2:2529134/1‑149 (MQ=255)
                                                                                gTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAAc                                                                      >  2:2741613/1‑149 (MQ=255)
                                                                                gTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAAc                                                                      <  2:253031/149‑1 (MQ=255)
                                                                                      gCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCa                                                                <  1:1551560/149‑1 (MQ=255)
                                                                                       cTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCAc                                                               <  2:3903156/149‑1 (MQ=255)
                                                                                        tGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCAcg                                                              >  2:3667161/1‑149 (MQ=255)
                                                                                           ccTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGa                                                           <  1:3067400/149‑1 (MQ=255)
                                                                                             tGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGAtt                                                         >  1:489901/1‑149 (MQ=255)
                                                                                               gTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGAtttt                                                       <  1:1654787/149‑1 (MQ=255)
                                                                                                 aaTGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTttgt                                                     >  1:981441/1‑149 (MQ=255)
                                                                                                     aaTATCGTCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGAt                                                 <  1:2892131/149‑1 (MQ=255)
                                                                                                     aaTATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGAt                                                 <  2:1324929/149‑1 (MQ=255)
                                                                                                      atatCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATa                                                >  1:2826010/1‑149 (MQ=255)
                                                                                                        atCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAAt                                              <  1:1741804/149‑1 (MQ=255)
                                                                                                          cGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCa                                            <  2:1800657/149‑1 (MQ=255)
                                                                                                          cGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCa                                            <  1:2477508/149‑1 (MQ=255)
                                                                                                              aCTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGtt                                        <  2:1922866/149‑1 (MQ=255)
                                                                                                              aCTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGtt                                        >  2:3172176/1‑149 (MQ=255)
                                                                                                              aCTTTGACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAACCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGtt                                        <  2:3218477/149‑1 (MQ=255)
                                                                                                               cTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGttt                                       >  1:1194718/1‑149 (MQ=255)
                                                                                                               cTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGttt                                       >  1:2850899/1‑149 (MQ=255)
                                                                                                                ttttACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGtttt                                      <  1:3740327/149‑1 (MQ=255)
                                                                                                                ttttACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGtttt                                      >  2:4071434/1‑149 (MQ=255)
                                                                                                                   tACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACg                                   >  1:1400560/1‑149 (MQ=255)
                                                                                                                    aCCACGCACGGTTAGCGGTACTGACCATTTATAATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGc                                  <  1:1142134/149‑1 (MQ=255)
                                                                                                                     ccACTCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCt                                 >  1:559212/1‑149 (MQ=255)
                                                                                                                            cgGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCa                          <  2:2554523/149‑1 (MQ=255)
                                                                                                                              gTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAAt                        >  1:3823526/1‑149 (MQ=255)
                                                                                                                              gTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAAt                        >  1:3692944/1‑149 (MQ=255)
                                                                                                                               ttAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAAtt                       >  1:3306277/1‑149 (MQ=255)
                                                                                                                                tAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTc                      <  2:3657249/149‑1 (MQ=255)
                                                                                                                                tAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTc                      >  1:3945101/1‑149 (MQ=255)
                                                                                                                                tAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTc                      >  1:3945106/1‑149 (MQ=255)
                                                                                                                                tAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTc                      <  1:1564965/149‑1 (MQ=255)
                                                                                                                                 aGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTcc                     <  1:2909337/149‑1 (MQ=255)
                                                                                                                                 aGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTcc                     >  2:486746/1‑149 (MQ=255)
                                                                                                                                 aGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTcc                     >  1:2049836/1‑149 (MQ=255)
                                                                                                                                        cTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTCCGgcagca              >  1:111172/1‑149 (MQ=255)
                                                                                                                                        cTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTCCGgcagca              >  1:110193/1‑149 (MQ=255)
                                                                                                                                          gACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTCCGGCAGCATc            <  2:76578/149‑1 (MQ=255)
                                                                                                                                             cATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTCCGGCAGCATCCAg         <  1:3751062/149‑1 (MQ=255)
                                                                                                                                               tttATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTCCGGCAGCATCCAGAt       >  2:707817/1‑149 (MQ=255)
                                                                                                                                               tttATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTCCGGCAGCATCCAGAt       >  2:3661757/1‑149 (MQ=255)
                                                                                                                                                  atcatcTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTCCGGCAGCATCCAGATCtt    <  2:2620519/149‑1 (MQ=255)
                                                                                                                                                  atcatcTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTCCGGCAGCATCCAGATCtn    <  1:3519496/149‑2 (MQ=255)
                                                                                                                                                   tcatcTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTCCGGCAGCATCCAGATCTTc   >  1:2450917/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCGCGCCGCATCCGCCAACAGTCAATCGAATTATCAGAAGCCATTAACGCCGCTGATATTCCGCCATTCCCGCTTTGCAGTCCACGCTGACCTGGTAATGAATATCGGCACTTTTACCACGCACGGTTAGCGGTACTGACCATTTATCATCTTTACCCTGAATGTCCTGCAAACTGACCCAAGCCACGGGATCGGCCTGACCGACAATTTTTTGATCATCTGCCCAACGCGCCACGCGATTTTGTTGATAATCACGTTTTACGCTCGCGGCAATTCCGGCAGCATCCAGATCTTCA  >  CP009273/1924247‑1924543

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: